Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 269,158 C→T intergenic (‑198/+131) ykfA ← / ← perR CP4‑6 prophage; putative GTP‑binding protein YkfA/putative transcriptional regulator PerR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913269,1580CT80.9% 571.7 / 124.4 226intergenic (‑198/+131)ykfA/perRCP4‑6 prophage; putative GTP‑binding protein YkfA/putative transcriptional regulator PerR
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base T (88/94);  total (107/119)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-02

Only 100 of 241 total aligned reads displayed.

CAATCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATAAAGTTACT  >  NC_000913/269059‑269255
                                                                                                   |                                                                                                 
cAATCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAgcc                                                                                                  >  2:357788/1‑101 (MQ=255)
 aaTCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAgtcg                                                                                                 <  2:357174/101‑1 (MQ=255)
  aTCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAgtcgc                                                                                                >  2:357792/1‑101 (MQ=255)
  aTCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAgtcgc                                                                                                >  1:357790/1‑101 (MQ=255)
   tCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAgtcgcc                                                                                               >  2:357794/1‑101 (MQ=255)
    ccTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGtcgccg                                                                                              <  1:357324/101‑1 (MQ=255)
     ctctGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGt                                                                                             <  1:357216/101‑1 (MQ=255)
      tctGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTg                                                                                            <  1:357402/101‑1 (MQ=255)
      tctGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTg                                                                                            <  2:357302/101‑1 (MQ=255)
       ctGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGa                                                                                           >  2:357797/1‑101 (MQ=255)
       ctGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGa                                                                                           <  1:357312/101‑1 (MQ=255)
        tGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGAt                                                                                          <  1:356980/101‑1 (MQ=255)
        tGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGAt                                                                                          >  1:357802/1‑101 (MQ=255)
        tGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGAt                                                                                          >  1:357800/1‑101 (MQ=255)
        tGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGAt                                                                                          <  2:357346/101‑1 (MQ=255)
           gcTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATaca                                                                                       >  2:357808/1‑101 (MQ=255)
           gcTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATaca                                                                                       >  2:357806/1‑101 (MQ=255)
             ttGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACAcc                                                                                     >  2:357811/1‑101 (MQ=255)
               gAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCtt                                                                                   <  2:357259/101‑1 (MQ=255)
                aGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTa                                                                                  >  2:357813/1‑101 (MQ=255)
                 ggCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTaa                                                                                 <  2:357516/101‑1 (MQ=255)
                   cAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTaaaa                                                                               >  2:357815/1‑101 (MQ=255)
                     gCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAAtc                                                                             <  2:357593/101‑1 (MQ=255)
                       tGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAAtctc                                                                           >  2:357818/1‑101 (MQ=255)
                            gACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAAccc                                                                      <  2:357439/101‑1 (MQ=255)
                             aCTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCa                                                                     <  1:357523/101‑1 (MQ=255)
                             aCTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCa                                                                     <  1:357415/101‑1 (MQ=255)
                              ctgctgGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAg                                                                    >  1:357821/1‑101 (MQ=255)
                               tgctgGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGc                                                                   >  1:357823/1‑101 (MQ=255)
                               tgctgGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGc                                                                   <  2:357358/101‑1 (MQ=255)
                                 ctgGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCaa                                                                 <  1:357354/101‑1 (MQ=255)
                                 ctgGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCaa                                                                 <  1:357281/101‑1 (MQ=255)
                                  tgGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCaaa                                                                >  2:357825/1‑101 (MQ=255)
                                   ggCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAg                                                               <  2:357623/101‑1 (MQ=255)
                                   ggCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAg                                                               >  1:357826/1‑101 (MQ=255)
                                   ggCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAg                                                               <  1:357500/101‑1 (MQ=255)
                                    gCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGt                                                              <  2:357383/101‑1 (MQ=255)
                                     cTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGtt                                                             <  1:357594/101‑1 (MQ=255)
                                      tCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGttt                                                            >  1:357829/1‑101 (MQ=255)
                                       cAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTc                                                           <  2:357551/101‑1 (MQ=255)
                                       cAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTc                                                           >  1:357831/1‑101 (MQ=255)
                                          gCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGa                                                        <  2:357538/101‑1 (MQ=255)
                                           cAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAg                                                       <  2:357564/101‑1 (MQ=255)
                                           cAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAg                                                       <  2:357207/101‑1 (MQ=255)
                                            aatgaatgAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGc                                                      >  1:357836/1‑101 (MQ=255)
                                             atgaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGcc                                                     <  2:357419/101‑1 (MQ=255)
                                             atgaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGcc                                                     <  1:357472/101‑1 (MQ=255)
                                             atgaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGcc                                                     <  2:357409/101‑1 (MQ=255)
                                              tgaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCcg                                                    >  2:357837/1‑101 (MQ=255)
                                              tgaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCcg                                                    >  2:357839/1‑101 (MQ=255)
                                              tgaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCcg                                                    >  1:357841/1‑101 (MQ=255)
                                               gaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCcgc                                                   <  2:357457/101‑1 (MQ=255)
                                               gaatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCcgc                                                   <  2:357690/101‑1 (MQ=255)
                                                aatgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCcgcg                                                  >  1:357843/1‑101 (MQ=255)
                                                 atgaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGa                                                 <  2:357386/101‑1 (MQ=255)
                                                   gaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGc                                               <  1:357495/101‑1 (MQ=255)
                                                   gaGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGc                                               <  1:357520/101‑1 (MQ=255)
                                                     gTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCaa                                             >  2:357846/1‑101 (MQ=255)
                                                      ttATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCaaa                                            <  2:357504/101‑1 (MQ=255)
                                                       tataATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCaaaa                                           <  2:357599/101‑1 (MQ=255)
                                                        ataataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAg                                          >  1:357852/1‑101 (MQ=255)
                                                        ataataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAg                                          >  1:357850/1‑101 (MQ=255)
                                                         taataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGt                                         >  1:357855/1‑101 (MQ=255)
                                                         taataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAACCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGt                                         >  2:357856/1‑101 (MQ=255)
                                                          aataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTg                                        >  1:357859/1‑101 (MQ=255)
                                                           ataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTga                                       >  1:357861/1‑101 (MQ=255)
                                                           ataGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTga                                       <  2:357352/101‑1 (MQ=255)
                                                              gcagcaTTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGagag                                    >  2:357862/1‑101 (MQ=255)
                                                                agcaTTTACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGcc                                  >  1:357865/1‑101 (MQ=255)
                                                                   aTTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTc                               >  2:357869/1‑101 (MQ=255)
                                                                    tttACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTca                              >  1:357871/1‑101 (MQ=255)
                                                                    tttACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTca                              >  2:357872/1‑101 (MQ=255)
                                                                      tACTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTcaca                            <  2:357474/101‑1 (MQ=255)
                                                                       aCTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACaa                           <  1:357347/101‑1 (MQ=255)
                                                                       aCTAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACaa                           <  2:357524/101‑1 (MQ=255)
                                                                         tAACAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACaaaa                         <  2:357800/101‑1 (MQ=255)
                                                                          aaCAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaat                        >  1:357876/1‑101 (MQ=255)
                                                                          aaCAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaat                        <  1:357371/101‑1 (MQ=255)
                                                                          aaCAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaat                        <  1:357513/101‑1 (MQ=255)
                                                                           aCAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaata                       <  1:357396/101‑1 (MQ=255)
                                                                            cAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaataa                      >  2:357878/1‑101 (MQ=255)
                                                                            cAGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaataa                      <  1:357453/101‑1 (MQ=255)
                                                                             aGGGATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaataat                     <  1:357556/101‑1 (MQ=255)
                                                                             aGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaataat                     >  1:357880/1‑101 (MQ=255)
                                                                              gggATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATg                    <  2:357481/101‑1 (MQ=255)
                                                                              gggATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATg                    >  2:357883/1‑101 (MQ=255)
                                                                                gATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCt                  >  2:357887/1‑101 (MQ=255)
                                                                                gATTTATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCt                  >  1:357886/1‑101 (MQ=255)
                                                                                  tttATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTa                <  2:357158/101‑1 (MQ=255)
                                                                                    tATTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGt              >  1:357891/1‑101 (MQ=255)
                                                                                     aTTGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGta             <  1:357581/101‑1 (MQ=255)
                                                                                      ttGAGAGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCGAAAGTGAGAGCCTTCACAAAATAATGCTTAGtaa            >  1:357893/1‑101 (MQ=255)
                                                                                        gagaGTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGtaata          >  1:357896/1‑101 (MQ=255)
                                                                                         agagTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGtaataa         >  2:357897/1‑101 (MQ=255)
                                                                                          gagTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATaaa        >  2:357900/1‑101 (MQ=255)
                                                                                           agTATGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATAAAg       >  1:357901/1‑101 (MQ=255)
                                                                                            gTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATAAAGt      <  2:357671/101‑1 (MQ=255)
                                                                                            gTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATAAAGt      <  1:357462/101‑1 (MQ=255)
                                                                                              aTGAGTCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATAAAGTTa    <  1:357483/101‑1 (MQ=255)
                                                                                                   ccgccgTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAaataat                     <  1:357715/79‑1 (MQ=255)
                                                                                                   |                                                                                                 
CAATCCTCTGCGCTTGAGGCAGCTGCGCGACTGCTGGCTCAGGCAATGAATGAGTTATAATAGCAGCATTTACTAACAGGGATTTATTGAGAGTATGAGCCGCCGTGATACACCTTAAAATCTCAACCCAGCAAAGTTTCGGAGCCGCGAGCAAAAGTGAGAGCCTTCACAAAATAATGCTTAGTAATAAAGTTACT  >  NC_000913/269059‑269255

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: