Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,757,000 G→T 5.5% E435* (GAA→TAA)  pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,0000GT 5.5% 282.9 / 2.0 81E435* (GAA→TAA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base T (3/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGT  >  NC_000913/1756864‑1757134
                                                                                                                                        |                                                                                                                                      
aTACTTCCCGGATGCAACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGTGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAAc                                                                                                                                     <  2:776357/140‑1 (MQ=255)
   cgtCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTgg                                                                                                                                  <  1:465129/138‑1 (MQ=255)
   cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAATAACTgg                                                                                                                                  >  2:128144/1‑140 (MQ=255)
   cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAATAACTgg                                                                                                                                  >  2:128825/1‑140 (MQ=255)
   cTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTgg                                                                                                                                  <  1:454748/140‑1 (MQ=255)
      cccGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGGACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGAGCACTGCTACGGATGATTTCTACCGTCTGGGTCAAGAACTGGctc                                                                                                                               >  1:146177/1‑140 (MQ=255)
                 ccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt                                                                                                                    <  1:666417/140‑1 (MQ=255)
                  caTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTc                                                                                                                   <  1:1507597/140‑1 (MQ=255)
                     ccTGGCCCTGACCTCCAACGACACATCGGCTCAGCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCTGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTgg                                                                                                                <  2:125322/140‑1 (MQ=255)
                         gCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcaca                                                                                                            >  1:723653/1‑140 (MQ=255)
                           aCTCACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcacaca                                                                                                          <  1:89988/140‑1 (MQ=255)
                               accaccAACGAATAAACGGCTCAGCAGTTGGTTCTGAGCAAAGGCGTTGTGCCGCAGCTTTTTAAAGAGATCACTTCTACTGTTGATTTCTACCGTCTGGGTAAAGATCTGTCTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                      <  2:1028601/140‑1 (MQ=255)
                               accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAATGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                      >  1:1156856/1‑140 (MQ=255)
                               accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                      >  2:1686549/1‑140 (MQ=255)
                               accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                      <  1:773055/140‑1 (MQ=255)
                               accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                      >  2:168653/1‑140 (MQ=255)
                               accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAgg                                                                                                      >  1:648033/1‑140 (MQ=255)
                                          aaaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAg                                                                                           <  2:630630/140‑1 (MQ=255)
                                          aaaaaCGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAg                                                                                           <  2:634535/140‑1 (MQ=255)
                                                  cTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGtt                                                                                   >  2:947532/1‑140 (MQ=255)
                                                    catcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTc                                                                                 <  1:596189/140‑1 (MQ=255)
                                                    catcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGATTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTc                                                                                 <  1:1686549/140‑1 (MQ=255)
                                                       caGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTTCTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                              <  2:32547/140‑1 (MQ=255)
                                                       caGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                              <  1:1344122/140‑1 (MQ=255)
                                                       caGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                              <  2:901614/140‑1 (MQ=255)
                                                       caGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                              <  2:904120/140‑1 (MQ=255)
                                                            ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCAc                                                                         <  2:988117/140‑1 (MQ=255)
                                                            ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCAc                                                                         <  1:1512367/140‑1 (MQ=255)
                                                            ggTACCGACCTAAGGCGTTGTGCCGCAGCTAGTTACACACAACACTTCTACTGATGATTTCTACCCTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTGGTAGTTATGGTTTCTGGTGCAc                                                                         <  2:1680599/140‑1 (MQ=255)
                                                                 tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTa                                                                    <  1:1386703/140‑1 (MQ=255)
                                                                 tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTa                                                                    <  1:370467/140‑1 (MQ=255)
                                                                 tGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTa                                                                    <  1:392578/140‑1 (MQ=255)
                                                                   aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCAGGGCTACTGTTGATGTCGTCCGTCGGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGAACACAAAGGTGACGTTGTAGTTATGGTTTCTGGGGCACTGGTAcc                                                                  >  2:686812/1‑140 (MQ=255)
                                                                   aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                                  >  2:268424/1‑140 (MQ=255)
                                                                   aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                                  >  2:1506060/1‑140 (MQ=255)
                                                                   aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                                  >  1:635903/1‑140 (MQ=255)
                                                                   aGCAAAGGCGTTGTGACGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCAGCAGAGCGGTCTGGAACACATAGGTTACTTTGTAGTTATGGTTTCTGGTGCACTGGTAcc                                                                  >  2:1249171/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  1:53005/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  1:498117/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  2:1515887/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  1:1072958/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  2:239143/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  2:239163/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  2:1208917/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  2:1188924/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  1:1156828/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGAGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCg                                                                 >  1:498127/1‑140 (MQ=255)
                                                                    gCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACGTCTACTGATGTTTTCTAAAGTCTGGGTAAATAACTGGCTCTGCATAGCTGTCTGGCACACAAGGGTGTCGTTGTAGGTATGGGTTCAGGTGCACTGGTAACg                                                                 >  2:1644048/1‑140 (MQ=255)
                                                                        aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCg                                                             >  2:232676/1‑140 (MQ=255)
                                                                        aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGGAGGTAGGGGTTCTGGTGCACTGGTCCCGAGCg                                                             >  2:840270/1‑140 (MQ=255)
                                                                        aGGCGTTCAGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTAGTAGTTAAGGTGACTGCTGCACTGGTACCGAGCg                                                             >  2:232719/1‑140 (MQ=255)
                                                                          gCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                           <  1:637623/140‑1 (MQ=255)
                                                                          gCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                           <  1:77970/140‑1 (MQ=255)
                                                                             ttGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCact                                                        >  2:863357/1‑140 (MQ=255)
                                                                               gtgCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTTCCGAGCGGCactac                                                      >  1:251065/1‑140 (MQ=255)
                                                                                  ccGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTaa                                                   <  1:1060240/140‑1 (MQ=255)
                                                                                     cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                                <  1:516949/140‑1 (MQ=255)
                                                                                     cAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAacac                                                >  1:263388/1‑140 (MQ=255)
                                                                                           gTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAATAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                          <  1:128825/140‑1 (MQ=255)
                                                                                           gTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAATAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                          <  1:128144/140‑1 (MQ=255)
                                                                                                agagATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                     <  1:1593904/140‑1 (MQ=255)
                                                                                                   gaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                  <  1:613778/140‑1 (MQ=255)
                                                                                                   gaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                  <  2:1145380/140‑1 (MQ=255)
                                                                                                   gaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                  <  2:439279/140‑1 (MQ=255)
                                                                                                   gaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                                  <  2:622398/140‑1 (MQ=255)
                                                                                                    aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGt                                 >  1:673624/1‑140 (MQ=255)
                                                                                                     tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTTACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTc                                <  1:168653/140‑1 (MQ=255)
                                                                                                     tCACTTCTACTGAGGATTTCTACCGTCTGGGTAAAGTACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTCGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTc                                <  2:1675284/140‑1 (MQ=255)
                                                                                                                 atgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCTCTACTAACACCGCATCTGTTCACGTCCTGTAATATTGc                    >  1:1125015/1‑140 (MQ=255)
                                                                                                                   gatTTCTACCGTCTTGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCtt                  <  2:158430/140‑1 (MQ=255)
                                                                                                                   gatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCtt                  >  2:1474501/1‑140 (MQ=255)
                                                                                                                   gagttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGTCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTCCTAACACCGCTTCTTTTCACATCCTGTTATATTGCtt                  <  2:158423/137‑1 (MQ=255)
                                                                                                                      ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTTGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTtg               >  1:1163857/1‑140 (MQ=255)
                                                                                                                      ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTGTGGTTTCTGGTGCACTGGTACCGAGCGGc                                                           <  1:1246345/96‑1 (MQ=255)
                                                                                                                      ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTGTGGTTTCTGGTGCACTGGTACCGAGCGGc                                                           >  2:1246345/1‑96 (MQ=255)
                                                                                                                             gTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaatta        <  2:669347/140‑1 (MQ=255)
                                                                                                                             gTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaatta        <  2:670935/140‑1 (MQ=255)
                                                                                                                             gTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaatta        <  2:669196/140‑1 (MQ=255)
                                                                                                                              tCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                           >  2:793466/1‑88 (MQ=255)
                                                                                                                              tCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGc                                                           <  1:793466/88‑1 (MQ=255)
                                                                                                                                   gTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTga  >  1:888656/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATACTTCCCGGATGCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGT  >  NC_000913/1756864‑1757134

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: