Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 977,748 T→G F475V (TTT→GTT)  mukB → chromosome condensin MukBEF, ATPase and DNA‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913977,7480TG98.9% 356.3 / ‑3.5 92F475V (TTT→GTT) mukBchromosome condensin MukBEF, ATPase and DNA‑binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (44/47);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

TGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGC  >  NC_000913/977600‑977894
                                                                                                                                                    |                                                                                                                                                  
tGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGgtt                                                                                                                                                  >  2:125455/1‑151 (MQ=255)
 gCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTg                                                                                                                                                 >  1:541370/1‑151 (MQ=255)
      aGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGc                                                                                                                                              >  1:980798/1‑149 (MQ=255)
      aGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAg                                                                                                                                            <  1:474490/151‑1 (MQ=255)
      aGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAg                                                                                                                                            <  2:554598/151‑1 (MQ=255)
        aaCTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGc                                                                                                                                          >  2:1041509/1‑151 (MQ=255)
        aaCTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGc                                                                                                                                          >  2:349002/1‑151 (MQ=255)
         aCTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCt                                                                                                                                         <  1:1087857/151‑1 (MQ=255)
               ccATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAg                                                                                                                                   <  2:1102793/151‑1 (MQ=255)
                     gCCGTACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCtg                                                                                                                                <  2:1019986/148‑1 (MQ=255)
                     gCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTAAGCAGGCTTATCAGCtggt                                                                                                                              <  2:935856/150‑1 (MQ=255)
                          aCTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCtggtggtg                                                                                                                          >  1:11664/1‑149 (MQ=255)
                           cTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTCAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGc                                                                                                                        >  1:630102/1‑150 (MQ=255)
                              aaCCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATc                                                                                                                    >  1:587335/1‑151 (MQ=255)
                                 cgccgACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAAc                                                                                                                 <  2:1114305/151‑1 (MQ=255)
                                   cagACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACgg                                                                                                               >  1:831935/3‑151 (MQ=255)
                                     gACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGcc                                                                                                             >  2:696935/1‑151 (MQ=255)
                                       cTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGccc                                                                                                            <  1:1104502/150‑1 (MQ=255)
                                               ccgAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGcgcg                                                                                                   >  2:1126915/1‑151 (MQ=255)
                                                 gAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGt                                                                                                  >  2:1197312/1‑150 (MQ=255)
                                                  aaTGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTaa                                                                                                <  2:252628/151‑1 (MQ=255)
                                                  aaTGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTa                                                                                                 >  2:597893/1‑150 (MQ=255)
                                                  aaTGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTa                                                                                                 >  1:1037459/1‑150 (MQ=255)
                                                  aaTGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTa                                                                                                 <  1:730753/150‑1 (MQ=255)
                                                    tGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCTCGTAACg                                                                                              >  2:690086/1‑151 (MQ=255)
                                                     ggCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACg                                                                                              >  2:498617/1‑150 (MQ=255)
                                                     ggCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACg                                                                                              <  2:286686/150‑1 (MQ=255)
                                                     ggCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACg                                                                                              <  1:349745/150‑1 (MQ=255)
                                                     ggCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGa                                                                                             <  2:780372/151‑1 (MQ=255)
                                                      gCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACg                                                                                              >  1:1278155/1‑149 (MQ=255)
                                                       cTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAgg                                                                                           >  2:90320/1‑151 (MQ=255)
                                                        tGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGc                                                                                          >  2:306042/1‑151 (MQ=255)
                                                        tGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGc                                                                                          >  1:528203/1‑151 (MQ=255)
                                                         ggAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGc                                                                                          <  1:1126278/150‑1 (MQ=255)
                                                           aaaCCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTg                                                                                       <  1:299459/151‑1 (MQ=255)
                                                             aCCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTggg                                                                                     >  2:1286340/1‑151 (MQ=255)
                                                                ttCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATg                                                                                  >  1:606845/1‑151 (MQ=255)
                                                                 tCCAGGCGAAAGAGCTGGAAGTGACTGAAAAAATGCCCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGAt                                                                                   <  1:862412/149‑1 (MQ=255)
                                                                 tCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGAt                                                                                   <  2:276662/149‑1 (MQ=255)
                                                                 tCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGAt                                                                                   <  1:363194/149‑1 (MQ=255)
                                                                 tCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATg                                                                                  <  1:656721/150‑1 (MQ=255)
                                                                 tCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATg                                                                                  <  2:341134/150‑1 (MQ=255)
                                                                     ggCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGtcgct                                                                             <  2:733435/151‑1 (MQ=255)
                                                                        gAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGtcgctcgc                                                                          <  1:1102921/151‑1 (MQ=255)
                                                                        gAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGtcgctcg                                                                           >  1:68503/1‑150 (MQ=255)
                                                                        gAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGtcgctcg                                                                           >  1:671660/1‑150 (MQ=255)
                                                                         aaaGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTcgcg                                                                         <  2:908259/151‑1 (MQ=255)
                                                                             agCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACt                                                                     >  2:704828/1‑151 (MQ=255)
                                                                              gCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACt                                                                     >  1:719022/1‑150 (MQ=255)
                                                                              gCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTa                                                                    <  1:531243/151‑1 (MQ=255)
                                                                                 ggAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTAtt                                                                  <  2:1253006/150‑1 (MQ=255)
                                                                                 ggAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTAtt                                                                  <  2:644207/150‑1 (MQ=255)
                                                                                 ggAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTg                                                                 <  1:1039767/151‑1 (MQ=255)
                                                                                 ggAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTg                                                                 <  1:670525/151‑1 (MQ=255)
                                                                                  gAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTgc                                                                <  2:1078136/151‑1 (MQ=255)
                                                                                   aaGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTgcg                                                               >  1:1144593/1‑151 (MQ=255)
                                                                                   aaGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTgcg                                                               <  1:376801/151‑1 (MQ=255)
                                                                                       gACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGa                                                            >  1:739696/1‑150 (MQ=255)
                                                                                         cTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAgg                                                         >  2:1036774/1‑151 (MQ=255)
                                                                                          tGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAggg                                                        >  1:1114218/1‑151 (MQ=255)
                                                                                            aaaaaaTGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGt                                                      >  1:96191/1‑151 (MQ=255)
                                                                                              aaaaTGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCg                                                    <  1:623447/151‑1 (MQ=255)
                                                                                                 aTGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATc                                                 >  2:1299049/1‑151 (MQ=255)
                                                                                                      ctctctTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGt                                            >  1:1237996/1‑151 (MQ=255)
                                                                                                        ctctTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTc                                           >  2:116474/1‑150 (MQ=255)
                                                                                                            tGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCAc                                         <  2:11664/148‑1 (MQ=255)
                                                                                                            tGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCt                                       <  1:557896/150‑1 (MQ=255)
                                                                                                            tGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCt                                       <  2:1130559/150‑1 (MQ=255)
                                                                                                             gAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAATGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTgg                                     <  2:482879/151‑1 (MQ=255)
                                                                                                               gCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGc                                    <  1:317091/150‑1 (MQ=255)
                                                                                                               gCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGc                                    >  2:724116/1‑150 (MQ=255)
                                                                                                                   aaaaTGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGc                               <  2:256442/151‑1 (MQ=255)
                                                                                                                   aaaaTGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGc                               >  1:944826/1‑151 (MQ=255)
                                                                                                                    aaaTGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCa                              <  1:76389/151‑1 (MQ=255)
                                                                                                                     aaTGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGc                               >  2:645856/1‑149 (MQ=255)
                                                                                                                     aaTGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGc                               >  1:1132470/1‑149 (MQ=255)
                                                                                                                      aTGAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAgg                            >  2:355498/1‑151 (MQ=255)
                                                                                                                        gAGCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGtt                          <  2:274758/151‑1 (MQ=255)
                                                                                                                          gCATGGCGCAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACAGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTc                         <  2:260385/150‑1 (MQ=255)
                                                                                                                                cgcAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCt                  >  2:788088/1‑151 (MQ=255)
                                                                                                                                cgcAAACCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCt                  <  1:338844/151‑1 (MQ=255)
                                                                                                                                   aaaCCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCg               >  2:1104751/1‑151 (MQ=255)
                                                                                                                                   aaaCCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCg               <  1:1065117/151‑1 (MQ=255)
                                                                                                                                   aaaCCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCg               >  1:1277757/1‑151 (MQ=255)
                                                                                                                                   aaaCCGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCg               >  2:407483/1‑151 (MQ=255)
                                                                                                                                    atccGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCgg              <  1:1179688/149‑1 (MQ=255)
                                                                                                                                     aCCGCGCACAGCCAGGTTGATCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTCACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGGAGGGGTCGATCAGCGTCACCTGGCAGCGCAGGTTCAGCCGc                   >  1:547418/1‑145 (MQ=255)
                                                                                                                                      ccGCGCACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGAt            <  1:704802/151‑1 (MQ=255)
                                                                                                                                        gcgcACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGc          <  2:448754/151‑1 (MQ=255)
                                                                                                                                        gcgcACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGc          <  2:1208242/151‑1 (MQ=255)
                                                                                                                                         cgcACAGCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCg         <  1:1180391/151‑1 (MQ=255)
                                                                                                                                               gCCAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTa     >  1:376916/1‑149 (MQ=255)
                                                                                                                                                ccAGGTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGc  <  1:675459/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTATTGCGCGAAGGGGTCGATCAGCGTCACCTGGCAGAGCAGGTTCAGCCGCTGCGGATGCGATTAAGC  >  NC_000913/977600‑977894

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