Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 784,012 T→A 5.4% Q271L (CAG→CTG)  zitB ← Zn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913784,0120TA 5.4% 303.9 / 5.2 93Q271L (CAG→CTG) zitBZn(2(+))/Cd(2(+))/Ni(2(+))/Cu(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base T (44/44);  new base A (3/2);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

CTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATA  >  NC_000913/783864‑784161
                                                                                                                                                    |                                                                                                                                                     
cTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                      >  2:236991/1‑150 (MQ=255)
  cACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCATGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGAt                                                                                                                                                    <  2:917547/150‑1 (MQ=255)
  cACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGAt                                                                                                                                                    <  1:1370583/150‑1 (MQ=255)
  cACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGAt                                                                                                                                                    <  2:1996546/150‑1 (MQ=255)
    cgcgcgCTGTCCCATTTATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATcc                                                                                                                                                  >  2:852995/1‑150 (MQ=255)
    cgcgcgCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATcc                                                                                                                                                  <  2:241019/150‑1 (MQ=255)
       gcgcTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaac                                                                                                                                               <  1:1944727/150‑1 (MQ=255)
         gcTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaacaa                                                                                                                                             <  1:852995/150‑1 (MQ=255)
             tCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCa                                                                                                                                         >  2:890216/1‑150 (MQ=255)
              cccATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCat                                                                                                                                        >  2:875692/1‑150 (MQ=255)
              cccATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacct                                                                                                                                        <  2:1941890/150‑7 (MQ=255)
              cccATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacc                                                                                                                                         >  1:350488/1‑144 (MQ=255)
              cccATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCAGATCcaacactacc                                                                                                                                         <  2:350488/149‑6 (MQ=255)
                cATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcg                                                                                                                                      <  1:780077/150‑1 (MQ=255)
                cATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcg                                                                                                                                      >  1:1971124/1‑150 (MQ=255)
                      tGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcg                                                                                                                                >  1:1780475/1‑150 (MQ=255)
                      tGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcg                                                                                                                                >  1:125273/1‑150 (MQ=255)
                        gTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcgtg                                                                                                                              >  2:865520/1‑150 (MQ=255)
                         tgatgaTGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCgtgt                                                                                                                             >  1:1736561/1‑150 (MQ=255)
                         tgatgaTGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCgtgt                                                                                                                             >  2:1441832/1‑150 (MQ=255)
                            tgatgTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGc                                                                                                                          >  2:1779223/1‑150 (MQ=255)
                                gtgAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGa                                                                                                                      <  1:875692/150‑1 (MQ=255)
                                  gAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATc                                                                                                                    <  2:81939/150‑1 (MQ=255)
                                  gAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATc                                                                                                                    <  1:1478304/150‑1 (MQ=255)
                                  gAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATc                                                                                                                    <  1:927205/150‑1 (MQ=255)
                                   aaTGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                                                                   >  2:624366/1‑150 (MQ=255)
                                    aTGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAc                                                                                                                  <  2:1529369/150‑1 (MQ=255)
                                       gCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTg                                                                                                               >  2:687909/1‑150 (MQ=255)
                                         cAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGCGATCGTGTGGCGGGATCACCTGca                                                                                                             <  2:1470477/150‑1 (MQ=255)
                                             cacCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATg                                                                                                         <  1:73187/150‑1 (MQ=255)
                                             cacCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATg                                                                                                         <  2:1408731/150‑1 (MQ=255)
                                             cacCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATg                                                                                                         <  1:1953255/150‑1 (MQ=255)
                                               ccccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTTATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                       <  1:1208111/150‑1 (MQ=255)
                                               ccccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                       >  1:714390/1‑150 (MQ=255)
                                               ccccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                       <  1:1046141/150‑1 (MQ=255)
                                                 ccTTCATTAAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTACATCTGAATGGTGGCGTGCGCAATATGATAGTGATCCATAAGGTAGTGTTGGATATGATACAACAAGGCATCGTTATCGTGTGGCGGGATCACCTGCACATGCCAc                                                                                                     >  2:1011257/1‑150 (MQ=255)
                                                     cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtca                                                                                                 >  1:34768/1‑150 (MQ=255)
                                                     cattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCCCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtca                                                                                                 <  1:460908/150‑1 (MQ=255)
                                                         caGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCAc                                                                                             >  1:1735627/1‑150 (MQ=255)
                                                           ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACACGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCCCCg                                                                                           >  2:1066665/1‑150 (MQ=255)
                                                           ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCg                                                                                           <  2:969212/150‑1 (MQ=255)
                                                           ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCg                                                                                           >  2:117030/1‑150 (MQ=255)
                                                              ggCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCt                                                                                        >  2:250571/1‑150 (MQ=255)
                                                               gCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCtt                                                                                       >  1:800312/1‑150 (MQ=255)
                                                                   tCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCg                                                                                   <  1:1897301/150‑1 (MQ=255)
                                                                   tCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCg                                                                                   <  1:1779223/150‑1 (MQ=255)
                                                                    cTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                                                                  <  1:564994/150‑1 (MQ=255)
                                                                     tGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGcc                                                                                 <  1:253827/150‑1 (MQ=255)
                                                                     tGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGcc                                                                                 <  1:954226/150‑1 (MQ=255)
                                                                      ggCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCt                                                                                >  2:1679160/1‑150 (MQ=255)
                                                                       gCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTa                                                                               <  2:1896017/150‑1 (MQ=255)
                                                                        cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                                                                              >  2:1182064/1‑150 (MQ=255)
                                                                        cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                                                                              >  2:1190027/1‑150 (MQ=255)
                                                                         ccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAcc                                                                             <  1:250571/150‑1 (MQ=255)
                                                                            tGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                                                                          >  1:556004/1‑150 (MQ=255)
                                                                             gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcc                                                                         >  1:1604989/1‑150 (MQ=255)
                                                                             gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcc                                                                         <  1:624366/150‑1 (MQ=255)
                                                                                  ggTTGATATTCCATCTGAATGGGGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtaca                                                                    <  2:1889203/150‑1 (MQ=255)
                                                                                     tGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgt                                                                 >  2:341991/1‑150 (MQ=255)
                                                                                      gATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgta                                                                <  1:393574/150‑1 (MQ=255)
                                                                                          ttCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacatg                                                            >  2:1421525/1‑150 (MQ=255)
                                                                                          ttCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacatg                                                            >  2:1518803/1‑150 (MQ=255)
                                                                                           tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATgg                                                           <  1:1617370/150‑1 (MQ=255)
                                                                                            ccATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGt                                                          <  1:688075/150‑1 (MQ=255)
                                                                                            ccATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGt                                                          >  2:315764/1‑150 (MQ=255)
                                                                                               tCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGca                                                       <  1:236991/150‑1 (MQ=255)
                                                                                                cTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGcac                                                      >  2:423779/1‑150 (MQ=255)
                                                                                                     tgggggCGTGCTCAATCTGATAGTGATCCATCAGGGAGTGTTGGATCTGATCCAACAAGGCATCGTGCTCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGc                                                 <  2:1694276/150‑1 (MQ=255)
                                                                                                       gtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCgg                                               <  2:1736561/150‑1 (MQ=255)
                                                                                                       gtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCgg                                               <  1:1679160/150‑1 (MQ=255)
                                                                                                         ggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGAc                                             >  2:1283057/1‑150 (MQ=255)
                                                                                                           cGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACtt                                           >  1:707994/1‑150 (MQ=255)
                                                                                                              gCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCg                                        >  2:1983687/1‑150 (MQ=255)
                                                                                                               cTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCgg                                       >  2:1940854/1‑150 (MQ=255)
                                                                                                                 cAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGa                                     >  1:1088850/1‑150 (MQ=255)
                                                                                                                        aTAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCg                              <  1:1399012/150‑1 (MQ=255)
                                                                                                                         tAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCgg                             >  2:1030102/1‑150 (MQ=255)
                                                                                                                          aGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGc                            <  1:1011257/150‑1 (MQ=255)
                                                                                                                           gTGATCCATCAGGTAGTTTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATAACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCAATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGca                           >  2:144675/1‑150 (MQ=255)
                                                                                                                           gTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGca                           <  1:865520/150‑1 (MQ=255)
                                                                                                                                 cATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacg                     >  1:895966/1‑150 (MQ=255)
                                                                                                                                 cATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacg                     <  2:1735627/150‑1 (MQ=255)
                                                                                                                                        tAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttc              <  2:218583/150‑1 (MQ=255)
                                                                                                                                         cttgtTGGATCAGCTCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACCTGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca             >  2:1615047/3‑150 (MQ=255)
                                                                                                                                         cttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttca             >  1:764538/3‑150 (MQ=255)
                                                                                                                                            gtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttcagtt          <  2:125273/150‑1 (MQ=255)
                                                                                                                                            gtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttcagtt          >  2:1723875/1‑150 (MQ=255)
                                                                                                                                            gtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttcagtt          <  2:1780475/150‑1 (MQ=255)
                                                                                                                                                gATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGc      <  1:1441832/150‑1 (MQ=255)
                                                                                                                                                 aTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCg     >  1:507760/1‑150 (MQ=255)
                                                                                                                                                 aTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCg     >  1:831000/1‑150 (MQ=255)
                                                                                                                                                 aTCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCg     >  1:1357302/1‑150 (MQ=255)
                                                                                                                                                  tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGa    >  2:185751/1‑150 (MQ=255)
                                                                                                                                                    tGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGTATTTCCCGGCACATACGACTCTTCAGTTCCGCGata  >  2:1164402/1‑150 (MQ=255)
                                                                                                                                                    tGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGata  >  1:1564178/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCGATA  >  NC_000913/783864‑784161

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