Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,705 C→T 9.3% intergenic (+42/+18) djlC → / ← hscC co‑chaperone DjlC/chaperone protein HscC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,7050CT 9.3% 216.0 / 10.0 75intergenic (+42/+18)djlC/hscCco‑chaperone DjlC/chaperone protein HscC
Reads supporting (aligned to +/- strand):  ref base C (32/36);  new base T (5/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

AATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAG  >  NC_000913/681556‑681854
                                                                                                                                                     |                                                                                                                                                     
aaTGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                                                                                       <  2:717405/150‑1 (MQ=255)
 atgatgAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGCGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccg                                                                                                                                                      <  1:424392/150‑1 (MQ=255)
    atgaGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCAt                                                                                                                                                   >  1:1349377/1‑150 (MQ=255)
      gaGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATcc                                                                                                                                                 <  1:306821/150‑1 (MQ=255)
       aGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCACATCCg                                                                                                                                                >  1:421832/1‑150 (MQ=255)
        gggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGAGTTTACGACGCATCCgg                                                                                                                                               >  2:1963542/1‑150 (MQ=255)
         ggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGc                                                                                                                                              <  1:1886719/150‑1 (MQ=255)
           tGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAg                                                                                                                                            <  2:303312/150‑1 (MQ=255)
            gAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCGGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGc                                                                                                                                           >  1:1939737/1‑150 (MQ=255)
                aaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtg                                                                                                                                       <  1:545519/150‑1 (MQ=255)
                aaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCGGACGATGGGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCatcggcg                                                                                                                                       <  1:1261335/150‑7 (MQ=38)
                aaCCACAGAGCATTACCAGAATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCTTTTACGCCGCATCCGGCAGCCgtg                                                                                                                                       <  2:1581127/150‑1 (MQ=255)
                 aCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                                                                                      <  2:276589/150‑1 (MQ=255)
                          cATTACCAGCATTATCAAAATCAGTTTAACGTTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                                                                             <  1:948943/150‑1 (MQ=255)
                             tACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtctc                                                                                                                          >  1:297266/1‑150 (MQ=255)
                                cAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                                                                                       <  1:632644/150‑1 (MQ=255)
                                 aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGTATGCAGTGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                                      <  2:1139550/150‑1 (MQ=255)
                                 aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                                      <  2:922440/150‑1 (MQ=255)
                                   cATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCGTCCGGCAGCCGTGTTTATGGGATCTCAATGGc                                                                                                                    <  1:346700/150‑1 (MQ=255)
                                     ttATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTa                                                                                                                  <  2:1043743/150‑1 (MQ=255)
                                     ttATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTa                                                                                                                  <  1:1957268/150‑1 (MQ=255)
                                      tATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaa                                                                                                                 >  1:1289981/1‑150 (MQ=255)
                                       aTCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaaa                                                                                                                <  1:1769998/150‑1 (MQ=255)
                                            aaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                                                                                           <  2:66498/150‑1 (MQ=255)
                                             aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcg                                                                                                          >  2:1462652/1‑150 (MQ=255)
                                                 gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcg                                                                                                      >  2:894842/1‑150 (MQ=255)
                                                 gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcg                                                                                                      >  2:467938/1‑150 (MQ=255)
                                                  tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtccctcgc                                                                                                     >  1:1024706/1‑150 (MQ=255)
                                                    tAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                                                                                   <  1:1610342/150‑1 (MQ=255)
                                                    tAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                                                                                   >  2:1404943/1‑150 (MQ=255)
                                                    tAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                                                                                   >  1:1659751/1‑150 (MQ=255)
                                                      aCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgc                                                                                                 <  1:523122/150‑1 (MQ=255)
                                                           gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgccg                                                                                            <  2:651804/150‑1 (MQ=255)
                                                              ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                                                                                         <  1:335180/150‑1 (MQ=255)
                                                                tGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                                                                                       >  2:1797864/1‑150 (MQ=255)
                                                                 gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTc                                                                                      >  1:1958171/1‑150 (MQ=255)
                                                                     ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCGGCCGCCGAACTTCATcc                                                                                  <  2:1939737/150‑1 (MQ=255)
                                                                         aTGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgc                                                                              <  1:169840/150‑1 (MQ=255)
                                                                            ctctCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                           >  1:1024789/1‑150 (MQ=255)
                                                                             tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGc                                                                          >  2:642911/1‑150 (MQ=255)
                                                                               tcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa                                                                        >  2:129769/1‑150 (MQ=255)
                                                                                cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat                                                                       >  1:822986/1‑150 (MQ=255)
                                                                                       aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt                                                                <  2:632446/150‑1 (MQ=255)
                                                                                       aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt                                                                <  2:1349377/150‑1 (MQ=255)
                                                                                        tCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg                                                               >  2:31059/1‑150 (MQ=255)
                                                                                          gcgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCGTCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGAc                                                             <  2:492639/150‑1 (MQ=255)
                                                                                          gcgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGAc                                                             <  1:1060108/150‑1 (MQ=255)
                                                                                           cgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACt                                                            >  1:803405/1‑150 (MQ=255)
                                                                                                      aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCa                                                 >  2:747539/1‑150 (MQ=255)
                                                                                                      aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCa                                                 <  1:1137390/150‑1 (MQ=255)
                                                                                                      aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCa                                                 >  2:129891/1‑150 (MQ=255)
                                                                                                       gCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCAt                                                >  2:153830/1‑150 (MQ=255)
                                                                                                        cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATc                                               >  1:240050/1‑150 (MQ=255)
                                                                                                        cTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGGTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATc                                               <  1:413019/150‑1 (MQ=255)
                                                                                                         tGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCa                                              >  1:16289/1‑150 (MQ=255)
                                                                                                                gtgtTTTTTGCGCCGATGCCGGATGCAGCGGTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCgg                                       <  1:395908/150‑1 (MQ=255)
                                                                                                                 tgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGt                                      <  1:1612530/150‑1 (MQ=255)
                                                                                                                  gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTg                                     <  1:1797864/150‑1 (MQ=255)
                                                                                                                  gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTg                                     >  2:220083/1‑150 (MQ=255)
                                                                                                                  gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTg                                     <  1:1660379/150‑1 (MQ=255)
                                                                                                                   ttttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGa                                    >  1:1465744/1‑150 (MQ=255)
                                                                                                                    tttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGAt                                   >  1:1644676/1‑150 (MQ=255)
                                                                                                                          cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGGTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCg                             <  2:1258919/150‑1 (MQ=255)
                                                                                                                          cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCg                             >  1:1822844/1‑150 (MQ=255)
                                                                                                                          cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCg                             >  2:1428621/1‑150 (MQ=255)
                                                                                                                          cgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCg                             >  1:1617766/1‑150 (MQ=255)
                                                                                                                            ccGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGc                           >  1:317482/1‑150 (MQ=255)
                                                                                                                             cGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGcc                          <  1:467938/150‑1 (MQ=255)
                                                                                                                              gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCg                         <  1:129891/150‑1 (MQ=255)
                                                                                                                                 gCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCCCCCGGCCGAtg                      >  1:620078/1‑149 (MQ=255)
                                                                                                                                 gCCGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTAAGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtc                      >  2:502146/1‑150 (MQ=255)
                                                                                                                                 gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtc                      <  1:220083/150‑1 (MQ=255)
                                                                                                                                  ccGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtct                     >  2:1316544/1‑150 (MQ=255)
                                                                                                                                  ccGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtct                     >  2:1309733/1‑150 (MQ=255)
                                                                                                                                  ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtct                     <  2:1465744/150‑1 (MQ=255)
                                                                                                                                   cGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtctc                    >  1:334065/1‑150 (MQ=255)
                                                                                                                                                tACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCa       <  1:747539/150‑1 (MQ=255)
                                                                                                                                                  cgccgcATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCAcc     >  1:624946/1‑150 (MQ=255)
                                                                                                                                                   gccgcATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCAccc    >  2:1361854/1‑150 (MQ=255)
                                                                                                                                                   gccgcATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCAccc    >  2:208358/1‑150 (MQ=255)
                                                                                                                                                    ccgcATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCa   <  2:1958171/150‑1 (MQ=255)
                                                                                                                                                    ccgcATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCa   <  1:642911/150‑1 (MQ=255)
                                                                                                                                                     cgcATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAg  >  2:221500/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAG  >  NC_000913/681556‑681854

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: