Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,579,672 C→T 6.1% A57T (GCG→ACG)  symE ← toxic peptide regulated by antisense sRNA symR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,579,6720CT 6.1% 222.0 / 3.8 66A57T (GCG→ACG) symEtoxic peptide regulated by antisense sRNA symR
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base T (2/2);  total (31/35)
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.98e-01

CACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAA  >  NC_000913/4579526‑4579821
                                                                                                                                                  |                                                                                                                                                     
cACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGt                                                                                                                                                    <  1:942317/150‑1 (MQ=255)
cACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGt                                                                                                                                                    <  2:938192/150‑1 (MQ=255)
  cTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGc                                                                                                                                                  >  1:444507/1‑150 (MQ=255)
   tCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGcc                                                                                                                                                 >  1:98029/1‑150 (MQ=255)
   tCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGcc                                                                                                                                                 >  1:688896/1‑150 (MQ=255)
   tCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGcc                                                                                                                                                 >  1:680383/1‑150 (MQ=255)
   tCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGcc                                                                                                                                                 >  1:495653/1‑150 (MQ=255)
      aaTAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGTCGTGCCAGt                                                                                                                                              >  1:892224/1‑150 (MQ=255)
      aaTAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGTCGTGCCAGt                                                                                                                                              >  1:770761/1‑150 (MQ=255)
       aTAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTg                                                                                                                                             <  2:607561/150‑1 (MQ=255)
       aTAAACGCCTGCACCTGCTTTTGATTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGATGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTg                                                                                                                                             <  2:1226833/150‑1 (MQ=255)
           aCGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCaa                                                                                                                                         <  2:558066/150‑1 (MQ=255)
               cTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAAcc                                                                                                                                     <  2:1020727/150‑1 (MQ=255)
                   aCCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACcggcg                                                                                                                                 >  1:319688/1‑150 (MQ=255)
                     cTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCggcggc                                                                                                                               >  1:1019843/1‑150 (MQ=255)
                        cTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTc                                                                                                                            >  2:552333/1‑150 (MQ=255)
                         ttttGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGACATTGGCAAAACCGGCGGCTTcc                                                                                                                           <  2:834896/150‑1 (MQ=255)
                                      ccGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCtt                                                                                                              >  1:1082489/1‑150 (MQ=255)
                                        gACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTCCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTc                                                                                                            >  2:103501/1‑150 (MQ=255)
                                         aCAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGTCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCa                                                                                                           >  1:267448/1‑150 (MQ=255)
                                           aGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAgg                                                                                                         >  1:1261027/1‑150 (MQ=255)
                                           aGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAgg                                                                                                         <  2:87549/150‑1 (MQ=255)
                                             tttGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGt                                                                                                       <  1:549842/150‑1 (MQ=255)
                                                      cTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCggg                                                                                              <  1:1108563/150‑1 (MQ=255)
                                                        ggCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGaa                                                                                            <  2:152121/150‑1 (MQ=255)
                                                             aGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACgg                                                                                       <  2:1366628/150‑1 (MQ=255)
                                                             aGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACgg                                                                                       <  1:160362/150‑1 (MQ=255)
                                                             aGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACgg                                                                                       <  1:311102/150‑1 (MQ=255)
                                                             aGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTCGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACgg                                                                                       <  1:1016002/150‑1 (MQ=255)
                                                               cGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGTCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCt                                                                                     <  2:892224/150‑1 (MQ=255)
                                                               cGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGTCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCt                                                                                     <  2:770761/150‑1 (MQ=255)
                                                               cGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCt                                                                                     <  1:955397/150‑1 (MQ=255)
                                                                gACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTg                                                                                    >  1:504024/1‑150 (MQ=255)
                                                                     cATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATc                                                                               >  2:311594/1‑150 (MQ=255)
                                                                      aTCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATcc                                                                              >  1:490713/1‑150 (MQ=255)
                                                                       tCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCg                                                                             <  2:653220/150‑1 (MQ=255)
                                                                            tCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAg                                                                        >  2:528274/1‑150 (MQ=255)
                                                                            tCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAg                                                                        <  2:38973/150‑1 (MQ=255)
                                                                               ctctCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGa                                                                     <  2:649487/150‑1 (MQ=255)
                                                                                tctcCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGACTACGGCTGTAATCCGGGTAGCGAc                                                                    >  2:1299420/1‑150 (MQ=255)
                                                                                  tcctcGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTACTCCGGGTAGCGACTc                                                                  >  1:335258/1‑150 (MQ=255)
                                                                                     tcGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCa                                                               >  2:1283645/1‑150 (MQ=255)
                                                                                     tcGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCa                                                               >  2:940386/1‑150 (MQ=255)
                                                                                     tcGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCa                                                               >  2:304196/1‑150 (MQ=255)
                                                                                     tcGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCa                                                               >  2:940385/1‑150 (MQ=255)
                                                                                                  gtTGGGGGGTGAGCAAAATACAGCCTTCCATCATTTTTACATCTCCCGCGGTGCCTGTGCAAAAACCGGCGGCTTCCAGCCAGTGACCTTACTGGGTGTTGGCGGTAATACGGCTTTAATCGGGTTAGCGACTCGCATAACCGACGGTGa                                                  <  1:1235691/150‑1 (MQ=255)
                                                                                                    tGGGCGGTGAGTAAAATACAGCCTTCCATAACTTTGACATCTACCGCCGTGCCTGTGGCAAACCGGGCGGCTTCCAGCCAGTGACCTTTCTGGGTTATGGCGGTATTACGGCTTTATTCCGGGTAGCGACTCGCATAACCGACGGTGACa                                                <  1:1250982/150‑1 (MQ=255)
                                                                                                         ggTGTGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCTGTGGCAAAACCGGCGGCTACCAGCCACTGACCTTTCTGGGTGTTGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACg                                           <  1:910060/150‑1 (MQ=255)
                                                                                                         ggTGAGGACAATACCGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGTAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACg                                           <  1:910049/150‑1 (MQ=255)
                                                                                                                    tACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCg                                >  1:848889/1‑150 (MQ=255)
                                                                                                                       aGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCgggg                             <  2:1078843/150‑1 (MQ=255)
                                                                                                                         ccTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGgag                           >  2:1090704/1‑150 (MQ=255)
                                                                                                                          cTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGgaga                          <  2:916388/150‑1 (MQ=255)
                                                                                                                              cATCACTTTGACATCCACCGCCGTGCCAGTGGCAAACCCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTc                      <  1:340920/150‑1 (MQ=255)
                                                                                                                                tCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTg                    >  1:598496/1‑150 (MQ=255)
                                                                                                                                   cTTTGCCATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCCGGGTGAGGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATACCCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCtt                 <  1:822378/150‑1 (MQ=255)
                                                                                                                                    tttGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTc                <  2:1396290/150‑1 (MQ=255)
                                                                                                                                    tttGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTc                >  1:502038/1‑150 (MQ=255)
                                                                                                                                    tttGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCAACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTc                <  1:1210946/150‑1 (MQ=255)
                                                                                                                                     ttGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCg               >  1:1248961/1‑150 (MQ=255)
                                                                                                                                      tGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACCTGACGGTGATTTGCCGGGGAGACTTCTGCGTCGa              >  1:625438/1‑150 (MQ=255)
                                                                                                                                         cATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTTATGGCTGTAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACg           <  1:862583/150‑1 (MQ=255)
                                                                                                                                         cATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACg           <  2:341001/150‑1 (MQ=255)
                                                                                                                                          aTCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACgg          <  1:155381/150‑1 (MQ=255)
                                                                                                                                             tACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACGGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTtg       >  1:358011/1‑150 (MQ=255)
                                                                                                                                              aCCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTtgt      <  1:336182/150‑1 (MQ=255)
                                                                                                                                               ccgccgTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTtgtg     <  2:1139613/150‑1 (MQ=255)
                                                                                                                                               ccgccgTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTtgtg     >  2:1358918/1‑150 (MQ=255)
                                                                                                                                                cgccgTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGGTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGc    >  1:1119645/1‑150 (MQ=255)
                                                                                                                                                  ccgTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACCGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTGATTTGCCGGGGAGACTTCTGCGTCGAACGGTTGTGCaa  >  1:824493/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCGACTGCATCAGTTCGCTCTCCTCGGCGGCGGGTGGTTGGGCGGTGAGGACAATACAGCCTTCCATCACTTTGACATCTACCGCCGTGCCAGTGGCAAAACCGGCGGCTTCCAGCCACTGACCTTTCAGGGTGATGGCGGGAATACGGCTGTAATCCGGGTAGCGACTCGCATAACCGACGGTGACATGACGGTTATTTGCCGGGGAGACTTCTGCTTCGAACGGTTGTGCAA  >  NC_000913/4579526‑4579821

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