Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,581 C→T V65V (GTC→GTT araC → DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,5810CT96.6% 261.7 / ‑3.5 90V65V (GTC→GTTaraCDNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (52/34);  total (54/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGC  >  NC_000913/70335‑70822
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGTTTGAGGCCAACGGTTATTTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCTGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                                                                                                                                                                                                                                                    >  1:224758/1‑251 (MQ=255)
         tGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGGAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCggggg                                                                                                                                                                                                                                           >  1:241931/1‑249 (MQ=255)
          agACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAAAGCCGATTTAGGCCCACGGGTTTTTCGGTTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTTTCAATCTCACCATTCGCGGGCAGGGGGTGGGGAAAAATCAGGGACGAGAATTTTTTTGCCGACCgggggg                                                                                                                                                                                                                                          >  1:45825/2‑248 (MQ=255)
            aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata                                                                                                                                                                                                                                       <  2:945614/251‑1 (MQ=255)
                   gTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGCAAGGGTATATTCTCAATCTCCCCATTCGCGGGCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGAAATTTTGCt                                                                                                                                                                                                                               >  1:976612/1‑251 (MQ=255)
                      tcttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGTGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGt                                                                                                                                                                                                                             <  1:749800/250‑1 (MQ=255)
                                  tgcgggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGCTTACCGCCGATTGAGGCCAACGGTTATCTCGCTTTTTTTATCGACCGACCGCTGGGAATGCAAGGTTTTATTCTAGAACTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCCGGGACAGGAATTTTTTTGCCGACCGGGTGCTATTTTGCTGTTTCAGCCAgg                                                                                                                                                                                                                  >  1:446667/4‑248 (MQ=255)
                                  ggcgggAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGCCCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCCCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAgg                                                                                                                                                                                                                  >  1:271110/4‑248 (MQ=255)
                                          tATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATtcatc                                                                                                                                                                                                       >  1:757277/1‑251 (MQ=255)
                                                      ggCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCCACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGATATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGgtcgtc                                                                                                                                                                                           >  2:401930/1‑251 (MQ=255)
                                                        atGAAGCGCAAAATGATCCCATGCTGCCGGGATACTCGATTAATGCCAATCTGGTGGCGGGTTTAGCGCCGAATGAAGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGGATGACAGGTTATATTCTCAATATCACCATTCGCGGGCAGGGGGGGGTGAAAAATCAGGGACGCGAATTTGTTTGCCGCCCGGGTGGTATTTTGCTGTTCCCGCCAGGCGAGATTCATCACTCCGGTCGTCAt                                                                                                                                                                                         >  2:594239/2‑251 (MQ=255)
                                                          gAAGCGCAAAATTACCCCCTGCTGCCGGGATACTCGTTTAATTCCCATCTGGTGGCGGGTTGAACGCCGATTGAGGCAATCGTTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTGGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCAt                                                                                                                                                                                         <  2:45825/249‑1 (MQ=255)
                                                          gAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGTGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATc                                                                                                                                                                                        <  1:698269/250‑1 (MQ=255)
                                                                   aaTGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTcg                                                                                                                                                                              >  1:163488/1‑251 (MQ=255)
                                                                      gATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCTCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGGGATATTTTGCTGTTCCCGCCAGGGGAGATTCATCACTACGGTCGTCATCCGGAGACTcgcg                                                                                                                                                                            >  2:72119/1‑250 (MQ=255)
                                                                         ccccTGCTGCCGGGATACTCGTTTATTGCCCATCTGGGGGCGTGTTTAACGCCGATTGAGGCCAACGTTTATCTCGATTTTTTTTTCGACCGACCGCGGGGAATGAAAGTTTATATTCCCAATCGCCCCTTTCGCGTTCATGGGTTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGCCCGGGTGTTATTTGGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGCATg                                                                                                                                                                       <  2:924787/251‑1 (MQ=255)
                                                                         ccccTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTTCCTCGATTTTGTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCACACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATg                                                                                                                                                                       <  2:241931/251‑1 (MQ=255)
                                                                         ccccGGCTGTCGGGATACTCTTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACGGACCGCTGGGATTGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAAGCGGGGACGAGAATTTGTTTGCCGACCGGTTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATg                                                                                                                                                                       <  1:518670/251‑1 (MQ=255)
                                                                             tgctgcCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATggcat                                                                                                                                                                   >  1:723482/1‑248 (MQ=255)
                                                                                    gggATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATggaatacccagg                                                                                                                                                            >  2:99942/1‑241 (MQ=255)
                                                                                     ggATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGCCCGCTGGGGATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGTAAAATAAGGGACGAGTATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGCTTCATCACTACGGTCGTCTTCCGGAGGCTCGGAGATGGTATCACCAGTg                                                                                                                                                          >  2:370804/1‑251 (MQ=255)
                                                                                       aTACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGGCCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGGGAGATTCATCACTTCGGTCGTCATCCGGAGGCTCGCGAATGGTATCAACAGTggg                                                                                                                                                        >  2:841665/1‑251 (MQ=255)
                                                                                         aCTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTTAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGCTTCATCACTACGGTCGTCATCCGTGGGATCGCGAATgggatacacagggggtt                                                                                                                                                       >  2:197332/1‑236 (MQ=255)
                                                                                         aCTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGGTATTTTGCTGTTCCCGCCAGGCGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCCCCAGTggggt                                                                                                                                                      >  2:366821/1‑249 (MQ=255)
                                                                                         aCTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGtt                                                                                                                                                      >  1:503590/1‑251 (MQ=255)
                                                                                                 aatgCCCATCTGGTGGCGGGTTTAACGCCGCTTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCAATACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGGGGGCTTACTTTCg                                                                                                                                              >  1:665596/4‑251 (MQ=255)
                                                                                                 aatgCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAAGGAAAGGTTATATTCTCAATCTCACCATCCGCGGTCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTTTTTGCCGGCCGTGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTGGGGAATGGTATCACCAGTGGGTTTACTTTc                                                                                                                                               <  1:383798/247‑1 (MQ=255)
                                                                                                  atgCCCATCTGGAGGCGGGTTTAACGCCGATTGAGGGCAACGGTTATCTCGGTTTTTTTATCGACCGAACGATGGGAATGAAAGGTTATATTCTCAATCTCATCATTCGCGGTCAGGGGGTGGGGAAAAAACAGGAGCGTGAATTAGTTTGCCGACCGGGGGAGATTTTGCTTTTCCCGCACAGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGAGTATGGTATCACCCAGGGGTTTACTTTcat                                                                                                                                             >  2:861334/3‑249 (MQ=255)
                                                                                                    gccagtctggggtggggTTTACCGCCTATGGAGGCCAACGGTTATCACGTTTTTTTTTCCGCCCGCCGGGGGGGAATGAAAGGTTATATTCTCAATCTCACCATGCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGATTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTcc                                                                                                                                           <  2:600877/237‑1 (MQ=255)
                                                                                                    gCCCATCGGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTc                                                                                                                                            <  1:537717/250‑1 (MQ=255)
                                                                                                            ggtggCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCcgcgcg                                                                                                                                      <  1:161512/248‑1 (MQ=255)
                                                                                                                gggggTTTAACGCCGTGGGGGGCCAACGGTTAACTCGTTTTTTGTAACGCCCGCCCGCTGGGAATGAAAGTTTATATTCTCAACCTCACCATTCGCGGGCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTGTTGGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTg                                                                                                                               <  2:124726/249‑1 (MQ=255)
                                                                                                                  gggTTTAACGCCGATTGCGGCCAACGGTTATCTCGATTTTTTTATCGACAGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGGTTACTTTCGTCCGCGCGCCTACTGGc                                                                                                                             >  2:734044/1‑251 (MQ=255)
                                                                                                                         aCGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCCCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGAACGAGAATTTGTTTTCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCGTAGGGAGGCTCGCGTAAGGTATAACCAGTGGGTTTACTTTCGTCCGGGCGCCACCTGGCATGaagg                                                                                                                      >  2:76337/1‑249 (MQ=255)
                                                                                                                          cGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATgg                                                                                                                     <  1:841665/251‑1 (MQ=255)
                                                                                                                             cGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGGGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCtt                                                                                                                  >  2:281242/1‑251 (MQ=255)
                                                                                                                                   agtcaaaCGTTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCACACCATTCGCGGTCAGGGGGTGTTGAAAAATCAGGGACGAGAATTGGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGAGTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTgg                                                                                                            <  2:241085/246‑1 (MQ=255)
                                                                                                                                      ccAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGAGGTGAAAAATAAGGGACGAGGGTTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCTGAGGCTCGCGAATAGTATCACCAGTGGGTTTACTTTCGTACGCGCGCCTACTGGCATGACTGGCTTAGCTGGc                                                                                                           >  1:958637/1‑249 (MQ=255)
                                                                                                                                        aaCGGTTATCTCGATTTTTTTCTCGACTGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACAATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTACCTGGCCGt                                                                                                        >  2:417148/1‑250 (MQ=255)
                                                                                                                                        aaCGGGTATTTCGATTTTTTTATCGACCGCCTCCTGGGAAGTAAAGGTTGAATTCTCAAACTCTCCTTTCGGGGTCGTGGGTTGGTGAAAAATCAGGGACGTGATTTTGTTGGCCGCCAGGGTGAATTTTAGCTGTCCCCGCCAGGAGAGATTCATCACTACGGTCGCCATCCGGAGGCTCGCGTATGGTATCACCATTGGGTTTTCGTTCGCCCGCGCGCCTACTGGCATGAATTGCTTAACTGGCCgcc                                                                                                       <  2:579184/251‑3 (MQ=255)
                                                                                                                                          aggTTAACTAGATTTTTTAATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGGGAAAAATAAGGGAAGAGAATTTATTTGACGACCGGGTGATATTTTGCTGTTACCGCCAGGAGAGATTAATCACTACGGTCGACATCCGGAGGCGCGCGAATGGTATCACCAGAGGGTTTACTTTCGTACGCGGGCATACTGTCTAGCATGGCTTAACTGGCAGTCa                                                                                                      >  2:162372/2‑250 (MQ=255)
                                                                                                                                           ggTTATCTCGATTTTTTAATCAACCGACCGCTGGGAATGAAGGGTTATATTCTCAATCTAACCATTCGCGGTCAGGGGGTGGTGAAACATCCGGGACGAGCATTTGTTTGCCGACCGGGTGATATTATGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCCAGCGCATGGTATCAACAGTGGGTTTACTTTCGTCCGCGCGCCGACTGGCATGAATGGCTTCTCTGGCCGTCaa                                                                                                     >  2:620090/1‑250 (MQ=255)
                                                                                                                                              tATCTCGATTTTTTTATCGAACGACCGCTGGGAATGAAAAGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGCGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGACTACTGGCATGAATGAGTTAACAGGCCGTCAATAtt                                                                                                 >  2:885808/1‑251 (MQ=255)
                                                                                                                                                   cGATTTTTCAATCGACCGACCGCTTGGAATGAACGGTTATATTCGCAATTCATACATTCGCGGTTAGGTGGTGGTGAAAAATAAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTg                                                                                               <  1:284173/248‑1 (MQ=255)
                                                                                                                                                    ggtttttttATCGACCGACCGCTGGGTATGAAAGGTTATATTCTCAATCTCATCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGCGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCAGCCAGGAGAGATTCCTCACTACGGTCGTCATCCGGAGGCTAGCGAATGGTATCACCAGATGGTTTACGTTCGTCCGAGCGGCTACTGGCATGAAACGCTTAACTGACCGTCACTATTTgaca                                                                                            >  2:727419/3‑247 (MQ=255)
                                                                                                                                                        tttttATGGCCCGACCGTGGGGAAGGAAAGGTTATATTCTCAATCTTCCCATCGGCGGTGGGGGGTGGGTGAAAAACAAGGTAGGAGAATTTGTTTGCCGCCGGGGTGATATTTTTCTTGTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATAc                                                                                        <  2:235800/250‑1 (MQ=255)
                                                                                                                                                        tttttATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATAc                                                                                        <  1:221433/250‑1 (MQ=255)
                                                                                                                                                            tagcggccgaccgCTGGGATTGAAAGGTTATATTCTCAATCTCCCCATTCTCGGTCGGGGGGTGGTGAAAAAGCAGGGACGGGAATTTGTGTGCCAACCGGTTGATATTTTGCTGTTCCCGCCAGGAGAGATTAATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCATTGGGTTTACTTTCGTCCGCGCGCCTACTTGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACgggg                                                                                    <  1:767088/244‑2 (MQ=255)
                                                                                                                                                             aTCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACggggtt                                                                                  >  1:734419/1‑248 (MQ=255)
                                                                                                                                                                 accgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATAAGGGgttcttt                                                                              >  1:435245/1‑251 (MQ=255)
                                                                                                                                                                  ccgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTGGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctttc                                                                             <  1:99942/251‑1 (MQ=255)
                                                                                                                                                                   cgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCg                                                                            <  1:823378/251‑1 (MQ=255)
                                                                                                                                                                    gaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTTTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGc                                                                           >  1:220021/1‑251 (MQ=255)
                                                                                                                                                                           gggaagaaaaGGTTATATACTCAATCTCCCCGTTCGCGGTCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATg                                                                    <  1:370804/244‑1 (MQ=255)
                                                                                                                                                                                tGAAAGGGTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGAAGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAATGGGTTTACTTTCGTCCGCGCGCATCCTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGAGTTCTTTCGCCCGGACGAAgcg                                                               >  2:771795/1‑251 (MQ=255)
                                                                                                                                                                                       ttATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTAATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGCATGGCTTACCTGGCCGTCAATATTTGCCAATCCGGGGTTCTTTCGCCCGGATGAAGCGGACCAGc                                                        >  2:593926/1‑251 (MQ=255)
                                                                                                                                                                                        tataTTCTCAATCTCACAATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTAATCAATACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGGGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGATTGGCTTAAATGGCCGTCAAAATTTGCCACTACGGGATTCTTTAGCCCGGGTGAAgtcaccaagcc                                                       >  2:643800/1‑241 (MQ=255)
                                                                                                                                                                                                cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTc                                                >  1:230456/1‑250 (MQ=255)
                                                                                                                                                                                                 cattttccccATCCGCGGTCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGTGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCAttt                                                 <  1:748518/241‑1 (MQ=255)
                                                                                                                                                                                                 aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGCCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTc                                                >  1:401758/1‑249 (MQ=255)
                                                                                                                                                                                                     tcACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAc                                          >  1:458258/1‑251 (MQ=255)
                                                                                                                                                                                                      cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGTGGCTCGCGACTGGTATCACCAGAGGGTTTACTTTCGCCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTTCCAATACAGGGTTCTTTCGCCCGGATGAAGCGAACAAGCCGCGTTTccactaac                                         >  2:130801/1‑244 (MQ=255)
                                                                                                                                                                                                      cACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCCCCAGTGGGTTTCCTTTCGTCCGCGCGCCTACTGGCATGACTGGCTTAACTGGCCGTCCATATTTGCCAATACGGGGGTTTTTCGCCCCGATGATGCGCACCAGCCGCCTTTCAGCGAcc                                         >  2:580851/1‑251 (MQ=255)
                                                                                                                                                                                                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGTGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGACATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCt                                        >  1:822782/1‑251 (MQ=255)
                                                                                                                                                                                                       aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCt                                        >  1:208750/1‑251 (MQ=255)
                                                                                                                                                                                                            ttCGCGGTAAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                                     <  1:774079/250‑1 (MQ=255)
                                                                                                                                                                                                              tggaggtgggggtgtgggaaaaaaCCAGTGAGGAGAATTTTTTTGCCACCGGGGTGATTTTTTGTTTTTCCCGCCAGGAGAGATTCATCACTATGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAAGGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                                    <  2:749589/231‑1 (MQ=255)
                                                                                                                                                                                                                 ggTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGc                                >  1:707232/1‑250 (MQ=255)
                                                                                                                                                                                                                  gTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAAAACGGGGTTCTTTCGCCCGGCTGAAGCGCAACAGCCGCATTTCAGCGACCTGTTTGGGCa                               >  2:737934/1‑250 (MQ=255)
                                                                                                                                                                                                                    ggggggggggTGAAAAGTCAGAGGTGAGAATTTGTTTGCCGGCCGGGGGATATTGTGCTGTTCCCGCCAGGAGAGATTCAGCACTACGGTCGTCCTCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt                            <  2:441589/249‑1 (MQ=255)
                                                                                                                                                                                                                    cAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCCCCAGCCGCATTTCAGCGACCTGTTTGGGCAAAt                            >  1:100365/1‑251 (MQ=255)
                                                                                                                                                                                                                    cAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTAGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGAGCACCAGCCGCATTTCAGCGACCTGTTTGGGCGAAt                            >  2:862187/1‑251 (MQ=255)
                                                                                                                                                                                                                     aGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGGGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATCTTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATc                           >  2:889129/1‑251 (MQ=255)
                                                                                                                                                                                                                       gggttggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGTGGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCTCCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAt                         <  1:307309/247‑1 (MQ=255)
                                                                                                                                                                                                                        gggTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAtt                        >  1:372619/1‑251 (MQ=255)
                                                                                                                                                                                                                          gtggtgAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCACGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCCGCCGCATTTCAGCGACCTGTTTGGGAAAATCATTa                       >  2:969186/1‑250 (MQ=255)
                                                                                                                                                                                                                            ggtgAAAAATCAGGGACGAGAATTTGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACg                    <  1:324811/251‑1 (MQ=255)
                                                                                                                                                                                                                             gggaaaaaaTCAGGGACGAAAATTTGTTTGCCGCCCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACg                    <  1:163689/247‑1 (MQ=255)
                                                                                                                                                                                                                             gtgAAAAATCAGGGACGAGAATTTGTTGGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGc                   <  1:62768/251‑1 (MQ=255)
                                                                                                                                                                                                                               tgAAAAATCAGGGACGCGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAAAACGGGGTTCTTTCGCCCGGATGAAGCGCAACAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGcc                  >  2:922921/1‑251 (MQ=255)
                                                                                                                                                                                                                                gAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCcg                                                                                                                                          >  1:419026/1‑130 (MQ=255)
                                                                                                                                                                                                                                gAAAAAGCAGGGACGAGAATTGGTTTGCCGACCGGTTGATATTTGGCGGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCcg                                                                                                                                          <  2:419026/130‑1 (MQ=255)
                                                                                                                                                                                                                                aaaaaaaCAGGGACAAGAATTTTTTTCCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCCCTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCg                 <  2:734419/250‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaaaTCAGGGACGAGAATTTGTTTGCCGACCGGGTGTTATTTTGTTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCgg                <  1:719777/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaTCAGGGCGGAAAATTTTTTTGCCGACCGGTTGTTTGTTTGTTGTCCCCGCCAGGAGAGATTCATCACTTCGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGc              <  2:329620/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTCAGGGACGCGAATTTGATTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGCGATTAATCACTACGGTCGTCATCCGGAGGCTCGAGAATGGTAACACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGGTCTTTCGCCCGGATGCAGCGCACAAGCCTCCTTTCCGCGACCTGGTTGGACAAATCATTAAAGACGGGCa             >  2:110829/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTCTGGGCAAATCATTAACGCCGGGCAAgg          >  2:667915/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGTTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAggg         <  1:580851/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gagaATTTTTTTGCCGACCGGTTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAggg   <  1:112102/251‑1 (MQ=255)
                                                                                                                                                                                                                                               agaATTTGTTTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGgc  >  1:750493/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
GATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGC  >  NC_000913/70335‑70822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: