Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,078,306:1 +G 6.3% intergenic (+199/+243) yeeW → / ← yoeD CP4‑44 prophage; uncharacterized protein;Phage or Prophage Related/pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,078,3061.G 6.3% 451.7 / 23.5 111intergenic (+199/+243)yeeW/yoeDCP4‑44 prophage; uncharacterized protein;Phage or Prophage Related/pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base . (58/46);  new base G (5/2);  total (63/48)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 114 total aligned reads displayed.

TAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACA  >  NC_000913/2078157‑2078451
                                                                                                                                                      |                                                                                                                                                 
aaGTGCCATAAACCACAGAAACAGCACTGTTGTCCTGTATAAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCTCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACa                                                                                                                                                    <  1:216343/149‑1 (MQ=255)
 aGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑g                                                                                                                                                  <  1:1076567/150‑1 (MQ=255)
  gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑Gc                                                                                                                                                 <  1:1136134/150‑1 (MQ=255)
  gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACACTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑ga                                                                                                                                                 >  1:1337097/1‑149 (MQ=255)
    gCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTg                                                                                                                                               >  1:400031/1‑150 (MQ=255)
    gCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTg                                                                                                                                               >  1:1599467/1‑150 (MQ=255)
      cATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGc                                                                                                                                             >  1:734231/1‑150 (MQ=255)
         tAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAG‑CTGTATTCAGGCTGGCGtc                                                                                                                                          <  1:1266704/150‑2 (MQ=255)
           aCCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAG‑CTGTATTCAGGCTGGCGtcag                                                                                                                                        <  2:905093/150‑4 (MQ=255)
            ccACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTc                                                                                                                                       <  2:781057/150‑1 (MQ=255)
             cacaGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCa                                                                                                                                      >  1:362488/1‑150 (MQ=255)
                aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCAttt                                                                                                                                   <  2:1226918/150‑1 (MQ=255)
                aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTATGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCAttt                                                                                                                                   >  1:336197/1‑150 (MQ=255)
                 gAAACAGCACTGTTGTCCTTTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTta                                                                                                                                  <  1:1445201/150‑1 (MQ=255)
                     cAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATaaa                                                                                                                              >  1:1656120/1‑150 (MQ=255)
                      aGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAAc                                                                                                                             <  1:469386/150‑1 (MQ=255)
                       gCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACa                                                                                                                            <  2:70501/150‑1 (MQ=255)
                        cACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTCGCGTTTTCATTTATAAACAg                                                                                                                           >  1:616325/1‑150 (MQ=255)
                                 ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGACAGGaaaa                                                                                                                  >  2:1262727/1‑150 (MQ=255)
                                 ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGaaaa                                                                                                                  >  2:616847/1‑150 (MQ=255)
                                 ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTACATTTATAAACAGAAAGGaaaa                                                                                                                  >  1:564602/1‑150 (MQ=255)
                                  cTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAAc                                                                                                                 >  1:1732662/1‑150 (MQ=255)
                                  cTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAAc                                                                                                                 <  1:1726376/150‑1 (MQ=255)
                                   tGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAAcc                                                                                                                >  2:74113/1‑150 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                           <  1:1357966/150‑1 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                           >  1:290348/1‑150 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAAGCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                           <  2:48224/150‑1 (MQ=255)
                                         aaGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataat                                                                                                          >  1:659799/1‑150 (MQ=255)
                                          aGCAGGCAATAGAATATGTCCTTACCCTGCGTTATCTGCAGCTTATGAACAAGAACAGTTGAATCTCAATCATCAGGCTTCATTTCTCGTGAGGCCAGCGTGATTACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGATAACCGataata                                                                                                         <  2:1669373/150‑1 (MQ=255)
                                          aGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataata                                                                                                         <  1:901563/150‑1 (MQ=255)
                                           gCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAtaataa                                                                                                        <  2:1018346/150‑1 (MQ=255)
                                             aGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatg                                                                                                      >  1:629579/1‑150 (MQ=255)
                                              ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                     >  1:800651/1‑150 (MQ=255)
                                              ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                     <  2:414409/150‑1 (MQ=255)
                                                  agagAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatg                                                                                                 >  2:761920/1‑150 (MQ=255)
                                                   gagaAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATAGA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatgg                                                                                                >  2:432652/1‑150 (MQ=255)
                                                    agaagaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggattggc                                                                                               >  1:315665/1‑149 (MQ=255)
                                                       agaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGAGAATCTCACTCAACAGGCATCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGc                                                                                            <  1:59506/150‑1 (MQ=255)
                                                        gaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGcc                                                                                           <  1:500379/150‑1 (MQ=255)
                                                        gaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCAGAATACA‑GCAGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAAAATAATGGAATGGACGcc                                                                                           >  2:1632666/1‑150 (MQ=255)
                                                          tGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGcccc                                                                                         >  1:525870/1‑150 (MQ=255)
                                                                tACCCTGCGCCATCTGCAGCTTCTGCCCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAt                                                                                   >  1:105655/1‑150 (MQ=255)
                                                                tACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAt                                                                                   >  1:101159/1‑150 (MQ=255)
                                                                   ccTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGaa                                                                                >  2:913448/1‑150 (MQ=255)
                                                                    cTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAg                                                                               <  1:1300649/150‑1 (MQ=255)
                                                                         ccATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACg                                                                          >  2:515029/1‑150 (MQ=255)
                                                                          cATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAACACCGAGAATAATGGAATGGACGCCCCGCCCATGAAGCAACgg                                                                         >  1:714090/1‑150 (MQ=255)
                                                                             cTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAg                                                                      >  1:485331/1‑150 (MQ=255)
                                                                              tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                     >  1:1534147/1‑150 (MQ=255)
                                                                              tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                     >  2:687601/1‑150 (MQ=255)
                                                                                       tGCACAAGAAACGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGAGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGTAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTg                                                            <  1:844996/150‑1 (MQ=255)
                                                                                          aCCAGAACCGTTGAATCTAACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         <  2:1442364/150‑1 (MQ=255)
                                                                                          aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         >  1:1176167/1‑150 (MQ=255)
                                                                                          aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         <  1:761920/150‑1 (MQ=255)
                                                                                            cAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGAaaaa                                                       >  2:792293/1‑150 (MQ=255)
                                                                                             aGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAc                                                      >  1:1579793/1‑150 (MQ=255)
                                                                                              gAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAcc                                                     <  1:1075415/150‑1 (MQ=255)
                                                                                              gAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAcc                                                     <  2:659799/150‑1 (MQ=255)
                                                                                                aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                                   >  2:1172509/1‑150 (MQ=255)
                                                                                                aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCACGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                                   >  1:540755/1‑150 (MQ=255)
                                                                                                 ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:1587654/1‑150 (MQ=255)
                                                                                                 ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTATTCATTTATAAACAGAAAGGAAAACCGATAAAAATGGAATGGACGCCCCGACCAAGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:698572/1‑150 (MQ=255)
                                                                                                  cGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAATGGTAGGATATTACTGATAAACCTTcc                                                 >  2:1454666/1‑150 (MQ=255)
                                                                                                  cGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTcc                                                 <  1:337031/150‑1 (MQ=255)
                                                                                                   ggTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCa                                                <  1:1281657/150‑1 (MQ=255)
                                                                                                    gTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGTTAGGATATTACTGATAAACCTTCCAg                                               <  1:74113/150‑1 (MQ=255)
                                                                                                    gTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTATCATTTATAAACAGAAAGGAAAACGGATAATAATGGAATGGACGCACCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAg                                               <  2:270155/150‑1 (MQ=255)
                                                                                                         tctcACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAg                                          <  1:1731580/150‑1 (MQ=255)
                                                                                                         tctcACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAg                                          <  2:1337097/150‑1 (MQ=255)
                                                                                                          atcactcaACATGCTTCATTTCTCCTGACGCCTGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAgg                                         <  1:1294142/149‑1 (MQ=255)
                                                                                                            cactcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAGGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTg                                       >  2:653074/1‑150 (MQ=255)
                                                                                                                 aaCAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTa                                  <  2:564602/150‑1 (MQ=255)
                                                                                                                   cAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATg                                >  2:1082616/1‑150 (MQ=255)
                                                                                                                    aGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATgg                               <  1:170705/150‑1 (MQ=255)
                                                                                                                      gCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTTCGTTATGGcc                             >  2:443747/1‑150 (MQ=255)
                                                                                                                       cTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCg                            <  2:1552667/150‑1 (MQ=255)
                                                                                                                        ttCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGt                           >  1:1606419/1‑150 (MQ=255)
                                                                                                                         tCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGCCGCCCCGCCCATGAAGCAACGGTAGGTTATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGtt                          >  2:1143814/1‑150 (MQ=255)
                                                                                                                           aTTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGtttt                        <  2:800651/150‑1 (MQ=255)
                                                                                                                              tctcCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCTTGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     <  2:1695335/150‑1 (MQ=255)
                                                                                                                              tctcCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     >  1:470491/1‑150 (MQ=255)
                                                                                                                              aatgaaaaCGCCAG‑CTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAACCCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACATTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     >  1:117715/8‑150 (MQ=255)
                                                                                                                              aatgaaaaCGCCAG‑CTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     >  1:1592232/8‑150 (MQ=255)
                                                                                                                               ctcCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCAAGTGGAGGTGCGTTATGGCCGTTTTAACg                    >  2:132786/1‑150 (MQ=255)
                                                                                                                                 gaaaaCGCCAG‑CTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGCCGCCCCGACCATGAAGCAACGGTAGGATATTAATGATAAACCTTCCAGTGGAGGTGCGGTATGGCCGTATAAACGAt                  >  2:333398/5‑150 (MQ=255)
                                                                                                                                 ccTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt                  <  2:422596/150‑1 (MQ=255)
                                                                                                                                 ccTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAGGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt                  >  1:811990/1‑150 (MQ=255)
                                                                                                                                   tGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                <  1:569892/150‑1 (MQ=255)
                                                                                                                                   tGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGAAGCCCCGACCATGAAGAAACGGTAGGATAATACTGATTAACATTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                >  2:1145460/1‑150 (MQ=255)
                                                                                                                                   aaaCGCCAG‑CTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                >  1:589516/3‑150 (MQ=255)
                                                                                                                                    gACGCCAGCCTGAATAAA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAg               <  1:432652/150‑1 (MQ=255)
                                                                                                                                     aCGCCAG‑CTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGTATATTACTGATAAACCTTCCAGTCGAGGAGCGATATGGCCCTTTTAACGATTAGc              >  1:191356/1‑150 (MQ=255)
                                                                                                                                      cGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa             <  2:105655/150‑1 (MQ=255)
                                                                                                                                      cGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa             <  2:101159/150‑1 (MQ=255)
                                                                                                                                       gCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGCCCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCAt            >  2:784976/1‑150 (MQ=255)
                                                                                                                                            ccTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCATCGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGAcc       <  1:913448/150‑1 (MQ=255)
                                                                                                                                            ccTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTAATGGCCGTTTTAACGATTAGCATTGAcc       >  1:977805/1‑150 (MQ=255)
                                                                                                                                              tGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTg     >  1:542600/1‑150 (MQ=255)
                                                                                                                                              tGAATAAA‑GCTTGCGTTTTCATTTATAAACAGAAAGTAAAAAATATAATAATGGAATGGACGCCCCGACCATGAATGAACGGTAGGATATAACTGATAAACCTTCCAGTGGAGGTGAGTTATGGCCGTTTTAACGATTAGCATTGACCTg     <  1:1252731/150‑1 (MQ=255)
                                                                                                                                                aaTACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGAc   <  1:1032656/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACA‑GCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACA  >  NC_000913/2078157‑2078451

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: ‑