Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 273,534 Δ1 bp 5.0% pseudogene (688/1108 nt) ykfC → CP4‑6 prophage; conserved protein;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913273,5340A. 5.0% 336.0 / 9.2 99pseudogene (688/1108 nt)ykfCCP4‑6 prophage; conserved protein;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (49/45);  new base . (2/3);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 101 total aligned reads displayed.

GAAGGCTGTACGCCGCAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGT  >  NC_000913/273386‑273681
                                                                                                                                                    |                                                                                                                                                    
gaAGGCTGTACGCCGCAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAAca                                                                                                                                                      >  1:901088/1‑149 (MQ=255)
      tGTACGCCGCAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGc                                                                                                                                                <  2:360482/149‑1 (MQ=255)
       gTACGCCGCAGGAGCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCt                                                                                                                                               <  2:1443268/149‑1 (MQ=255)
          cgccgcAGGATCAGGGAGGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGCTGTCGGGCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGaa                                                                                                                                            <  1:1042916/149‑1 (MQ=255)
           gccgcAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAAt                                                                                                                                           >  2:1024748/1‑149 (MQ=255)
            ccgcAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATg                                                                                                                                          >  1:702094/1‑149 (MQ=255)
             cgcAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGa                                                                                                                                         <  1:917572/149‑1 (MQ=255)
             cgcAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGa                                                                                                                                         >  1:1416164/1‑149 (MQ=255)
                aGGATCAGTGACGCTCGTTTCATGACTCGGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGGGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGtt                                                                                                                                     <  2:1422906/149‑1 (MQ=255)
                 ggATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTc                                                                                                                                    <  1:741478/149‑1 (MQ=255)
                   agcAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGa                                                                                                                                  <  2:523830/147‑1 (MQ=255)
                     cAGTGACGCACGTTTCATGACTCTGCGGTGGAAAACCCTCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATc                                                                                                                                <  2:459641/149‑1 (MQ=255)
                        tGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAAt                                                                                                                             <  1:1365953/149‑1 (MQ=255)
                        tGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACTCATTCAGCATGATGTTCGAtaata                                                                                                                             <  2:598055/149‑5 (MQ=255)
                         gACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACTCATTCAGCATGATGTTCGAtaatag                                                                                                                            <  1:485320/149‑6 (MQ=255)
                         gACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACTCATTCAGCATGATGTTCGAtaatag                                                                                                                            <  1:1464992/149‑6 (MQ=255)
                              aCGTTTCATGACTCTGCTGTGGAAAACCATCAACGCGGGACATATCGATGTAGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGc                                                                                                                       >  1:1340350/1‑149 (MQ=255)
                              aCGTTTCATGACTCTGCTGTGGAAAACAATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCCGTGTTATATCGCCGCTATTATCGAACATCATTCTGAATGAGTTCGATCAATACCTGc                                                                                                                       >  2:676612/1‑149 (MQ=255)
                                        aCTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTa                                                                                                             >  2:184982/1‑149 (MQ=255)
                                         ctctGCTGTGGTAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTAc                                                                                                            <  2:326287/149‑1 (MQ=255)
                                         ctctGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTAc                                                                                                            <  2:588016/149‑1 (MQ=255)
                                          tctGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTAcc                                                                                                           >  1:579217/1‑149 (MQ=255)
                                          tctGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTAcc                                                                                                           >  2:1374570/1‑149 (MQ=255)
                                            tgctgTGGAAAACCATCAAGGCGGGTCATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGAt                                                                                                                                 <  2:1187396/125‑1 (MQ=255)
                                            tgctgTGGAAAACCATCAAGGCGGGTCATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGAt                                                                                                                                 >  1:1187396/1‑125 (MQ=255)
                                                gtgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTAATATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCg                                                                                                     >  1:1475621/1‑149 (MQ=255)
                                                gaggAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGACCg                                                                                                     >  1:939969/3‑149 (MQ=255)
                                                 tgGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCgg                                                                                                    <  2:1537704/149‑1 (MQ=255)
                                                   gAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGa                                                                                                  >  1:1321737/1‑149 (MQ=255)
                                                   gAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGa                                                                                                  >  1:462242/1‑149 (MQ=255)
                                                     aaaCCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATAGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGaaa                                                                                                <  1:186718/149‑1 (MQ=255)
                                                      aaCCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAATGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGATACCTGAGCGGGAAAg                                                                                               >  2:33463/1‑149 (MQ=255)
                                                       aCCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGc                                                                                              <  1:1255231/149‑1 (MQ=255)
                                                        ccATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGcc                                                                                             >  2:1403242/1‑149 (MQ=255)
                                                          aTCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAg                                                                                           >  1:130248/1‑149 (MQ=255)
                                                           tCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGa                                                                                          >  2:660006/1‑149 (MQ=255)
                                                           tCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGa                                                                                          >  2:870060/1‑149 (MQ=255)
                                                                  gggACATATCGATGTGGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATCTCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATc                                                                                   <  1:307223/149‑1 (MQ=255)
                                                                  gggACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATc                                                                                   <  2:699617/149‑1 (MQ=255)
                                                                   ggACATATCGATGTCGGTCTCTGTCGGGCGGCCAGTGAAGGTGTGCCACATGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACATTCATGATCGCTACATTATCGGGAAAGCCAGaaaaaaacc                                                                                  >  2:783810/1‑144 (MQ=255)
                                                                    gACATATCGATGTCGGGCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGGTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCgg                                                                                 <  1:318822/149‑1 (MQ=255)
                                                                     aCATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggt                                                                                <  2:1378098/149‑1 (MQ=255)
                                                                     aCATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggt                                                                                <  1:545196/149‑1 (MQ=255)
                                                                      cATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggtg                                                                               <  2:646137/149‑1 (MQ=255)
                                                                       atatCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCggtgg                                                                              >  2:1107951/1‑149 (MQ=255)
                                                                            gATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGAGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTg                                                                         >  2:734400/1‑149 (MQ=255)
                                                                              tGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGa                                                                       <  2:893805/149‑1 (MQ=255)
                                                                               gTCGGTCACTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGaa                                                                      <  2:275602/149‑1 (MQ=255)
                                                                                  ggTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATaa                                                                   <  2:584933/149‑1 (MQ=255)
                                                                                  ggTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATaa                                                                   <  1:1137581/149‑1 (MQ=255)
                                                                                  ggTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATaa                                                                   <  2:1540540/149‑1 (MQ=255)
                                                                                    tctcTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCTTGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACa                                                                 <  1:985816/149‑1 (MQ=255)
                                                                                    tctcTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCAAGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACa                                                                 <  1:189149/149‑1 (MQ=255)
                                                                                     ctctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAg                                                                >  1:243701/1‑149 (MQ=255)
                                                                                      tctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGt                                                               <  2:1060309/149‑1 (MQ=255)
                                                                                      tctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGt                                                               >  1:491630/1‑149 (MQ=255)
                                                                                      tctTTCGGGCGGCCAGTGAAGGTGAGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGAGGTACAGGAATAACAGt                                                               >  1:1145585/1‑149 (MQ=255)
                                                                                       ctTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTa                                                              <  2:1340350/149‑1 (MQ=255)
                                                                                           cGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCa                                                          >  2:1058262/1‑149 (MQ=255)
                                                                                                 gcCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACgggg                                                    >  2:20624/1‑149 (MQ=255)
                                                                                                 gcCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTCCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTCGTACTGGAATAACAGTATCGACCgggg                                                    >  1:819477/1‑149 (MQ=255)
                                                                                                   cAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGcc                                                  >  1:334700/1‑149 (MQ=255)
                                                                                                     gTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGTTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGa                                                >  2:1288110/1‑149 (MQ=255)
                                                                                                      tGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGATTGAGTTCGATCAATAGCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGaa                                               <  2:134362/149‑1 (MQ=255)
                                                                                                      tGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGaa                                               <  2:396491/149‑1 (MQ=255)
                                                                                                      tGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGaa                                               >  1:1160823/1‑149 (MQ=255)
                                                                                                         aGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTa                                            >  1:400757/1‑149 (MQ=255)
                                                                                                           gtgtGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTAcg                                          >  1:1266196/1‑149 (MQ=255)
                                                                                                            ggtgCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTAcgg                                         <  2:403513/148‑1 (MQ=255)
                                                                                                                ccACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGt                                     <  2:702094/149‑1 (MQ=255)
                                                                                                                 cacaGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTc                                    >  2:911206/1‑149 (MQ=255)
                                                                                                                 cacaGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGAGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTc                                    >  2:542175/1‑149 (MQ=255)
                                                                                                                   caGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCag                                  >  1:670730/1‑149 (MQ=255)
                                                                                                                    aGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCaga                                 <  1:1024748/149‑1 (MQ=255)
                                                                                                                    aGGGCGGGGTTATATCGCCGCTATTATCGAAAATCAGGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCaga                                 <  1:766081/149‑1 (MQ=255)
                                                                                                                      ggCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAgaga                               >  2:45592/1‑149 (MQ=255)
                                                                                                                             ttATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTgg                        <  1:184982/149‑1 (MQ=255)
                                                                                                                             ttATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTgg                        >  2:1476587/1‑149 (MQ=255)
                                                                                                                             ttATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTgg                        <  2:870474/149‑1 (MQ=255)
                                                                                                                             gtataGCGCCGCTATTATATAAAATCATGCTGAATGAGTTCGATCAATACCTTCCTGACCGCTACCTGACCGGGACAGCCAGAAAAGATCGGTGGTACTGTAATAACAGTATCCAACTGTGCCTAAGTACGGCGGTCATAGAAAACTgg                        <  2:237733/148‑1 (MQ=255)
                                                                                                                               atatCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCa                      >  2:724044/1‑149 (MQ=255)
                                                                                                                               atatCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCa                      >  2:111166/1‑149 (MQ=255)
                                                                                                                               atatCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCa                      >  2:1565962/1‑149 (MQ=255)
                                                                                                                                tatCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAg                     >  1:1015617/1‑149 (MQ=255)
                                                                                                                                  tCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTg                   >  1:1366077/1‑149 (MQ=255)
                                                                                                                                  tCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACCGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTg                   <  2:1416164/149‑1 (MQ=255)
                                                                                                                                   cgccgcTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTgg                  <  2:358451/149‑1 (MQ=255)
                                                                                                                                      cgcTATTATCGAACATCATGCTGAATGAGTGCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGaaa               <  2:130248/149‑1 (MQ=255)
                                                                                                                                        tattgatCGAACTCATTCAGCATGATGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAAcc             >  2:1414540/6‑149 (MQ=255)
                                                                                                                                         tattatCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAAccc            >  1:1200152/1‑149 (MQ=255)
                                                                                                                                         tattatCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAAccc            >  1:1206356/1‑149 (MQ=255)
                                                                                                                                         tattatCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAAccc            >  1:218114/1‑149 (MQ=255)
                                                                                                                                          attatCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCcg           <  1:724044/149‑1 (MQ=255)
                                                                                                                                             atCGAACTCATTCAGCATGATGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCgg        >  2:750059/1‑149 (MQ=255)
                                                                                                                                              tCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGt       >  2:396659/1‑149 (MQ=255)
                                                                                                                                              tCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGt       >  2:1557190/1‑149 (MQ=255)
                                                                                                                                                 aaCATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTg                   <  2:663566/134‑1 (MQ=255)
                                                                                                                                                 aaCATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTg                   >  1:663566/1‑134 (MQ=255)
                                                                                                                                                  aCATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCg   <  1:3190/149‑1 (MQ=255)
                                                                                                                                                   catcatGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCCGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCCAAGTACGGATGTCAGAGAAAACTGGCAGTGGAAACCCGCGGAGGCGt  >  2:173878/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGGCTGTACGCCGCAGGATCAGTGACGCACGTTTCATGACTCTGCTGTGGAAAACCATCAAGGCGGGACATATCGATGTCGGTCTCTTTCGGGCGGCCAGTGAAGGTGTGCCACAGGGCGGTGTTATATCGCCGCTATTATCGAACATCATGCTGAATGAGTTCGATCAATACCTGCATGAGCGCTACCTGAGCGGGAAAGCCAGAAAAGATCGGTGGTACTGGAATAACAGTATCCAACGGGGCCGAAGTACGGCGGTCAGAGAAAACTGGCAGTGGAAACCCGCGGTGGCGT  >  NC_000913/273386‑273681

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