Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 572,368 A→G M1M (ATG→GTG) † ylcH → DLP12 prophage; uncharacterized protein YlcH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913572,3680AG88.7% 205.0 / 13.4 71M1M (ATG→GTG) †ylcHDLP12 prophage; uncharacterized protein YlcH
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (32/31);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

AAAGCTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATC  >  NC_000913/572234‑572499
                                                                                                                                      |                                                                                                                                   
aaaGCTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTgcgc                                                                                                                                >  2:924676/1‑140 (MQ=255)
aaaGCTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTgcgc                                                                                                                                >  2:30673/1‑140 (MQ=255)
aaaGCTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTCTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTgcgc                                                                                                                                >  2:282222/1‑140 (MQ=255)
  aGCTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGcac                                                                                                                              >  2:364232/1‑140 (MQ=255)
    cTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACAt                                                                                                                            >  2:848351/1‑140 (MQ=255)
     ttCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCTCCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACAtg                                                                                                                           <  2:115029/140‑2 (MQ=255)
      tCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATAc                                                                                                                          <  2:62446/140‑1 (MQ=255)
        cAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACaa                                                                                                                        <  2:643857/140‑1 (MQ=255)
         agagagTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAAt                                                                                                                       >  2:777265/1‑140 (MQ=255)
               tAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAAc                                                                                                                 <  2:589580/140‑1 (MQ=255)
               tAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAAct                                                                                                                >  2:747676/1‑140 (MQ=255)
                    gtgtTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAActctct                                                                                                            <  2:394752/140‑1 (MQ=255)
                     tgtTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGGCTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCtctt                                                                                                           <  2:352535/140‑1 (MQ=255)
                     tgtTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCtctt                                                                                                           <  1:364232/140‑1 (MQ=255)
                      gtTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGTACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCtcttc                                                                                                          <  1:927689/140‑1 (MQ=255)
                       ttttAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAActctct                                                                                                            >  2:141018/1‑137 (MQ=255)
                       ttttAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAActctct                                                                                                            <  1:141018/137‑1 (MQ=255)
                       ttttAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCtcttct                                                                                                         >  1:380947/1‑140 (MQ=255)
                       ttttAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCtcttct                                                                                                         >  2:610446/1‑140 (MQ=255)
                                 gAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAg                                                                                               >  2:335645/1‑140 (MQ=255)
                                 gAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAg                                                                                               >  2:240383/1‑140 (MQ=255)
                                           aTCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAAt                                                                                     >  2:147968/1‑140 (MQ=255)
                                             cGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGc                                                                                   <  2:380947/140‑1 (MQ=255)
                                                 tttgattttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCt                                                                               >  2:808163/1‑140 (MQ=255)
                                                    gattttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCa                                                                            <  1:243064/140‑1 (MQ=255)
                                                    gattttgGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCa                                                                            <  1:610446/140‑1 (MQ=255)
                                                          gggAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCtt                                                                      <  1:848351/140‑1 (MQ=255)
                                                           ggAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCttt                                                                     <  2:459347/140‑1 (MQ=255)
                                                               aTCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCAt                                                                 >  1:25528/1‑140 (MQ=255)
                                                                tCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATc                                                                >  2:430920/1‑140 (MQ=255)
                                                                         cATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTg                                                       >  1:344609/1‑140 (MQ=255)
                                                                          aTAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGttt                                                        <  2:508321/138‑1 (MQ=255)
                                                                          aTAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGttt                                                        >  1:508321/1‑138 (MQ=255)
                                                                           tAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGAc                                                     >  2:957709/1‑140 (MQ=255)
                                                                           tAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGAc                                                     <  1:147968/140‑1 (MQ=255)
                                                                                aTGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTg                                                <  2:675549/140‑1 (MQ=255)
                                                                                  gTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCg                                              <  1:74173/140‑1 (MQ=255)
                                                                                     aCCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAg                                           <  1:30673/140‑1 (MQ=255)
                                                                                     aCCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAg                                           <  1:924676/140‑1 (MQ=255)
                                                                                      ccGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAgg                                          >  2:804637/1‑140 (MQ=255)
                                                                                       cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                         <  2:11605/140‑1 (MQ=255)
                                                                                       cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                         >  1:480363/1‑140 (MQ=255)
                                                                                       cGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAggg                                         >  2:874219/1‑140 (MQ=255)
                                                                                         gAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGa                                       >  1:31187/1‑140 (MQ=255)
                                                                                         gAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGa                                       <  1:430920/140‑1 (MQ=255)
                                                                                          aGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGaa                                      <  2:897515/140‑1 (MQ=255)
                                                                                          aGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGaa                                      >  1:897515/1‑140 (MQ=255)
                                                                                               gAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTg                                 >  2:185974/1‑140 (MQ=255)
                                                                                               gAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTg                                 <  2:567805/140‑1 (MQ=255)
                                                                                                aacaacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGa                                >  1:601913/1‑140 (MQ=255)
                                                                                                 acaacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGaa                               <  2:11619/140‑1 (MQ=255)
                                                                                                 acaacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGaa                               <  1:985185/140‑1 (MQ=255)
                                                                                                  caacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAAc                              <  2:569001/140‑1 (MQ=255)
                                                                                                   aacTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAAcc                             <  1:777265/140‑1 (MQ=255)
                                                                                                    acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGaa                               <  1:329812/137‑1 (MQ=255)
                                                                                                    acTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGaa                               >  2:329812/1‑137 (MQ=255)
                                                                                                     cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACctc                           >  1:641465/1‑140 (MQ=255)
                                                                                                     cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACctc                           >  2:595794/1‑140 (MQ=255)
                                                                                                     cTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACctc                           >  1:700784/1‑140 (MQ=255)
                                                                                                          gTGACTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTcaca                      >  2:600261/1‑140 (MQ=255)
                                                                                                             aCTTCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGa                   <  1:282222/140‑1 (MQ=255)
                                                                                                               ttCTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACg                 >  2:151119/1‑140 (MQ=255)
                                                                                                                 cTGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGc               <  2:601913/140‑1 (MQ=255)
                                                                                                                  tGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa              <  1:600261/140‑1 (MQ=255)
                                                                                                                  tGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa              <  1:185974/140‑1 (MQ=255)
                                                                                                                  tGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa              <  1:151119/140‑1 (MQ=255)
                                                                                                                  tGCGCTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa              <  2:31187/140‑1 (MQ=255)
                                                                                                                   gcgcTAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCAt             <  2:480363/140‑1 (MQ=255)
                                                                                                                       tAAACGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCaaa         >  2:584337/1‑140 (MQ=255)
                                                                                                                         aaCGGGGACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAAtt       >  2:957520/1‑140 (MQ=255)
                                                                                                                              ggACGTTTGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATc  >  2:864241/1‑140 (MQ=255)
                                                                                                                                   tttGTGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCa              <  2:466956/123‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
AAAGCTTCCAGAGAGTAAAAGTGTTTTAAGAAGGAGCAATTCTATCGATTTTGATTTTGGGAAATCAACACGGCATAATTATGTCACCGGAGCCTGAACAACTCCGGTGACTTCTGCGCTAAACGGGGACGTTTATGCGCACATACAATCCAAACTCTCTTCTCCCTTCACAGATGCAGAAATGCACCTGCAATTCTTTGCATCTAGCGTTTGACCTCTGCGGAGGGGAAGCGTGAACCTCTCACAAGACGGCATCAAATTACATC  >  NC_000913/572234‑572499

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