Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,021 T→A 7.6% intergenic (+37/‑28) lomR → / → stfR pseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,0210TA 7.6% 256.7 / 4.6 109intergenic (+37/‑28)lomR/stfRpseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (54/45);  new base A (4/6);  total (58/51)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCC  >  NC_000913/1428875‑1429165
                                                                                                                                                  |                                                                                                                                                
tGCAGTTTAACCCGACCGAATCCGAGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGTCGACTGGCGCACTGACGGTTTCATCGTGTGTGTCGGGTATAAGTGCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTc                                                                                                                                                <  1:1066625/149‑1 (MQ=255)
 gCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTACa                                                                                                                                               >  1:593354/1‑149 (MQ=255)
    gTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGc                                                                                                                                            <  1:1399211/149‑1 (MQ=255)
     tttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCg                                                                                                                                           >  2:1404498/1‑149 (MQ=255)
     tttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCg                                                                                                                                           >  1:261702/1‑149 (MQ=255)
     tttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCg                                                                                                                                           >  2:403806/1‑149 (MQ=255)
      ttAACCCTACCGAATACGTGGACATTGATATTGCTAATCAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTGCTGATTAGCCATGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                                                                                                                                          <  2:1149578/149‑1 (MQ=255)
      ttAACCCGACCGAATCCGTGGCCATTGATATTGCTTCTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                                                                                                                                          <  2:133418/149‑1 (MQ=255)
      ttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGACGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                                                                                                                                          <  2:1157666/149‑1 (MQ=255)
      ttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                                                                                                                                          <  1:1437169/149‑1 (MQ=255)
       tAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCggg                                                                                                                                         <  2:1478496/149‑1 (MQ=255)
          cccGACCGAATCCGTGGCCATTGATATTGCTTATGACGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCTGAACCGGCGGGCtg                                                                                                                                      <  1:1512457/149‑2 (MQ=255)
          cccGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCtt                                                                                                                                      >  2:121166/1‑149 (MQ=255)
           ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCttt                                                                                                                                     >  2:443298/1‑149 (MQ=255)
           ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCttt                                                                                                                                     >  1:552157/1‑149 (MQ=255)
           ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCttt                                                                                                                                     >  2:1566516/1‑149 (MQ=255)
             gaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCttttt                                                                                                                                   >  1:18216/1‑149 (MQ=255)
              accgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCTGAACCGGCGGGCtgtcat                                                                                                                                  <  1:969494/149‑6 (MQ=255)
                  aaTCCGTGGCCATTGATATTGCTTATTAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTAATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCACGTAACACAGTGATATGACAGCCCGCCTGTTCAGGCCGGCTTTTTTGTgg                                                                                                                              >  2:1467873/1‑149 (MQ=255)
                  aaTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTgg                                                                                                                              >  1:479813/1‑149 (MQ=255)
                  aaTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTgg                                                                                                                              <  1:791988/149‑1 (MQ=255)
                   aTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                                             <  1:1405356/149‑1 (MQ=255)
                   aTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                                             >  2:123261/1‑149 (MQ=255)
                   aTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGATCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                                                                                                                             <  1:147411/149‑1 (MQ=255)
                         ggCCATTGCTATGGCTTATGCAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGCGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAat                                                                                                                       <  1:211631/149‑1 (MQ=255)
                             aTTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATgg                                                                                                                   >  1:546724/1‑149 (MQ=255)
                                   atTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTCATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGaa                                                                                                             >  2:732133/1‑149 (MQ=255)
                                      gCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCCGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTaaata                                                                                                          >  2:16948/1‑147 (MQ=255)
                                       catATGAAGGCCCCGTCAGTGGCGACTGGCGCACTGACGGTTTCCTCGTGGGTGGCGGTTATAAGTTCTGATTAGCCAGGTAACACATTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGAt                                                                                                         <  1:189138/147‑1 (MQ=255)
                                       cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGAt                                                                                                         >  1:820827/1‑149 (MQ=255)
                                            gAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAg                                                                                                    <  2:800266/149‑1 (MQ=255)
                                            gAAGGACCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATAGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAg                                                                                                    <  1:871358/149‑1 (MQ=255)
                                                  cccGGCAGTGGCGACTGGCGCACTGACGGTATCATCGTGGGTGTCGGTTATAAGTTCTGAGTAGCCAGGTAACACAGAGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTATTGTGGGGTGAATATGGAAGTAAAGATTTCAGGTTTAc                                                                                              >  2:536997/1‑149 (MQ=255)
                                                   ccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGTGTGAATATGGCAGTAAAGATTTCAGGTGTACt                                                                                             <  2:516191/149‑1 (MQ=255)
                                                    cGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTg                                                                                            >  2:815536/1‑149 (MQ=255)
                                                        aGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAg                                                                                        >  1:92555/1‑149 (MQ=255)
                                                         gTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGa                                                                                       >  1:1125326/1‑149 (MQ=255)
                                                          tGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTAAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGAc                                                                                      >  2:1288111/1‑149 (MQ=255)
                                                           ggCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACg                                                                                     <  1:121166/149‑1 (MQ=255)
                                                           ggCGACTGGCGCACTGACGGTTTCATCGGGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACg                                                                                     <  1:405564/149‑1 (MQ=255)
                                                              gACTCGCGCACTGACCGTTTCATGCTGGTTGTCGGTTATAAGTTCTGAAGAGCCAGGTAACACAGTGTTAAGACAGCCCGCCGGTTCAGCCGGGCTTTTTTGTGGAGTGAATATGGCAGTAAAGATTTCAGGAGTACAGAAAgcctgcc                                                                                  >  1:1058506/1‑143 (MQ=255)
                                                               aCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGcac                                                                                 >  1:205079/1‑149 (MQ=255)
                                                                  ggCGCACTGACGGTTCCCTCCTTGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCTGGCTTTTTTGTGGGGTGAAGATGGCAGTAAAGATTTCAGGTGTCCTGAAAGACGGCACAgg                                                                              <  1:1022780/149‑1 (MQ=255)
                                                                   gcgcACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGa                                                                             >  1:1204616/1‑149 (MQ=255)
                                                                     gcACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGaaa                                                                           <  1:758759/149‑1 (MQ=255)
                                                                     gcACTGACGGTTTCATCGTGGGTGTCGGTTATAAGATCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCACTAAACATTTCAGGTGTACTGAAACACGGCACAGTaaa                                                                           >  2:967792/1‑149 (MQ=255)
                                                                       aCTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAAc                                                                         <  1:84017/149‑1 (MQ=255)
                                                                         ttaCGGTTTCAACGAGCGTGTCGGCAATAAGTTCAGATTAGCCAGGTAAGACAGTGTTATGACAGCCAGCCGGTTCAGGCGGGCTTTTATGTTGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCg                                                                       <  2:925098/147‑1 (MQ=255)
                                                                         tGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAAATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCg                                                                       >  1:920550/1‑149 (MQ=255)
                                                                         tGACGGTTTCATCGTGGGGGTCGGTGATAAGTTCTGATAAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGTGTGAATATGGCAGTAAAGTTTTCAGGTGTACTGAAAGACGGCACAGGAAAACCg                                                                       <  1:1161118/149‑1 (MQ=255)
                                                                           aCGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGt                                                                     <  1:123261/149‑1 (MQ=255)
                                                                                atCATAGTGGGTGTCGGTTATAAGTTCTGATTAGACAGGTAACACAGTGTTATGACAGCCAGCCGGTTCACGCGGGCGTTATTGTGGGGTGAATATGGAAGTAAAGATTTCAGGTGTACTAAAAGACGGAACATAAAAACCGGTACAGa                                                                >  2:693337/2‑149 (MQ=255)
                                                                                   aTCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACt                                                             <  1:403806/149‑1 (MQ=255)
                                                                                    tCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTg                                                            >  1:1354773/1‑149 (MQ=255)
                                                                                     cGTGTGAGTCGGTTATAAGTGCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGc                                                           <  1:279052/149‑1 (MQ=255)
                                                                                      gTGGGTGTCGGTTATAAGTGCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGATCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGca                                                          <  2:1097999/149‑1 (MQ=255)
                                                                                      gTGGGGGTCGGGTATAAGGTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGca                                                          <  1:143756/149‑1 (MQ=255)
                                                                                       tGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGcac                                                         >  1:616696/1‑149 (MQ=255)
                                                                                        ggggGTCGGTTATAAGTGCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGATCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGcaca                                                        <  1:329524/149‑1 (MQ=255)
                                                                                         ggTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACaa                                                       >  2:699974/1‑149 (MQ=255)
                                                                                         ggTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACaa                                                       >  2:443137/1‑149 (MQ=255)
                                                                                           ggtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATc                                                     <  2:1271246/148‑1 (MQ=255)
                                                                                           ggtCCGTTATAAGTTCTGATTAGCCTGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATc                                                     <  2:98092/148‑1 (MQ=255)
                                                                                             tCGGTTATAATTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGTCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCa                                                   <  2:980794/149‑1 (MQ=255)
                                                                                              cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAg                                                  <  2:1440086/149‑1 (MQ=255)
                                                                                                ggtataAGTTCTGATTAGGCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACTGAACTGCACAATCCAGCt                                                <  2:970157/147‑1 (MQ=255)
                                                                                                gTTTTAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGTGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACGGGTACAGAACTGCACAATCCAGcg                                                >  1:411203/1‑148 (MQ=255)
                                                                                                gTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGCTTCAGCCGGGCTTTTTTGTCGGGAGAATATGGCAGTAAAGATTTCAGGTGAACTGAAAGACGTCACAGGAAAACCGGTACAGAAATGCACAATCCAACt                                                >  2:1434258/1‑149 (MQ=255)
                                                                                                   ataAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCTCAGGTTCAGGCGGGAATTTTAGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAAACCGTACACAACTGCACTATCCAGCTGaa                                             >  1:853809/1‑149 (MQ=255)
                                                                                                    taAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGaaa                                            <  2:1409193/149‑1 (MQ=255)
                                                                                                    taAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGaaa                                            >  1:162916/1‑149 (MQ=255)
                                                                                                       gTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCa                                         >  1:653612/1‑149 (MQ=255)
                                                                                                        ttCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCaa                                        <  1:1174889/149‑1 (MQ=255)
                                                                                                           tGAATAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGATCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAATACGGCACAGGAAAACCGGTACATAACTGCACAATCCAACTGAAAGCAAAAc                                     >  1:1484097/1‑149 (MQ=255)
                                                                                                             aTTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGt                                   >  1:945466/1‑149 (MQ=255)
                                                                                                             aTTAGACAGGTAACACAGTGTTATGACAGCCCTCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGGGTACTGAAAGACGGCACAGGAAACCCGGTCAAGAACTGCACAATACAGCTGCAAGAAAAACGt                                   >  2:1038287/1‑149 (MQ=255)
                                                                                                              ttAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTa                                  >  2:768221/1‑149 (MQ=255)
                                                                                                                aGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc                                <  2:261702/149‑1 (MQ=255)
                                                                                                                 gCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACa                               <  1:236586/149‑1 (MQ=255)
                                                                                                                  ccAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAg                              <  1:443298/149‑1 (MQ=255)
                                                                                                                  ccAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGGGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAg                              >  1:1434494/1‑149 (MQ=255)
                                                                                                                        aaCACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGcacccc                        >  1:246535/1‑147 (MQ=255)
                                                                                                                        aaCACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGcaccac                        >  1:1219737/1‑149 (MQ=255)
                                                                                                                            caGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtg                    <  2:820827/149‑1 (MQ=255)
                                                                                                                             aGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtgg                   >  2:560388/1‑149 (MQ=255)
                                                                                                                               tgtTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtg                 <  2:459625/149‑1 (MQ=255)
                                                                                                                                 ttATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt               <  1:1363761/149‑1 (MQ=255)
                                                                                                                                 ttATGACAGCCCGCCGGTGCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt               <  1:1300047/149‑1 (MQ=255)
                                                                                                                                  tATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGgtggtggtg              >  2:909805/1‑149 (MQ=255)
                                                                                                                                   aTGCCAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATTTGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCGCAGGAAAACCCCTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGa             <  2:1058506/149‑1 (MQ=255)
                                                                                                                                    tGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGaa            <  2:701018/149‑1 (MQ=255)
                                                                                                                                     gACAGCCCTGAGGATCCGGCGGGCTTTTTTGTGGGGTGTTTATGGCAGTAAAGATTTCAGGTGTACTGAAAGTCGGCACAGGAACACCGGGAGAGACCAGCTCAATACAGCTGATAGCAAAACTTAACAGCACGACGGGGGTGGTGCac           >  2:1194878/1‑149 (MQ=255)
                                                                                                                                      aaaGCCCGCCTGAACCGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAaca          >  2:955438/3‑149 (MQ=255)
                                                                                                                                       cAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCCCGGTGGTGGTGAaccc         >  1:1309226/1‑147 (MQ=255)
                                                                                                                                        aGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACg        >  2:1110540/1‑149 (MQ=255)
                                                                                                                                           ccggcgGATCAGGCGGTCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACTGCACAGGAAAACCGGGACAGAACTGAACAATCCATCTTACAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTg     <  1:1467873/149‑1 (MQ=255)
                                                                                                                                            cgccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTgg    >  1:1412666/1‑149 (MQ=255)
                                                                                                                                             gccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGc   <  1:454603/149‑1 (MQ=255)
                                                                                                                                             gccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTAAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGc   >  2:1181818/1‑149 (MQ=255)
                                                                                                                                              ccgGTTCATGCGGGCATAATTGTGGGGAGAATATTGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCAAACGAAAACCGGTACAGAACTGCACCATCCAGCTGCACGCAAAACGTAACAGCACCACGGTGCAGGAGAACACGCTGGcc  >  2:1545168/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCC  >  NC_000913/1428875‑1429165

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: