Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 289,256 G→C 7.7% intergenic (‑94/+45) yagI ← / ← argF CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913289,2560GC 7.7% 228.8 / 9.1 72intergenic (‑94/+45)yagI/argFCP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base C (2/3);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

GTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCC  >  NC_000913/289109‑289400
                                                                                                                                                     |                                                                                                                                                
gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgt                                                                                                                                                 >  2:163479/1‑149 (MQ=255)
gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgt                                                                                                                                                 <  8:206799/149‑1 (MQ=255)
 ttGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgtg                                                                                                                                                <  7:281562/149‑1 (MQ=255)
 ttGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGAATAATGCTGGCATGTCCACGCTGGAATTAAGGCTAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgtg                                                                                                                                                <  7:139973/149‑1 (MQ=255)
   gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATTTCCACGCTGGAATTAAGGCGATAACGTACCTTTTTCGGGATGAAAAATCATTGAAATGGAGCCGGT‑GT‑TCTCCCTCTGCCTGTgg                                                                                                                                               <  2:265137/148‑1 (MQ=255)
   gAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTgg                                                                                                                                               <  5:227497/148‑1 (MQ=255)
      cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGgaga                                                                                                                                           <  1:163479/149‑1 (MQ=255)
       aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAATGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGagag                                                                                                                                          <  6:295571/149‑1 (MQ=255)
        ggTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCACAGGGAGAgg                                                                                                                                         <  6:375387/149‑1 (MQ=255)
        ggTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCGCCCTCTCCCACAGGGAGAgg                                                                                                                                         <  3:352959/149‑1 (MQ=255)
             aGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCgg                                                                                                                                    <  2:55485/149‑1 (MQ=255)
                  cTGCAACGCACGGTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                               <  1:160769/149‑1 (MQ=255)
                   tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                 >  2:160769/1‑146 (MQ=255)
                     cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAggg                                                                                                                            <  6:36977/149‑1 (MQ=255)
                       aCGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                          >  7:348327/1‑149 (MQ=255)
                        ctcACGTTCAAAAGACTTATTAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACTTAGCTTTTTCGGGATGAAAATTCAATGAATTGGCGCCGGT‑TT‑TCACCCTCTCCGTGTGGGAGAGGGCGGGGGTGAGGGCAc                                                                                                                         <  6:17425/147‑1 (MQ=255)
                         gCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAc                                                                                                                         >  1:15151/1‑148 (MQ=255)
                           aCGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCaa                                                                                                                      >  3:68006/1‑149 (MQ=255)
                             gTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcg                                                                                                                    >  1:233658/1‑149 (MQ=255)
                             gTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcg                                                                                                                    >  3:295926/1‑149 (MQ=255)
                                 aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgtg                                                                                                                                                <  6:224610/117‑1 (MQ=255)
                                 aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgtg                                                                                                                                                >  5:224610/1‑117 (MQ=255)
                                 aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgcg                                                                                                                  >  3:99530/1‑147 (MQ=255)
                                 aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgcg                                                                                                                  >  1:274595/1‑147 (MQ=255)
                                 aaCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAcgcg                                                                                                                  >  5:203519/1‑147 (MQ=255)
                                  aCAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAg                                                                                                               >  4:92309/1‑149 (MQ=255)
                                   cAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGc                                                                                                              >  7:194041/1‑149 (MQ=255)
                                        tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc                                                                                                         >  6:30741/1‑149 (MQ=255)
                                        tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc                                                                                                         >  1:170374/1‑149 (MQ=255)
                                         gAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCCCAACGCGCAGCGGAcc                                                                                                         >  2:174171/1‑148 (MQ=255)
                                          aatatgcGGCATAATGGTGGCATGTCCCCGATTGAATTACGGCGAAAACGTACCTTTTTCGTGATTAAAAATCATTCAAATGGACCCGGT‑GT‑TCTCCCTCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAccc                                                                                                        <  6:350807/142‑1 (MQ=255)
                                            taatCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGgaga                                                                                                                                           <  8:218967/111‑1 (MQ=255)
                                            taatCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGgaga                                                                                                                                           >  7:218967/1‑111 (MQ=255)
                                                  gCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcc                                                                                               >  3:195441/1‑149 (MQ=255)
                                                    aTAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgt                                                                                                                                                 <  3:230120/97‑1 (MQ=255)
                                                    aTAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCtgt                                                                                                                                                 >  4:230120/1‑97 (MQ=255)
                                                            ggCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                                                                                     >  2:221239/1‑149 (MQ=255)
                                                               aTGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCAAGCGTTGCca                                                                                  >  6:123846/1‑149 (MQ=255)
                                                                 gcccACGCTGTAATTAAGGCTAAAACTTACTTTTTTGGGGTTCAAAAATCATTGAATTGGCGCCGGT‑GT‑TCTCCCGCCCCCTGGGGGTGTGGGCCGGGGGGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatc                                                                                <  4:295926/147‑1 (MQ=255)
                                                                    cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcc                                                                                               >  1:212278/1‑131 (MQ=255)
                                                                    cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatcat                                                                              >  7:374730/1‑148 (MQ=255)
                                                                    aaCGCTGTAATTAAGGCTAAAACGTACCTTTTTCGGGATTAAAAATCAATGAAATGGAGCCGGT‑TT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcc                                                                                               <  2:212278/130‑1 (MQ=255)
                                                                       ccttgAATTACGGCCAAAACGTACCTTTTGCGTGATCAAAAAGAAATCAAATGGAGCGGGT‑GG‑TCTCCCTCCCCCTGGGGTCGTGGGCGGGGGTGAGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc                                                                          <  1:157058/145‑1 (MQ=255)
                                                                                 aGGCGAAACCGTACCTTTTTCGGGATTAAAAATCAATGAAAGGGCGCCGGT‑GT‑GCTCCCTCTCCCTGGGGGAGAGGGCCGGGGGGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATc                                                                 <  3:92309/148‑1 (MQ=255)
                                                                                 aGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCg                                                                >  4:321418/1‑149 (MQ=255)
                                                                                   gctaaaaCGTACCTTTTTCGGGATTAAAAATCAATTAATTGGCGCGGGT‑GT‑TCTCCCTCGCCCGGGGGTTGTGGGCCGGGGGGAGGGCACCAACGCGCAGGGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                              <  4:99530/146‑1 (MQ=255)
                                                                                   gCGAAAACGTACCTTTTTCGGGATTAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                              <  2:274595/149‑1 (MQ=255)
                                                                                    ctaaaaCGTACCTTTTTCGGGATGAAAAATCAATCAAAGGGACCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                             <  2:233658/147‑1 (MQ=255)
                                                                                    cGAAACCGTACCTTTTTCGGTATTAAAAATCAATTAAATGGAGCCGGG‑GT‑TCTCCCTCTCCCTGGGGGAGCGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                             <  2:280394/149‑1 (MQ=255)
                                                                                    cGAAACCGACCCTTTTTGGGGATGAAAAATCAATGAAAGGGACCCGGT‑GT‑GCTCCCTCTCCCTGGGGGAGAGGGCCGGGGGGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                             <  5:157200/149‑1 (MQ=255)
                                                                                    acaaacCGTACTTTTTGGGGTATCAAAAAACAAT‑AAATGGGCGCGGG‑GTGTCTCCCTCTCCGTGGGGGCGTGGGCGGGGGGGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAt                                                             <  5:30741/147‑1 (MQ=255)
                                                                                     gAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATg                                                            >  3:100977/1‑149 (MQ=255)
                                                                                         aCGTCCCTTTTGCGGGACTAAAAATCATTGAACGGCCCCGGGT‑GT‑GCTCCCTCGCCGTGGGGGAGTGGGCCGGGGGGTGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATc                                                        <  5:123846/149‑1 (MQ=255)
                                                                                           gTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAAt                                                                                                     >  3:208066/1‑102 (MQ=255)
                                                                                           gTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCg                                                      >  7:338445/1‑149 (MQ=255)
                                                                                             tcatttttCGGTGTTAAAAATACATTAAAGGGTGGCGGTGGT‑GCTCCCGCTCCGTGGGTGTTTGGGCGGGGGTGAGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCg                                                      <  8:348327/145‑1 (MQ=255)
                                                                                              ccTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGtt                                                   <  1:221239/149‑1 (MQ=255)
                                                                                               cTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCACAGGGAGAgg                                                                                                                                         >  5:262323/1‑62 (MQ=38)
                                                                                               cTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCACAGGGAGAgg                                                                                                                                         <  6:262323/62‑1 (MQ=38)
                                                                                                       ggATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCACCGGACCCAATCACTCCCCAACCGTTCCCATCATCACCGCCTTAATCTTATCCAACCGGTTTACCCCCTg                                          >  2:248362/1‑149 (MQ=38)
                                                                                                        gATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCACGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTgg                                         >  2:344650/1‑149 (MQ=255)
                                                                                                           gAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGACATCATCACCGCCTTAATCGTTTGCATCCGGGTTTCCGCCTGGTc                                       >  8:379154/1‑148 (MQ=255)
                                                                                                                    aTGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATg                             >  6:367559/1‑149 (MQ=255)
                                                                                                                       aaaTGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCa                                                                                                                >  7:126960/1‑63 (MQ=255)
                                                                                                                        aaTGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAac                                   >  7:205414/1‑139 (MQ=255)
                                                                                                                                 ggT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACt                >  7:265137/1‑149 (MQ=255)
                                                                                                                                   t‑gt‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGGTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGCTGCTCGCCGCCGACTc               >  3:201894/1‑148 (MQ=255)
                                                                                                                                         ctccctctccctGTGGGAGTGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAac           <  1:248362/148‑1 (MQ=255)
                                                                                                                                          accctctccgtGTGGTTGTGGGCGGGGGGGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAaca          <  1:174171/147‑1 (MQ=255)
                                                                                                                                            cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGTGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt       >  2:80017/1‑149 (MQ=255)
                                                                                                                                            cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGCATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt       >  1:121457/1‑149 (MQ=255)
                                                                                                                                            cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACAAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATACGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt       >  3:260085/1‑149 (MQ=255)
                                                                                                                                            cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt       >  7:149312/1‑149 (MQ=255)
                                                                                                                                            cctctcccaAAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCt       >  2:368124/1‑149 (MQ=255)
                                                                                                                                              tctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTgg                                         >  7:203392/1‑113 (MQ=255)
                                                                                                                                              tctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTgg                                         <  8:203392/113‑1 (MQ=255)
                                                                                                                                                 cgctGGGGGTGAGGGCGGGGGTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtcc  <  2:170374/147‑1 (MQ=255)
                                                                                                                                                  cctGTGGTTGTGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtc   <  2:15151/147‑1 (MQ=255)
                                                                                                                                                  cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                       <  8:155765/127‑1 (MQ=255)
                                                                                                                                                  cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                       >  7:155765/1‑127 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGT‑GT‑TCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCC  >  NC_000913/289109‑289400

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