Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 267,774 G→A 5.5% R78C (CGT→TGT)  yafZ ← CP4‑6 prophage; conserved protein YafZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,7740GA 5.5% 227.4 / 2.2 72R78C (CGT→TGT) yafZCP4‑6 prophage; conserved protein YafZ
Reads supporting (aligned to +/- strand):  ref base G (37/31);  new base A (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CCCGTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATACACTGGGCAC  >  NC_000913/267625‑267921
                                                                                                                                                     |                                                                                                                                                   
cccGTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCacg                                                                                                                                                     >  2:444947/1‑150 (MQ=255)
 ccGTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCacga                                                                                                                                                    <  2:845887/150‑1 (MQ=255)
 ccGTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGATTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCacga                                                                                                                                                    <  1:299063/150‑1 (MQ=255)
  cGTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCacgac                                                                                                                                                   >  2:799761/1‑150 (MQ=255)
   gTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCacgacg                                                                                                                                                  >  1:1146648/1‑150 (MQ=255)
    ttCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCAcgacga                                                                                                                                                 <  2:1147203/150‑1 (MQ=255)
      cTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGAcc                                                                                                                                               <  1:1004803/150‑1 (MQ=255)
      cTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGAcc                                                                                                                                               <  2:1326430/150‑1 (MQ=255)
        gACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCg                                                                                                                                             >  2:1388146/1‑150 (MQ=255)
        gACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCg                                                                                                                                             >  1:51315/1‑150 (MQ=255)
            aaCCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTc                                                                                                                                         <  1:1262108/150‑1 (MQ=255)
               cGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTcacg                                                                                                                                      >  1:1206622/1‑150 (MQ=255)
                     aaaTAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCAcgcg                                                                                                                                >  1:1132975/1‑150 (MQ=255)
                       aTAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCggg                                                                                                                              >  2:1530412/1‑150 (MQ=255)
                         aGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATTTGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTc                                                                                                                            <  1:433795/150‑1 (MQ=255)
                         aGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCGGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTc                                                                                                                            >  2:468287/1‑150 (MQ=255)
                           gcccGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTg                                                                                                                          <  1:1479129/149‑1 (MQ=255)
                             ccGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGAc                                                                                                                        >  2:1550134/1‑150 (MQ=255)
                             ccGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGAc                                                                                                                        >  1:50389/1‑150 (MQ=255)
                                   aCATCTGATACGAACTGGTTCCATCGTGAGAGATGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCaa                                                                                                                  <  1:1618977/150‑1 (MQ=255)
                                      tCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGa                                                                                                               <  2:302632/150‑1 (MQ=255)
                                         gATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATg                                                                                                            >  2:1141667/1‑150 (MQ=255)
                                         gATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCAGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATg                                                                                                            >  2:70026/1‑150 (MQ=255)
                                           tACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGc                                                                                                          >  1:1447841/1‑150 (MQ=255)
                                             cGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTg                                                                                                        <  1:1051863/150‑1 (MQ=255)
                                             cGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTg                                                                                                        >  1:966983/1‑150 (MQ=255)
                                               aaCTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGa                                                                                                      <  1:1105225/150‑1 (MQ=255)
                                                 cTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAg                                                                                                    <  1:799761/150‑1 (MQ=255)
                                                   ggTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGcc                                                                                                  >  2:1038649/1‑150 (MQ=255)
                                                   ggTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGcc                                                                                                  >  2:3954/1‑150 (MQ=255)
                                                        cATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTtctc                                                                                             <  1:1550134/150‑1 (MQ=255)
                                                              gagaGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTa                                                                                       <  2:1206622/150‑1 (MQ=255)
                                                                  gTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCt                                                                                   >  2:641901/1‑150 (MQ=255)
                                                                  cttGAGTAGAATAAATTCCTGAACCTGTTTACCGGTGATCTGACATTCCCGCCGCAGACGCAGCATATTATTTGTATCTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCt                                                                                   >  2:698814/2‑150 (MQ=21)
                                                                    tGAGAAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCACGCCGCAGACGCAGCATATGATTAGTATGTTCACGACGACCCGGGTCACGCCCGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTATCGCTGTAGCCTGt                                                                                 >  1:1595091/1‑150 (MQ=21)
                                                                     gAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTc                                                                                >  1:970961/1‑150 (MQ=255)
                                                                     gAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACAACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGta                                                                                >  2:768071/1‑149 (MQ=255)
                                                                     gAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACAACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTc                                                                                >  2:773896/1‑150 (MQ=255)
                                                                          gAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGACTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCa                                                                           >  2:94798/1‑150 (MQ=255)
                                                                           aataatTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAg                                                                          <  1:476313/150‑1 (MQ=255)
                                                                              aatTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGga                                                                       <  1:153276/150‑1 (MQ=255)
                                                                               atTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGgag                                                                      >  1:263606/1‑150 (MQ=255)
                                                                                tttCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGgaga                                                                     >  1:1456364/1‑150 (MQ=255)
                                                                                tttCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGgaga                                                                     <  2:1190912/150‑1 (MQ=255)
                                                                                   ccGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGAtgg                                                                  <  1:1388146/150‑1 (MQ=255)
                                                                                    cGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGAtggt                                                                 >  1:1605333/1‑150 (MQ=255)
                                                                                     ggCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGAtggtg                                                                >  1:1544232/1‑150 (MQ=255)
                                                                                      gCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGAtggtgg                                                               <  2:741909/150‑1 (MQ=255)
                                                                                       cagcTGTTTACCGGGGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTggg                                                              <  2:485316/147‑1 (MQ=255)
                                                                                           tGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGtata                                                          >  1:817938/1‑150 (MQ=255)
                                                                                           tGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGCGAGATGGTGGGtata                                                          >  2:490986/1‑150 (MQ=255)
                                                                                            gTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGtatat                                                         <  1:857882/150‑1 (MQ=255)
                                                                                             tttACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGtatata                                                        >  2:173275/1‑150 (MQ=255)
                                                                                              ttACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAg                                                       <  1:1483436/150‑1 (MQ=255)
                                                                                              ttACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAg                                                       >  2:859018/1‑150 (MQ=255)
                                                                                               tACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAgg                                                      <  1:1530412/150‑1 (MQ=255)
                                                                                                 ccGGGGATCTGCCCGTCCCGCCGCAGACGCAGCATATGCTTTGTATGTGCACGACGACCCGGGGCACGCACTCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGta                                                    <  2:871768/150‑1 (MQ=21)
                                                                                                      gATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACg                                               >  1:24523/1‑150 (MQ=255)
                                                                                                         cTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTc                                            >  2:244906/1‑150 (MQ=255)
                                                                                                           gCCCTTCCCGCCGCAGCCGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCCCGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCAc                                          >  2:1387549/1‑150 (MQ=255)
                                                                                                           gCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCAc                                          <  2:1550617/150‑1 (MQ=255)
                                                                                                            cccTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACt                                         <  2:51315/150‑1 (MQ=255)
                                                                                                             ccTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTa                                        >  1:593074/1‑150 (MQ=255)
                                                                                                                tCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACgg                                     >  1:458499/1‑150 (MQ=255)
                                                                                                                      cgcAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCg                               <  1:1342946/150‑1 (MQ=255)
                                                                                                                       gcAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCgt                              >  2:883200/1‑150 (MQ=255)
                                                                                                                        cAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCgtg                             <  1:3954/150‑1 (MQ=255)
                                                                                                                          gACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGtt                           <  2:1447841/150‑1 (MQ=255)
                                                                                                                            cGCAGCATATGCTTTGTATGTTCACAACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTg                         <  1:768071/150‑1 (MQ=255)
                                                                                                                            cGCAGCATATGCTTTGTATGTTCACAACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTg                         <  1:773896/150‑1 (MQ=255)
                                                                                                                             gcagcaTATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGt                        <  1:444947/150‑1 (MQ=255)
                                                                                                                                    atGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACt                 <  2:24523/150‑1 (MQ=255)
                                                                                                                                      gCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGa               >  2:570600/1‑150 (MQ=255)
                                                                                                                                         ttGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATa            >  1:1028935/1‑150 (MQ=255)
                                                                                                                                          tGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATac           >  1:505526/1‑150 (MQ=255)
                                                                                                                                          tGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATac           <  1:70026/150‑1 (MQ=255)
                                                                                                                                                 tCACGACGACCCGTTTCACGCACTCGGGTCTGACAGGCAAAGAATGGCTGGAAGACTTCTCGCTGTAGGCTGTCGATCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTTAATACACTGGGc    <  1:1470135/150‑1 (MQ=21)
                                                                                                                                                 tCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATACACTGGGc    >  2:1112856/1‑150 (MQ=255)
                                                                                                                                                   acgacgACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATACACTGGGCAc  <  1:659121/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCGTTCTGACAAACCGCACGAAATAGTCCCGGCAACATCTGATACGAACTGGTTCCATCGTGAGAGTTGAGTAGAATAATTTCCGGCACCTGTTTACCGGTGATCTGCCCTTCCCGCCGCAGACGCAGCATATGCTTTGTATGTTCACGACGACCCGGGTCACGCACGCGGGTCTGACAGGCAAAGAATGGCTGGAAGCCTTCTCGCTGTAGGCTGTCGAGCAGGGAGATGGTGGGTATATAGGTATAACGCTCACTACGGGACTCGTGTTTGTCCTCACTGAATACACTGGGCAC  >  NC_000913/267625‑267921

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