Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,078,297 Δ1 bp 6.2% intergenic (+190/+252) 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,2970C. 6.2% 342.9 / 16.8 113intergenic (+190/+252)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 C (59/47);  new base . (5/2);  total (64/49)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 116 total aligned reads displayed.

CAGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACC  >  NC_000913/2078150‑2078446
                                                                                                                                                   |                                                                                                                                                      
cAGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCt                                                                                                                                                      >  2:569892/1‑150 (MQ=255)
cAGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGACAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGctg                                                                                                                                                      >  1:16904/1‑148 (MQ=255)
 aGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTg                                                                                                                                                     >  1:627395/1‑150 (MQ=255)
 aGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTg                                                                                                                                                     <  2:1373865/150‑1 (MQ=255)
   ggTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGaa                                                                                                                                                   <  2:326254/150‑1 (MQ=255)
     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)
         gTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACACTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAga                                                                                                                                            >  1:1337097/1‑149 (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)
                                        ccTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTACATTTATAAACAGAAAGGaaaa                                                                                                             >  1:564602/1‑150 (MQ=255)
                                         cTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAc                                                                                                            <  1:1726376/150‑1 (MQ=255)
                                          tGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAAcc                                                                                                           >  2:74113/1‑150 (MQ=255)
                                               tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      >  1:290348/1‑150 (MQ=255)
                                               tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      <  1:1357966/150‑1 (MQ=255)
                                               tAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAAGCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataa                                                                                                      <  2:48224/150‑1 (MQ=255)
                                                aaGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGataat                                                                                                     >  1:659799/1‑150 (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)
                                                           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)
                                                                                  aTCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGt                                                                   >  2:560764/1‑150 (MQ=255)
                                                                                    cTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAg                                                                 >  1:485331/1‑150 (MQ=255)
                                                                                     tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                >  2:687601/1‑150 (MQ=255)
                                                                                     tGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAgg                                                                >  1:1534147/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)
                                                                                                 aCCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATa                                                    >  1:1176167/1‑150 (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:1587654/1‑150 (MQ=255)
                                                                                                        ccGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTc                                             >  1:1627242/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)
                                                                                                                 atcactcaACATGCTTCATTTCTCCTGACGCCTGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAgg                                    <  1:1294142/149‑1 (MQ=255)
                                                                                                                   cactcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCAGGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTg                                  >  2:653074/1‑150 (MQ=255)
                                                                                                                     ctcaACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCg                                <  1:560764/150‑1 (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)
                                                                                                                              cTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCg                       >  2:759198/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)
                                                                                                                                     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:331550/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)
                                                                                                                                                   |                                                                                                                                                      
CAGGGTATAGTGCCATTAACCACAGAAACAGCACTGTTGTCCTGTATTAAGCAGGCAAGAGAAGATGTCCTTACCCTGCGCCATCTGCAGCTTCTGCACCAGAACCGGTGAATCTCACTCAACAGGCTTCATTTCTCCTGACGCCAGCCTGAATACAGCTGGCGTTTTCATTTATAAACAGAAAGGAAAACCGATAATAATGGAATGGACGCCCCGACCATGAAGCAACGGTAGGATATTACTGATAAACCTTCCAGTGGAGGTGCGTTATGGCCGTTTTAACGATTAGCATTGACC  >  NC_000913/2078150‑2078446

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: