Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,130,037 G→T 17.5% A31E (GCA→GAA)  flgM ← anti‑sigma factor for FliA (sigma 28)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,0370GT17.5% 183.5 / 7.8 96A31E (GCA→GAA) flgManti‑sigma factor for FliA (sigma 28)
Reads supporting (aligned to +/- strand):  ref base G (7/71);  new base T (16/2);  total (23/73)
Fisher's exact test for biased strand distribution p-value = 4.76e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCA  >  NC_000913/1129802‑1130274
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
tGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTtcc                                                                                                                                                                                                                                                   >  2:449596/1‑238 (MQ=255)
  gTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCTTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGTTTTTgcctc                                                                                                                                                                                                                                                 >  2:522317/1‑236 (MQ=255)
  gTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGCCGGTTGGTTTTtcccc                                                                                                                                                                                                                                                 >  2:365883/1‑236 (MQ=255)
  gTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGATTTCCTTGCGCGTCGCTTAACGACACCATTGTGCTGGTGGAGGCGGTTGGTTTTgtcgc                                                                                                                                                                                                                                                 >  1:86233/1‑238 (MQ=255)
       ggATCTCTGCCAGACGTGTCATACGGTCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTCCCGGGTTGCATCAGTTTTGCTTCCGCGCCGCTCAACGTCCCACTGGGGCTGGTGGAGGCGGTTGTTTTTGTCGCCCGTCt                                                                                                                                                                                                                                           >  1:229513/1‑239 (MQ=255)
          tctcctcaattcgtgcctttcGATCAGTTTCGTCTGAAGTGTGGCTGTGCTTCGTTGATCAGCGAGTCGGCAATTTTGCTCGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACGCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGtt                                                                                                                                                                                                                                       <  1:967943/219‑1 (MQ=255)
                catggcgttgtaTACGATCTTTTACTCTGCAAGTTTTGCTGCGCTTCGTTTATCCGCGCATCGGCCATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGCCACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTAc                                                                                                                                                                                                                               <  1:293746/228‑1 (MQ=255)
                     cggtttctaagATCAGTTTCCTTGCAAGTTTTGGTGCGCTTCGTTGGCAAGCGGATCGGCAATTTTGCCGGTTTCCATTTTTAGTTCACCGTTACGAATCGGCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCgg                                                                                                                                                                                                                            <  1:907550/227‑1 (MQ=255)
                        gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTTCCGCCCCGCCGGTCGTTaaccgccgc                                                                                                                                                                                                                        >  2:293746/1‑230 (MQ=255)
                        gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCCGTTCGTTTCCggggggc                                                                                                                                                                                                                       >  2:411457/1‑234 (MQ=255)
                        gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGATTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGAGGAGGCGGTTTCTTTTTCTGCCCGGCGGTTCGTTTCCgggcggt                                                                                                                                                                                                                       >  2:737278/1‑234 (MQ=255)
                        gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAAATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCCCCCTGGTGCTGGTGGAGGCGGGTGGTTTTTGCGCCCCGccgctccgttcgggggggg                                                                                                                                                                                                                          >  1:1278364/1‑221 (MQ=255)
                                 tCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTTCCCCCCGGCGGTTCGTTTACCGGGCGTCacagggtt                                                                                                                                                                                                             >  2:161356/1‑232 (MQ=255)
                                      tttctctggcggtttttctggtcTTTGTTTATCCGCGCTTTGGCTATTTGGGCGGGTTCCATTTTTAGTTACCCGTTACGAATCGTCAGTTTTAACGCTTCGACACGTTAAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGGCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcgcg                                                                                                                                                                                                        <  1:442254/217‑1 (MQ=255)
                                      tACTCTTGAAGTCTTGGTTCTGTTCATTGATCAGCGCATCGGCAATGTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATAGCACTGCTGCCGGGTTGCAGCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcgcg                                                                                                                                                                                                        <  1:787885/238‑1 (MQ=255)
                                      tACTCTGCAAGTCTTGCTGCTCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcgcg                                                                                                                                                                                                        <  2:176684/238‑1 (MQ=255)
                                      gtcagtgtcAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcg                                                                                                                                                                                                          <  2:31784/229‑1 (MQ=255)
                                       tacttggCAAGTCTTTCCGCGTTTGGTTGATCAGCCGATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGAAACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcgcg                                                                                                                                                                                                        <  1:107194/232‑1 (MQ=255)
                                        ctctGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCAACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCgg                                                                                                                                                                                                     <  2:1279397/239‑1 (MQ=255)
                                        ctctGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCgg                                                                                                                                                                                                     <  2:435903/239‑1 (MQ=255)
                                          tctGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTGTTTTCCGCCCGCGTGGTCGTTTCCGGGGCGGCAGTGGGTTTCGCgggg                                                                                                                                                                                                     >  2:510030/1‑239 (MQ=255)
                                          tctGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGGGCTGGTGGAGTCGGTTTTTTTTTCCGCCCCGCTGTTCCTTTCCCGGCGCGGACGTGGtttgtcgggg                                                                                                                                                                                                     >  1:206520/1‑232 (MQ=255)
                                            tGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGGGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGAGGTTGTTTTTTGTGCCCGGCGGTTCGTTGCCGGCGCGtacggggggtagggcgggtt                                                                                                                                                                                                   >  1:93976/1‑220 (MQ=255)
                                            tGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTTCCGCCCGGCCGTTCGTTACCGGcccgcccgggggtttccgcgggtt                                                                                                                                                                                                   >  2:351513/1‑219 (MQ=255)
                                            tGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCCGTTCGTTTCCCGGGCGGCAGTGGGTTTGCGCGGGTg                                                                                                                                                                                                  >  2:211056/1‑239 (MQ=255)
                                            tGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGGTTTTTCCGCCCGGCTGTTTCTTTACCgggggtcccggggtttcgccgggtg                                                                                                                                                                                                     >  2:147969/1‑216 (MQ=255)
                                            tGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGTTTTTGCCCTCCCGCCGTTCGTTACCGGCGCGGCCGGGGGTTCGCGCGGCTg                                                                                                                                                                                                  >  2:558772/1‑239 (MQ=255)
                                             ggaaGTCTTGCTGCGCTTCGTTGATCAGAGCATCGGCAATTTTGCCGGTGTCCATTTTTATTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGa                                                                                                                                                                                                 <  2:647488/237‑1 (MQ=255)
                                             cgaaGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTTGCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGa                                                                                                                                                                                                 <  2:1014210/237‑1 (MQ=255)
                                               cttggttGCTCCGCTTCGGGTATCAGCCCCACTGGAATTTTGCCGGTGTCCTTTTTTAGTTTCCCGGTCCGAATCGCCAGTTTTAACGGTTCGACCCGTTTAAGTATGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGTTGTAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTATCGGCGACTCAGTGGTTTCGCGCTGTTGAAc                                                                                                                                                                                               <  1:510030/234‑1 (MQ=255)
                                                aGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGATGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGCGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTATCGGCGCGGCAGTGGGTTTGCCCGGGTGGACg                                                                                                                                                                                              >  2:793071/1‑239 (MQ=255)
                                                   cTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGGTTTTTCCGCCCGGCCGTTCGTTTCCGGCCGGGTGGGGGGTTTcgcgggtgtgaccgg                                                                                                                                                                                             >  2:103475/1‑227 (MQ=255)
                                                    ttGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGGTTTTTTCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTGTTCCGCGGGTGAACggg                                                                                                                                                                                            >  2:1280554/1‑236 (MQ=255)
                                                    ttGCTGAGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGTTTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTTGCTTAACCTAGCACTGTGGCGGGCGGAGTGGGTTGTTTTTTTCGCCCCGCCGGTGGTTACCGTCCCGTGAGTggggtaccgcgggtgttcggggg                                                                                                                                                                                          >  1:207724/1‑218 (MQ=255)
                                                      ggtgcGCTTCGTTTATCAGCGCATTGGCATGTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCtt                                                                                                                                                                                        <  1:912691/237‑1 (MQ=255)
                                                      gctgcGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGTTTTTTCCGCCCGGCTGTTCCTTTACGGGGCGCCAGGGGTGTCCGGCGGTTgtatgggggtt                                                                                                                                                                                        >  2:790051/1‑229 (MQ=255)
                                                      gctgcGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGCCGGTTGTTTTTTCCGCCCGGCTGTTCGTTTCCCGCCCGTCCCGGGTTTCGCGCGGGTTAAccgggcgt                                                                                                                                                                                        >  2:732334/1‑232 (MQ=255)
                                                       ctgcGTTTCGTTGATCAGCGCAGCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                       <  2:1035084/239‑1 (MQ=255)
                                                         tgccttcgtttgtccgggctttcgCAATTTTTCCGGTTGCCATTTTTAGTTCACCGTTACGAATCGGCAGTTTTAACGCTTCGACACGTTCAAGATTGATAGCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACa                                                                                                                                                                                     <  1:453117/216‑1 (MQ=255)
                                                         tcgtgtgtgtttttcggccATCGGGAATTTTGCCGGGTTCCATGTTTAGGTCCTCGTTACGAAGCGGCAGTTTTAACCGTTCGGTACGGTCAAGATTGATATCACTGCTTCCGGGTTGCATCAGTTTTGCGTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACa                                                                                                                                                                                     <  1:104542/221‑1 (MQ=255)
                                                         gcgcTTGGTTGATCAGGGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGAGTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACa                                                                                                                                                                                     <  2:231675/239‑1 (MQ=255)
                                                          cgcTTCGTTGATCAGCGTATCGGCAATTTTGCCGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAg                                                                                                                                                                                    <  2:219873/239‑1 (MQ=255)
                                                                         cgcATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCCGTTTGTTTCCGGCGCGCCAGGGGGTTTGCGCGGGTGGACCGGGGTTTCAGGCtttccaggggaga                                                                                                                                                                     >  2:352071/1‑228 (MQ=255)
                                                                         cgcATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGGTTTTTCCGCCCGGCTGTTCGGTTCCGGCGCGCCACGGGTTTTGCGCGGGTTAAAGGGGCTTCACGGCGTCGGAGGGGaaa                                                                                                                                                                    >  2:1266413/1‑238 (MQ=255)
                                                                          gcATCGGATTTTTTTCCGGTGTCTCTTTTTTGTTTACCTTTACGAATCGCGAGTTTTAACGCTTCGACACGGTCAAGATTGATATCGCTGCGGCCGGGTTGCAGCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGt                                                                                                                                                                   <  1:161356/239‑1 (MQ=255)
                                                                            ccgcgatgttttGCGGGTTCCTTTTTTGATTCACCCATACGAATCGCCGGTTTTAACTCTTCGGCACGTTCAAGATTGATAACCGTGCTGGCGGGTTGCATCAGTGTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTTGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTAACGGCGCTTCAGTGGTTTCGCGCGGTTGAACGGGGCTTACAGGCTTCAGAGGGGAAGAGCg                                                                                                                                                                <  1:666514/231‑1 (MQ=255)
                                                                                tcatttttGCCTGTTCCCTTTTTTTGTTCCCCGTTTCGAATCGCCAGTTTTAACGCTTCCACAGGTTCAAGATTGATATCACTGCTGCCGGGTTGCAGCAGTTTTGCGTGCGAGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATc                                                                                                                                                             <  1:1182476/235‑1 (MQ=255)
                                                                                   tattgggcGGTGTCCATTTTTAGTTCCCCCTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCaa                                                                                                                                                           <  2:842250/231‑1 (MQ=255)
                                                                                       tGCCGGTTGCCATTTTTAGTTCCACCTTACAAATCGCCAGTTTTAAATGCGCGACACGTTTAAGATTGATATCACTGCGGCCGGGTTGCATCAGTTTTGGTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACCGGCTTCAGAGGCGAAGTGCGATCAATACTc                                                                                                                                                      <  1:1282619/239‑1 (MQ=255)
                                                                                        ttctgtttCCATTTTTTGTCTACTGTTTCGACACGCCAGTTTTTACGCTTGGAGTCGTGTGAGATTGTGATCCATGCTGCCGGGTTGGATCAGGTTTGCTTTCTCGTTGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGGCGCCCGGCTGTTCGTTACCGGCGCTTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTTAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  1:558772/232‑1 (MQ=255)
                                                                                        tgcgggcgccATTTTTTGTTGCCCGTTATGAATCGCTTGTTTTAACGCTTCGTCACTTTCAAGATTGATACCGCTGCTGGCGGGTTGCAGCTGTTTAGATTGCGCGTCGCTTAGCGTCACACTGGTGCGTGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTAGTTTCCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  1:1182863/231‑1 (MQ=255)
                                                                                        ggcGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGAAACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  2:1180879/237‑1 (MQ=255)
                                                                                        gcctgttcccttttttAGTTTCCCGGTTACAATCGCCAGTTTTAACGGTTCGACGCGTTCAAGATTGATACCCCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCGGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  1:411457/228‑1 (MQ=255)
                                                                                        gCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  2:229513/239‑1 (MQ=255)
                                                                                        gCCGGGTTCCATTTTTAGTTCCCCGTTACGAATCACCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  2:1020722/239‑1 (MQ=255)
                                                                                        gCCGGGTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTGTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                     <  2:891739/239‑1 (MQ=255)
                                                                                          ctttgtCCCTTTTTAGTTTACCGTTACGAATCGCAAGTTTTAACCCTTCGACCTGTTGAAGGTTTGTACCCCGGGTGCCGGGTTGCATCGAGTTTGCTTGCGCGTCGCTTAGCGGCACACTGGTGCTGGGGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTGGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGTCGAAGTGCGATCTATACTCATg                                                                                                                                                   <  1:790051/236‑1 (MQ=255)
                                                                                             tgatcttttttCGTTCCCCGTTCCGCATCGCTAGTTTTAACGCTTCTACCCGTTCAAGATTGGTATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTTTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGt                                                                                                                                                 <  1:1011344/232‑1 (MQ=255)
                                                                                               tttctttttttttttcccttttcGAATCGCCGGTTTTAACGCTTCGACCCGTTCAAGATTGATATCTCGGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGttt                                                                                                                                               <  1:908804/217‑1 (MQ=255)
                                                                                               ttctttttttgtttccctcttACGATTCGTCAGTTTTAACGCTTCGGCTCTTTCAAGCTTGATATCCGTGCTGTCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACTCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGGTGTTCGTTACCGGCGCGTTAGTGGTTTCGCGCGGTTGATCGGTGATTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGttt                                                                                                                                               <  1:163740/220‑1 (MQ=255)
                                                                                               tCCTTTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACTCGTTCAAGGTTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCCTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTa                                                                                                                                              <  2:639019/239‑1 (MQ=255)
                                                                                                  aTTTTGAGTTCACCGTTACGAATCGCCATTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTGTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTc                                                                                                                                           <  2:86233/239‑1 (MQ=255)
                                                                                                        aGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCAtt                                                                                                                                     <  2:706984/239‑1 (MQ=255)
                                                                                                              tcccttACGATTCGCCAGTTTTAAAGCTTCGACACGTTCAAGATTGATATCACTGGTGCCGGGTTGCACTAGTTTTGCTTGCGCGGCGCTTAACGTCTCACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                <  1:1114753/235‑1 (MQ=255)
                                                                                                              ccctttccgattcctcgtttttAACGGTTCTGCTCGTGAAAGTTTGTTTTCCCTGCTGCCGGTTTGCATCAGTTTTTCTTGCGCGTCGCTTAGCGTCACACTGGTGCTTGTGGAGGCGGTTGTTTTTGCCGCTCGTCTGGTCGTTACCGGCGCGTCAGTGGTTTCGCGAGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTACTCCTCATTGAagg                                                                                                                                <  1:737278/221‑1 (MQ=255)
                                                                                                                cttttcgaagtctcAGTTTTAACCCTTCGACACGTTCAAGATTGGTATCCCGGCGGCCGGGTTGGATCAGTTTTGCTTGCGCGTCGCTTAACGTTACACTGGTGCGGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGgcg                                                                                                                              <  1:799607/226‑1 (MQ=255)
                                                                                                                cGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGgcg                                                                                                                              <  2:206520/239‑1 (MQ=255)
                                                                                                                cGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGgcg                                                                                                                              <  2:857889/239‑1 (MQ=255)
                                                                                                                cGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGgcg                                                                                                                              <  2:199680/239‑1 (MQ=255)
                                                                                                                  tttcGACTCGCCAGTTTTTACGCTTGTAACCGTTCAAGATGTTAATCGCTGGTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCt                                                                                                                            <  1:147969/236‑1 (MQ=255)
                                                                                                                   tACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCtt                                                                                                                           <  2:1278364/239‑1 (MQ=255)
                                                                                                                        actgCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATc                                                                                                                      <  2:272817/236‑1 (MQ=255)
                                                                                                                            cctgTTTTAACGTTTGCACACTTGCAAGATTGGTATCCCTGCTGCCGGGTTGCTTCAGTTTTTTTTGCGGGTCGCTTAACGTCGCACTGGTGGTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCAtgt                                                                                                                  <  1:1226327/236‑1 (MQ=255)
                                                                                                                             tgttttttaggtttCGACACGTTCAAGATTGATTTCCCTGCATCCGGGTTGCATTGGTTTTGCTTGCGCGTCGCTTAACTTCACACTGGTGCTGGTGGAGTCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCAagtg                                                                                                                 <  1:423474/227‑1 (MQ=255)
                                                                                                                                  tttaCGCTTCGGCATTTTCAAGATTTATATCACTGCTTGCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCt                                                                                                            <  1:487340/236‑1 (MQ=255)
                                                                                                                                     aCGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCttat                                                                                                         <  2:246605/239‑1 (MQ=255)
                                                                                                                                              cacTTTCAAGATTGATATCCCTGGTGCCGGGTGTGCACGTTTTTGCTTGGGGGGCGCTTAACGTCCCACTGGTGCTGGGGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCa                                                                                                <  1:667337/239‑1 (MQ=255)
                                                                                                                                                         ttGATATCACTGCTGCCGTGTTGCATCTGTTTTGCTTGCTGGTCGCTTAACGTCACGCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACggg                                                                                       <  2:207724/237‑1 (MQ=255)
                                                                                                                                                         ttGATATCACTGCTGACGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTa                                                                                     <  2:847119/239‑1 (MQ=255)
                                                                                                                                                                 ccttctgctgGGTTGCGTCTTGTTTGCTGGGGCGCCGGTTTACGTCCCCATGGTGCGGGTGGAGGCGGTTGTTTTGGCCGCCCGGCTGTTCGTTGCCGGGGCGTCAGTGGGTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGAGGGGTAATCTTTaa                                                                             <  1:1266413/235‑1 (MQ=255)
                                                                                                                                                                 aCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCaa                                                                                                                                                           <  2:925178/161‑1 (MQ=255)
                                                                                                                                                                 aCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGGTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTTAGGGCGCTTTTTTCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                             >  2:1221662/1‑239 (MQ=255)
                                                                                                                                                                         cGGGTTTCATCAGTTTTGCTTGCGCGTCGCTTAACGTCCCACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                     <  2:93976/239‑1 (MQ=255)
                                                                                                                                                                         cGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTa                                                                      <  2:813826/238‑1 (MQ=255)
                                                                                                                                                                          cggttttgttctgTTTTTCTTGCGGGTTGTTTAACGGAACACTGGTGCTGGGGGAGGCGTTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                     <  1:1280554/227‑1 (MQ=255)
                                                                                                                                                                             ttGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGGTTTTTGCGCCCGGCTGGTCGTTACCGGCGCGGCCGGGGGTTTGCCCGGGTTTAGGGGGCTTCCATGGTTTAGCGGGCACGGGGGGGCAAAATCCATGGGTATTTCCCCCTTGGTGGCGCTTTTTTCCTGTGTTGTTGTTTTTTCCGGAAGGGACGGGtgactttttaaacggttacggg                                                                 >  1:786696/1‑218 (MQ=255)
                                                                                                                                                                                   tCAGTTTTGCTTGCGTGGCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTtat                                                             <  2:98134/238‑1 (MQ=255)
                                                                                                                                                                                        tttGCTTGCGCGTGCCTTAACGTTACACTGGTGCTGGTGGAGGGGGTTGTTTTTGCCGCCCGGCTGTTCGTTTCCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATTCTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGa                                                         <  2:646125/237‑1 (MQ=255)
                                                                                                                                                                                                              tgtcgttggttttggttgGGGTGGTTGTTTTTTCCGCCCGGCTGTTCGTTACCGGCGCCTGAGTGTTTTCGGTCGGTTGAACGGTGCTTACAGGCTTAAGAGGCGAAGTGCGATCAATTCTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGTCAAGGGACGTGTAATCTTTAACAGCTTACAGTTTTATAAGAATATTCCCATCTGCATCAACAAcg                                 <  1:336522/227‑1 (MQ=255)
                                                                                                                                                                                                               cacgtcgtggcgggtggAGGCGTTTTTTTTTGCCGCTCTGCTTGTCGTTGCCGGGGCGTGGGTGGTTTTGCGCGGTTGAACGGTGCTTACAGGCATCAGGGGCGAAGTGCGATGAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTGTATAAGAATATTCCCATCTGCATCAACAAcga                                <  1:705601/227‑2 (MQ=255)
                                                                                                                                                                                                                 ctcggttggtggtggAGGGGGTTGTTTTTGCCGTCCGGGTGTTTGTTACCGGCGCGTCAGTGGTTTCGCGCTGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATTCCCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCCTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAAcgccg                              <  1:896561/230‑1 (MQ=255)
                                                                                                                                                                                                                  acTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                             <  2:454905/191‑1 (MQ=255)
                                                                                                                                                                                                                    tggtgctggtgGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTg                           <  2:786696/239‑1 (MQ=255)
                                                                                                                                                                                                                    tggggctggggGGGGGGGTTTTTTTGGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGGGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTg                           <  1:103475/239‑1 (MQ=255)
                                                                                                                                                                                                                          tggtgtGGGGGGTTGTTTTTGCCTCCCGGCTGTTTGTTACCGGCGCGTCAGTGGTTGCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                     <  1:352071/239‑1 (MQ=255)
                                                                                                                                                                                                                          tggtggAGGCGGTTGTTTTTGCCGCCCGGCTTTTCGTTACCGGTGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                     <  2:1195795/239‑1 (MQ=255)
                                                                                                                                                                                                                          gggtggGGGTGGTTGTTTGTGCCGCCGGGCTGTTCGTGTCGGGGGTGTTAGGGGTTTCGGGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATATTCATGTTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCGTACAGGTTTATAAGAATATTCCCATTTGCAGCAATACAGCCGGTGGCTAcc                     <  1:732334/238‑1 (MQ=255)
                                                                                                                                                                                                                           ggtggGGGGGTTTTTTTTTTCCGCCCGGCTTTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                     <  1:1221662/238‑1 (MQ=255)
                                                                                                                                                                                                                                           tctGCCGCCCGGGTTGTTGTTACTGGCCCTGTAATGGTTTCGCGCGGTGGAATGGTGTTTACCGGCTGCAGGGGTCGAGGAGGATGAATCCCCATGGTTTATTCCTGATTGATGTAGCTTTTATAATGTGTTGCTTATTTATCGGCCAGGGACGGGTTATCTTTAACAGCCTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATgcg    <  1:449369/237‑1 (MQ=255)
                                                                                                                                                                                                                                             tGCCGCCCGGCTGTTGGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCa  <  2:1188929/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
TGGTCGAGGATCTCTGCAAGACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCA  >  NC_000913/1129802‑1130274

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