Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,078,304 A→T 6.1% intergenic (+197/+245) 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,3040AT 6.1% 316.7 / 7.7 114intergenic (+197/+245)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 A (59/48);  new base T (5/2);  total (64/50)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 115 total aligned reads displayed.

TATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACA  >  NC_000913/2078155‑2078451
                                                                                                                                                     |                                                                                                                                                    
tataGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTAACCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATa                                                                                                                                                      >  1:1146739/1‑150 (MQ=255)
  aaGTGCCATAAACCACAGAAACAGCACTGTTGTCCTGTATAAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCTCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACa                                                                                                                                                    <  1:216343/149‑1 (MQ=255)
   aGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAg                                                                                                                                                  <  1:1076567/150‑1 (MQ=255)
    gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGc                                                                                                                                                 <  1:1136134/150‑1 (MQ=255)
      gCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTg                                                                                                                                               >  1:400031/1‑150 (MQ=255)
      gCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTg                                                                                                                                               >  1:1599467/1‑150 (MQ=255)
        cATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGCATCTCCCTCAACAGGGTTCATTTCTCCTGACGCAAGCCTGAATACAGCTGGc                                                                                                                                             >  2:13716/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)
                  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)
                                 gtCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGAGTTTTCATTTATAAACAGAAAGGaa                                                                                                                    >  1:862702/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: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                                                                                                     <  2:414409/150‑1 (MQ=255)
                                                ggCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatgg                                                                                                     >  1:800651/1‑150 (MQ=255)
                                                     gagaAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATAGAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggaatgg                                                                                                >  2:432652/1‑150 (MQ=255)
                                                      agaagaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATaatggattggc                                                                                               >  1:315665/1‑149 (MQ=255)
                                                          gaTGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCAGAATACAGCAGGCGTTTTCATTTATAAACAGAAAGGAAAACCGAAAATAATGGAATGGACGcc                                                                                           >  2:1632666/1‑150 (MQ=255)
                                                            tGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGcccc                                                                                         >  1:525870/1‑150 (MQ=255)
                                                             gTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGCATCTAAATCAACAGGCTTCTTTTCTCCTGACGCCGGCCTGATGACCGCTGGCGTTTTGTGGTATAAACTGAAAGGAAAACCGATAATAATGGAATGGACGcccct                                                                                        >  1:251228/1‑149 (MQ=255)
                                                                cTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGAcc                                                                                     >  1:663467/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)
                                                                             aTCTTCAGATTCTGCACCAGAACCGTTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCATCCTGAATACAGCTGGCTTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATTGACGCCCCGACCATGAAGCAACGGt                                                                        <  1:13716/150‑1 (MQ=255)
                                                                             aTCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGt                                                                        >  2:560764/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)
                                                                                      ttCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTa                                                               <  1:449074/150‑1 (MQ=255)
                                                                                         tGCACAAGAAACGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGAGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGTAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTg                                                            <  1:844996/150‑1 (MQ=255)
                                                                                            aCCAGAACCGTTGAATCTAACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         <  2:1442364/150‑1 (MQ=255)
                                                                                            aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                         <  1:761920/150‑1 (MQ=255)
                                                                                              cAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGAaaaa                                                       >  2:792293/1‑150 (MQ=255)
                                                                                               aGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAc                                                      >  1:1579793/1‑150 (MQ=255)
                                                                                                gAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAcc                                                     <  2:659799/150‑1 (MQ=255)
                                                                                                gAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAAcc                                                     <  1:1075415/150‑1 (MQ=255)
                                                                                                  aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                                   >  2:1172509/1‑150 (MQ=255)
                                                                                                  aCCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCACGACCATGAAGCAACGGTAGGATATTACTGATAAACCtt                                                   >  1:540755/1‑150 (MQ=255)
                                                                                                   ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:1627242/1‑150 (MQ=255)
                                                                                                   ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:1587654/1‑150 (MQ=255)
                                                                                                   ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTATTCATTTATAAACAGAAAGGAAAACCGATAAAAATGGAATGGACGCCCCGACCAAGAAGCAACGGTAGGATATTACTGATAAACCTTc                                                  >  1:698572/1‑150 (MQ=255)
                                                                                                    cGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTcc                                                 <  1:337031/150‑1 (MQ=255)
                                                                                                     ggTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCa                                                <  1:1281657/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)
                                                                                                                ctcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCg                                     <  1:560764/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)
                                                                                                                         cTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCg                            >  2:759198/1‑150 (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)
                                                                                                                                 ctcCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCAAGTGGAGGTGCGTTATGGCCGTTTTAACg                    >  2:132786/1‑150 (MQ=255)
                                                                                                                                 ctcCTGACGCCACCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACg                    <  2:1124966/150‑1 (MQ=255)
                                                                                                                                   ccTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt                  <  2:422596/150‑1 (MQ=255)
                                                                                                                                   ccTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAGGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGAt                  >  1:811990/1‑150 (MQ=255)
                                                                                                                                     tGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTa                <  1:569892/150‑1 (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:101159/150‑1 (MQ=255)
                                                                                                                                        cGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa             <  2:105655/150‑1 (MQ=255)
                                                                                                                                        cGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCa             >  2:331550/1‑150 (MQ=255)
                                                                                                                                         gCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGCCCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCAt            >  2:784976/1‑150 (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)
                                                                                                                                                   aTACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACa  >  1:1231807/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACCTGACA  >  NC_000913/2078155‑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: