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.

TAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACA  >  NC_000913/2078157‑2078451
                                                                                                                                                      |                                                                                                                                                 
aaGTGCCATAAACCACAGAAACAGCACTGTTGTCCTGTATAAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCTCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACa                                                                                                                                                    <  1:216343/149‑1 (MQ=255)
 aGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAg                                                                                                                                                  <  1:1076567/150‑1 (MQ=255)
  gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGc                                                                                                                                                 <  1:1136134/150‑1 (MQ=255)
  gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACACTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAga                                                                                                                                                 >  1:1337097/1‑149 (MQ=255)
    gCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTg                                                                                                                                               >  1:400031/1‑150 (MQ=255)
    gCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTg                                                                                                                                               >  1:1599467/1‑150 (MQ=255)
      cATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGc                                                                                                                                             >  1:734231/1‑150 (MQ=255)
         tAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCTGTATTCAGGCTGGCGtc                                                                                                                                          <  1:1266704/150‑2 (MQ=255)
           aCCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCTGTATTCAGGCTGGCGtcag                                                                                                                                        <  2:905093/150‑4 (MQ=255)
            ccACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTc                                                                                                                                       <  2:781057/150‑1 (MQ=255)
             cacaGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCa                                                                                                                                      >  1:362488/1‑150 (MQ=255)
                aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCAttt                                                                                                                                   <  2:1226918/150‑1 (MQ=255)
                aGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTATGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCAttt                                                                                                                                   >  1:336197/1‑150 (MQ=255)
                 gAAACAGCACTGTTGTCCTTTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTta                                                                                                                                  <  1:1445201/150‑1 (MQ=255)
                     cAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATaaa                                                                                                                              >  1:1656120/1‑150 (MQ=255)
                      aGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAAc                                                                                                                             <  1:469386/150‑1 (MQ=255)
                       gCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACa                                                                                                                            <  2:70501/150‑1 (MQ=255)
                        cACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTCGCGTTTTCATTTATAAACAg                                                                                                                           >  1:616325/1‑150 (MQ=255)
                                 ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGACAGGaaaa                                                                                                                  >  2:1262727/1‑150 (MQ=255)
                                 ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGaaaa                                                                                                                  >  2:616847/1‑150 (MQ=255)
                                 ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTACATTTATAAACAGAAAGGaaaa                                                                                                                  >  1:564602/1‑150 (MQ=255)
                                  cTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAc                                                                                                                 >  1:1732662/1‑150 (MQ=255)
                                  cTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAc                                                                                                                 <  1:1726376/150‑1 (MQ=255)
                                   tGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAcc                                                                                                                >  2:74113/1‑150 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                           <  1:1357966/150‑1 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                           >  1:290348/1‑150 (MQ=255)
                                        tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAAGCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                           <  2:48224/150‑1 (MQ=255)
                                         aaGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataat                                                                                                          >  1:659799/1‑150 (MQ=255)
                                          aGCAGGCAATAGAATATGTCCTTACCCTGCGTTATCTGCAGCTTATGAACAAGAACAGTTGAATCTCAATCATCAGGCTTCATTTCTCGTGAGGCCAGCGTGATTACAGCTGGCGTTTTCATTTATAAACAGAAAGGATAACCGataata                                                                                                         <  2:1669373/150‑1 (MQ=255)
                                          aGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataata                                                                                                         <  1:901563/150‑1 (MQ=255)
                                           gCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAtaataa                                                                                                        <  2:1018346/150‑1 (MQ=255)
                                             aGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatg                                                                                                      >  1:629579/1‑150 (MQ=255)
                                              ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                     >  1:800651/1‑150 (MQ=255)
                                              ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                     <  2:414409/150‑1 (MQ=255)
                                                  agagAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatg                                                                                                 >  2:761920/1‑150 (MQ=255)
                                                   gagaAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATAGAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatgg                                                                                                >  2:432652/1‑150 (MQ=255)
                                                    agaagaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggattggc                                                                                               >  1:315665/1‑149 (MQ=255)
                                                       agaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGAGAATCTCACTCAACAGGCATCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGc                                                                                            <  1:59506/150‑1 (MQ=255)
                                                        gaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGcc                                                                                           <  1:500379/150‑1 (MQ=255)
                                                        gaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCAGAATACAGCAGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAAAATAATGGAATGGACGcc                                                                                           >  2:1632666/1‑150 (MQ=255)
                                                          tGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGcccc                                                                                         >  1:525870/1‑150 (MQ=255)
                                                                tACCCTGCGCCATCTGCAGCTTCTGCCCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAt                                                                                   >  1:105655/1‑150 (MQ=255)
                                                                tACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAt                                                                                   >  1:101159/1‑150 (MQ=255)
                                                                   ccTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGaa                                                                                >  2:913448/1‑150 (MQ=255)
                                                                    cTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAg                                                                               <  1:1300649/150‑1 (MQ=255)
                                                                         ccATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACg                                                                          >  2:515029/1‑150 (MQ=255)
                                                                          cATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAACACCGAGAATAATGGAATGGACGCCCCGCCCATGAAGCAACgg                                                                         >  1:714090/1‑150 (MQ=255)
                                                                             cTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAg                                                                      >  1:485331/1‑150 (MQ=255)
                                                                              tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                     >  1:1534147/1‑150 (MQ=255)
                                                                              tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                     >  2:687601/1‑150 (MQ=255)
                                                                                       tGCACAAGAAACGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGAGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGTAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTg                                                            <  1:844996/150‑1 (MQ=255)
                                                                                          aCCAGAACCGTTGAATCTAACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         <  2:1442364/150‑1 (MQ=255)
                                                                                          aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         >  1:1176167/1‑150 (MQ=255)
                                                                                          aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         <  1:761920/150‑1 (MQ=255)
                                                                                            cAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGAaaaa                                                       >  2:792293/1‑150 (MQ=255)
                                                                                             aGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAc                                                      >  1:1579793/1‑150 (MQ=255)
                                                                                              gAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAcc                                                     <  1:1075415/150‑1 (MQ=255)
                                                                                              gAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAcc                                                     <  2:659799/150‑1 (MQ=255)
                                                                                                aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                                   >  2:1172509/1‑150 (MQ=255)
                                                                                                aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCACGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                                   >  1:540755/1‑150 (MQ=255)
                                                                                                 ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:1587654/1‑150 (MQ=255)
                                                                                                 ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTATTCATTTATAAACAGAAAGGAAAACCGATAAAAATGGAATGGACGCCCCGACCAAGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:698572/1‑150 (MQ=255)
                                                                                                  cGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAATGGTAGGATATTACTGATAAACCTTcc                                                 >  2:1454666/1‑150 (MQ=255)
                                                                                                  cGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTcc                                                 <  1:337031/150‑1 (MQ=255)
                                                                                                   ggTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCa                                                <  1:1281657/150‑1 (MQ=255)
                                                                                                    gTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGTTAGGATATTACTGATAAACCTTCCAg                                               <  1:74113/150‑1 (MQ=255)
                                                                                                    gTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTATCATTTATAAACAGAAAGGAAAACGGATAATAATGGAATGGACGCACCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAg                                               <  2:270155/150‑1 (MQ=255)
                                                                                                         tctcACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAg                                          <  1:1731580/150‑1 (MQ=255)
                                                                                                         tctcACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAg                                          <  2:1337097/150‑1 (MQ=255)
                                                                                                          atcactcaACATGCTTCATTTCTCCTGACGCCTGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAgg                                         <  1:1294142/149‑1 (MQ=255)
                                                                                                            cactcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAGGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTg                                       >  2:653074/1‑150 (MQ=255)
                                                                                                                 aaCAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTa                                  <  2:564602/150‑1 (MQ=255)
                                                                                                                   cAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATg                                >  2:1082616/1‑150 (MQ=255)
                                                                                                                    aGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATgg                               <  1:170705/150‑1 (MQ=255)
                                                                                                                      gCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTTCGTTATGGcc                             >  2:443747/1‑150 (MQ=255)
                                                                                                                       cTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCg                            <  2:1552667/150‑1 (MQ=255)
                                                                                                                        ttCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGt                           >  1:1606419/1‑150 (MQ=255)
                                                                                                                         tCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGCCGCCCCGCCCATGAAGCAACGGTAGGTTATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGtt                          >  2:1143814/1‑150 (MQ=255)
                                                                                                                           aTTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGtttt                        <  2:800651/150‑1 (MQ=255)
                                                                                                                              tctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCTTGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     <  2:1695335/150‑1 (MQ=255)
                                                                                                                              tctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     >  1:470491/1‑150 (MQ=255)
                                                                                                                              aatgaaaaCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAACCCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACATTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     >  1:117715/8‑150 (MQ=255)
                                                                                                                              aatgaaaaCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAAc                     >  1:1592232/8‑150 (MQ=255)
                                                                                                                               ctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCAAGTGGAGGTGCGTTATGGCCGTTTTAACg                    >  2:132786/1‑150 (MQ=255)
                                                                                                                                 gaaaaCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGCCGCCCCGACCATGAAGCAACGGTAGGATATTAATGATAAACCTTCCAGTGGAGGTGCGGTATGGCCGTATAAACGAt                  >  2:333398/5‑150 (MQ=255)
                                                                                                                                 ccTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt                  <  2:422596/150‑1 (MQ=255)
                                                                                                                                 ccTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAGGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt                  >  1:811990/1‑150 (MQ=255)
                                                                                                                                   tGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                <  1:569892/150‑1 (MQ=255)
                                                                                                                                   tGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGAAGCCCCGACCATGAAGAAACGGTAGGATAATACTGATTAACATTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                >  2:1145460/1‑150 (MQ=255)
                                                                                                                                   aaaCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                >  1:589516/3‑150 (MQ=255)
                                                                                                                                    gACGCCAGCCTGAATAAAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAg               <  1:432652/150‑1 (MQ=255)
                                                                                                                                     aCGCCAGCTGTATTCAGGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGTATATTACTGATAAACCTTCCAGTCGAGGAGCGATATGGCCCTTTTAACGATTAGc              >  1:191356/1‑150 (MQ=255)
                                                                                                                                      cGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa             <  2:105655/150‑1 (MQ=255)
                                                                                                                                      cGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa             <  2:101159/150‑1 (MQ=255)
                                                                                                                                       gCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGCCCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCAt            >  2:784976/1‑150 (MQ=255)
                                                                                                                                            ccTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCATCGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGAcc       <  1:913448/150‑1 (MQ=255)
                                                                                                                                            ccTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTAATGGCCGTTTTAACGATTAGCATTGAcc       >  1:977805/1‑150 (MQ=255)
                                                                                                                                              tGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTg     >  1:542600/1‑150 (MQ=255)
                                                                                                                                              tGAATAAAGCTTGCGTTTTCATTTATAAACAGAAAGTAAAAAATATAATAATGGAATGGACGCCCCGACCATGAATGAACGGTAGGATATAACTGATAAACCTTCCAGTGGAGGTGAGTTATGGCCGTTTTAACGATTAGCATTGACCTg     <  1:1252731/150‑1 (MQ=255)
                                                                                                                                                aaTACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGAc   <  1:1032656/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACA  >  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: