Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,606,883 G→A E145K (GAG→AAG)  zntA → zinc, cobalt and lead efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,606,8830GA99.2% 470.3 / ‑3.9 121E145K (GAG→AAG) zntAzinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (62/58);  total (63/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 121 total aligned reads displayed.

ATTCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGGAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTAT  >  NC_000913/3606736‑3607025
                                                                                                                                                   |                                                                                                                                              
aTTCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGc                                                                                                                                             >  3:580572/1‑151 (MQ=255)
   cGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAg                                                                                                                                           <  1:409083/150‑1 (MQ=255)
    gTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGtt                                                                                                                                         <  2:354979/151‑1 (MQ=255)
    gTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGtt                                                                                                                                         >  3:576114/1‑151 (MQ=255)
     tGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGtt                                                                                                                                         >  3:567013/1‑150 (MQ=255)
       cacaAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCa                                                                                                                                       >  4:387197/1‑150 (MQ=255)
         caAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAatc                                                                                                                                    <  3:397731/151‑1 (MQ=255)
         caAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAatc                                                                                                                                    >  4:117938/1‑151 (MQ=255)
            gTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAatcatc                                                                                                                                 >  3:1332543/1‑151 (MQ=255)
              tGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCg                                                                                                                               >  2:617341/1‑151 (MQ=255)
               gAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATcc                                                                                                                                <  1:114584/149‑1 (MQ=255)
               gAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCg                                                                                                                               <  3:1401487/150‑1 (MQ=255)
                aaTCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGtt                                                                                                                             <  3:998915/151‑1 (MQ=255)
                aaTCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGtt                                                                                                                             <  4:252969/151‑1 (MQ=255)
                  tCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCg                                                                                                                           >  4:282907/1‑151 (MQ=255)
                     gcgcTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCg                                                                                                                           <  3:230366/148‑1 (MQ=255)
                     gcgcTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCg                                                                                                                           <  1:50682/148‑1 (MQ=255)
                     gcgcTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGc                                                                                                                        <  1:387201/151‑1 (MQ=255)
                          gCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTg                                                                                                                   >  2:1007282/1‑151 (MQ=255)
                            aaaaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGc                                                                                                                 <  4:126174/151‑1 (MQ=255)
                              aaaGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGt                                                                                                               >  4:28982/1‑151 (MQ=255)
                                   aGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATc                                                                                                          <  3:178074/151‑1 (MQ=255)
                                    ggCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATcg                                                                                                         >  2:1085361/1‑151 (MQ=255)
                                     gCTATTCCCTGCGCGATGAACGGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATcgc                                                                                                        <  2:1143529/151‑1 (MQ=255)
                                     gCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATcgc                                                                                                        <  4:466754/151‑1 (MQ=255)
                                      cTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGTAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCCGTTCAATCATCAGATCGGGCAACTGGCGTTTATcgcg                                                                                                       >  2:1135450/1‑151 (MQ=255)
                                       tATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATcgcg                                                                                                       >  1:260371/1‑150 (MQ=255)
                                           cccTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCa                                                                                                   <  2:1042922/150‑1 (MQ=255)
                                            ccTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACg                                                                                                 >  3:138746/1‑151 (MQ=255)
                                               gcgcGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTg                                                                                              <  2:294951/151‑1 (MQ=255)
                                               gcgcGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTg                                                                                              <  3:915466/151‑1 (MQ=255)
                                                cgcgATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTgg                                                                                             >  2:217104/1‑151 (MQ=255)
                                                cgcgATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTgg                                                                                             <  2:810854/151‑1 (MQ=255)
                                                 gcgATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGt                                                                                            <  2:202984/151‑1 (MQ=255)
                                                  cgATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGtt                                                                                           >  2:92481/1‑151 (MQ=255)
                                                   gATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTg                                                                                          >  3:1277322/1‑151 (MQ=255)
                                                    aTGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTgg                                                                                         >  2:369244/1‑151 (MQ=255)
                                                      gAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGc                                                                                       >  1:21306/1‑151 (MQ=255)
                                                       aaCAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCt                                                                                      >  3:634400/1‑151 (MQ=255)
                                                         cAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGt                                                                                    <  2:1082708/151‑1 (MQ=255)
                                                              cgccgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCg                                                                               <  4:383867/151‑1 (MQ=255)
                                                              cgccgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCg                                                                               >  1:233991/1‑151 (MQ=255)
                                                                 cgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGAtt                                                                            <  4:315135/151‑1 (MQ=255)
                                                                 cgAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGAtt                                                                            <  2:1171611/151‑1 (MQ=255)
                                                                  gaagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTg                                                                           >  2:98795/1‑151 (MQ=255)
                                                                  gaagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTg                                                                           >  3:116080/1‑151 (MQ=255)
                                                                  gaagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTg                                                                           <  4:470220/151‑1 (MQ=255)
                                                                  gaagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTg                                                                           <  1:249484/151‑1 (MQ=255)
                                                                   aagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGc                                                                          >  1:335325/1‑151 (MQ=255)
                                                                   aagaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGc                                                                          >  4:335985/1‑151 (MQ=255)
                                                                    agaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCt                                                                         >  4:255291/1‑151 (MQ=255)
                                                                     gaaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTc                                                                        <  2:590654/151‑1 (MQ=255)
                                                                      aaCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCg                                                                       >  2:768716/1‑151 (MQ=255)
                                                                       aCCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGt                                                                      <  2:803696/151‑1 (MQ=255)
                                                                         cGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCa                                                                    <  2:580541/151‑1 (MQ=255)
                                                                          gCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTc                                                                     >  2:1276136/1‑149 (MQ=255)
                                                                             aGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGc                                                                 >  3:565253/1‑150 (MQ=255)
                                                                              gCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCAGGTTGGGCTGTACCCGATTGCTCGTCAGGCAt                                                               <  1:190292/151‑1 (MQ=255)
                                                                               cATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCAtt                                                              >  1:459880/1‑151 (MQ=255)
                                                                               cATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCAtt                                                              <  2:1282543/151‑1 (MQ=255)
                                                                                         gAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGAt                                                     <  4:2861/150‑1 (MQ=255)
                                                                                         gAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATc                                                    >  1:49831/1‑151 (MQ=255)
                                                                                         gAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATc                                                    >  3:1404113/1‑151 (MQ=255)
                                                                                          aaaGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCa                                                   <  3:474802/151‑1 (MQ=255)
                                                                                           aaaagaATCTTCCGCTGATTACGCTAATCGTGATGATGGCAATCATCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACGACGCTGGTTAGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCaa                                                  <  3:1135504/148‑1 (MQ=255)
                                                                                             gagaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAAt                                                >  1:68126/1‑151 (MQ=255)
                                                                                             gagaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAAt                                                <  3:1276203/151‑1 (MQ=255)
                                                                                               gaATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAAt                                                <  4:118301/149‑1 (MQ=255)
                                                                                                aaTCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCg                                             >  4:164578/1‑151 (MQ=255)
                                                                                                aaTCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCg                                             >  2:388032/1‑151 (MQ=255)
                                                                                                 aTCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCgg                                            >  3:669913/1‑151 (MQ=255)
                                                                                                      ccGCTGTTTACGCTAATCGTGATGAGGGCAACCAGTTGGGTTCTGAAGAAGTTCAATCATCCTTTCGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCa                                          <  3:1096425/147‑1 (MQ=255)
                                                                                                      ccGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCt                                       >  4:46316/1‑151 (MQ=255)
                                                                                                       cGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTa                                      <  1:141709/151‑1 (MQ=255)
                                                                                                         cTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACt                                    >  1:53164/1‑151 (MQ=255)
                                                                                                          tGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACtt                                   >  1:116057/1‑151 (MQ=255)
                                                                                                           gATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACt                                    <  1:105738/149‑1 (MQ=255)
                                                                                                           gATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTc                                  <  1:147890/151‑1 (MQ=255)
                                                                                                           gATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTc                                  <  2:1191478/151‑1 (MQ=255)
                                                                                                               aCGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGcc                               <  1:229615/150‑1 (MQ=255)
                                                                                                               aCGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCa                              <  2:634374/151‑1 (MQ=255)
                                                                                                                 gCTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCAtt                            <  2:1045082/151‑1 (MQ=255)
                                                                                                                  cTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTg                           >  1:241897/1‑151 (MQ=255)
                                                                                                                  cTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTg                           >  3:573378/1‑151 (MQ=255)
                                                                                                                  cTAATCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTg                           >  2:731167/1‑151 (MQ=255)
                                                                                                                    aaTCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGa                          <  2:1124239/150‑1 (MQ=255)
                                                                                                                      tCGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAAc                       >  2:999119/1‑151 (MQ=255)
                                                                                                                       cGTGATGATGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAAcc                      >  3:49946/1‑151 (MQ=255)
                                                                                                                          gatgatGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCtt                    <  2:484082/150‑1 (MQ=255)
                                                                                                                          gatgatGGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCtt                    >  4:393561/1‑150 (MQ=255)
                                                                                                                             gatgGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTaa                  <  4:452622/149‑1 (MQ=255)
                                                                                                                              atgGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATg                >  1:453113/1‑150 (MQ=255)
                                                                                                                              atgGCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGa               <  3:369276/151‑1 (MQ=255)
                                                                                                                                ggCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTGCCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAg              <  1:361553/150‑1 (MQ=255)
                                                                                                                                ggCAATCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAg              <  3:1007324/150‑1 (MQ=255)
                                                                                                                                 gaaaTCAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCg            >  1:256819/3‑151 (MQ=255)
                                                                                                                                      cAGCTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGcc       >  1:195770/1‑151 (MQ=255)
                                                                                                                                         cTGGGGTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCt    >  2:563084/1‑151 (MQ=255)
                                                                                                                                           ggggTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCt    <  2:1332470/149‑1 (MQ=255)
                                                                                                                                           ggggTCTGAAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTAt  <  3:617368/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ATTCGTGCACAAGTTGAATCTGCGCTGCAAAAAGCAGGCTATTCCCTGCGCGATGAACAGGCCGCCGAAGAACCGCAAGCATCACGCCTGAAAGAGAATCTGCCGCTGATTACGCTAATCGTGATGATGGCAATCAGCTGGGGTCTGGAGCAGTTCAATCATCCGTTCGGGCAACTGGCGTTTATCGCGACCACGCTGGTTGGGCTGTACCCGATTGCTCGTCAGGCATTACGGTTGATCAAATCCGGCAGCTACTTCGCCATTGAAACCTTAATGAGCGTAGCCGCTAT  >  NC_000913/3606736‑3607025

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