Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,581 T→A 7.4% intergenic (+174/‑325) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,5810TA 7.4% 329.4 / 11.7 95intergenic (+174/‑325)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (40/48);  new base A (2/5);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATG  >  NC_000913/4235440‑4235729
                                                                                                                                             |                                                                                                                                                    
cACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCgg                                                                                                                                               >  2:1448603/1‑149 (MQ=255)
   gccgcATCTGGCAAGCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGTCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGCTTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAt                                                                                                                                            <  4:898109/149‑1 (MQ=255)
   gccgcATCCGGGAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAt                                                                                                                                            >  3:1146283/1‑149 (MQ=255)
   gccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAt                                                                                                                                            >  2:168854/1‑149 (MQ=255)
   gccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAt                                                                                                                                            <  3:1090469/149‑1 (MQ=255)
   gccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAt                                                                                                                                            <  4:1361072/149‑1 (MQ=255)
    ccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                                                                           <  3:1109006/149‑1 (MQ=255)
      gcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttt                                                                                                                                         >  1:839795/1‑149 (MQ=255)
        aTCAGGCATCGGTTGCCGGATGCGGCGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGCTTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttttt                                                                                                                                       >  1:388279/1‑149 (MQ=255)
        aTCAGGCATCGGTTGCCGGATGCGACGCGGTCGCGTCTTATCAGGCCTACACGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttttt                                                                                                                                       >  1:1012818/1‑149 (MQ=255)
         tCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttttt                                                                                                                                      <  3:1046011/149‑1 (MQ=255)
          ccGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTa                                                                                                                                     <  4:1416886/149‑1 (MQ=255)
          ccGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTa                                                                                                                                     <  1:458557/149‑1 (MQ=255)
           cGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTAt                                                                                                                                    >  3:841818/1‑149 (MQ=255)
             gCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACGTAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCa                                                                                                                                  <  1:1152141/149‑1 (MQ=255)
                         tGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCt                                                                                                                      >  3:309965/1‑149 (MQ=255)
                          gataagACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCa                                                                                                                                  <  1:1112577/131‑1 (MQ=255)
                            taagACGCGACCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                   >  4:794175/4‑149 (MQ=255)
                            tGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                   <  3:699970/149‑1 (MQ=255)
                                aCGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGt                                                                                                               <  3:1369559/149‑1 (MQ=255)
                                  gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa                                                                                                             <  2:99204/149‑1 (MQ=255)
                                  gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGCCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa                                                                                                             <  1:259235/149‑1 (MQ=255)
                                        gcgTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAAc                                                                                                       <  4:825037/149‑1 (MQ=255)
                                         cgTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACg                                                                                                      <  1:346501/149‑1 (MQ=255)
                                          gTCTTATCAGGCCTACAGGTCGATGCCGATATGTACCTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTAATGTCGGCATTTTTTTATCAGATAAATCCCCGTGTCTGTAATTTAACgg                                                                                                     >  4:1230740/1‑149 (MQ=255)
                                            cTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                   <  4:1365622/133‑1 (MQ=255)
                                            cTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                   >  3:1365622/1‑133 (MQ=255)
                                                tCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTc                                                                                                                  <  3:1456021/130‑1 (MQ=255)
                                                tCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTc                                                                                                                  >  4:1456021/1‑130 (MQ=255)
                                                  aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGa                                                                                                                                <  4:1453463/114‑1 (MQ=255)
                                                  aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGa                                                                                                                                >  3:1453463/1‑114 (MQ=255)
                                                  aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                             >  1:868646/1‑149 (MQ=255)
                                                  aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                             <  1:168854/149‑1 (MQ=255)
                                                  aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                             <  2:1297502/149‑1 (MQ=255)
                                                  aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                             <  1:172509/149‑1 (MQ=255)
                                                   ggCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            <  1:403133/149‑1 (MQ=255)
                                                   ggCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            >  2:1418228/1‑149 (MQ=255)
                                                   ggCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            <  1:821578/149‑1 (MQ=255)
                                                       tACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            <  2:226744/145‑1 (MQ=255)
                                                       tACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            >  1:226744/1‑145 (MQ=255)
                                                       tACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTg                                                                                        >  3:1397359/1‑149 (MQ=255)
                                                               gATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAAt                                                                                >  3:1387681/1‑149 (MQ=255)
                                                                aTGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATc                                                                               <  1:1330987/149‑1 (MQ=255)
                                                                    cGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaa                                                                           >  4:622753/1‑149 (MQ=255)
                                                                      atatGTACATCGAATTCAACAATTAATACCTAACCCGATTGATTAAATAACCTTAATAACAATGCTGACGTTATGTCGGCATTTTTTTATCAGATAAGGCCCCTTGGCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaaaa                                                                         <  2:1012818/149‑1 (MQ=255)
                                                                       tatGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAAt                                                                        <  3:1334808/149‑1 (MQ=255)
                                                                        atGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATa                                                                       <  2:199337/149‑1 (MQ=255)
                                                                          gTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGa                                                                     <  4:81408/149‑1 (MQ=255)
                                                                          gTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGa                                                                     <  1:1448603/149‑1 (MQ=255)
                                                                           tACATCGTATTCGGCAATTAATATATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAt                                                                    >  4:728092/1‑149 (MQ=255)
                                                                           tACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAGCCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAt                                                                    <  1:1205800/149‑1 (MQ=255)
                                                                            aCATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAtt                                                                   >  4:1287269/1‑149 (MQ=255)
                                                                              aTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAtttt                                                                 <  2:839795/149‑1 (MQ=255)
                                                                              aTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAtttt                                                                 >  1:607701/1‑149 (MQ=255)
                                                                               tCGTATTCGGCAACTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAttttt                                                                <  4:117026/149‑1 (MQ=255)
                                                                                       ggCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgttg                                                        >  2:578015/1‑149 (MQ=255)
                                                                                        gCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgttgt                                                       <  4:252695/149‑1 (MQ=255)
                                                                                        gCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgttgt                                                       <  3:1235505/149‑1 (MQ=255)
                                                                                          aaTTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGttt                                                     <  3:622753/149‑1 (MQ=255)
                                                                                            ttAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTAt                                                   >  4:1069159/1‑149 (MQ=255)
                                                                                                tACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCg                                               >  1:925972/1‑149 (MQ=255)
                                                                                                tACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCg                                               >  2:393116/1‑149 (MQ=255)
                                                                                                   aTAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGATATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTc                                                <  1:1387676/145‑1 (MQ=255)
                                                                                                   aTAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTc                                                >  2:1387676/1‑145 (MQ=255)
                                                                                                     aGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt                                                                                                                     >  1:641560/1‑74 (MQ=255)
                                                                                                     aGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt                                                                                                                     <  2:641560/74‑1 (MQ=255)
                                                                                                        acagttgattAAATAAACTTAATAAAAAAGCCGACGGTATGTTCGGATTTTTTTTACAGATAAATCCCCTTTGCTTTTATTTTACGGAAATTATTCCCTTAGGGTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATAc                                       <  1:618451/145‑1 (MQ=255)
                                                                                                        aCGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaaaa                                                                         >  3:270099/1‑115 (MQ=255)
                                                                                                        aCGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaaaa                                                                         <  4:270099/115‑1 (MQ=255)
                                                                                                                ttAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            >  3:1446257/1‑88 (MQ=255)
                                                                                                                ttAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATAc                                                                                            <  4:1446257/88‑1 (MQ=255)
                                                                                                                ttAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATAc                                       >  4:432543/1‑141 (MQ=255)
                                                                                                                ttAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATAc                                       <  3:432543/141‑1 (MQ=255)
                                                                                                                ttAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTa                               >  3:488205/1‑149 (MQ=255)
                                                                                                                   aaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATc                                                                                                                                   <  1:653916/46‑1 (MQ=255)
                                                                                                                   aaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATc                                                                                                                                   >  2:653916/1‑46 (MQ=255)
                                                                                                                   aaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAAtt                            <  4:488205/149‑1 (MQ=255)
                                                                                                                    aTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTg                           >  2:264077/1‑149 (MQ=255)
                                                                                                                       aCCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcccc                                                                                                                       <  4:864525/54‑1 (MQ=255)
                                                                                                                       aCCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcccc                                                                                                                       >  3:864525/1‑54 (MQ=255)
                                                                                                                         cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                   >  3:1340347/1‑56 (MQ=255)
                                                                                                                         cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                   <  4:1340347/56‑1 (MQ=255)
                                                                                                                         cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAttttt                                                                <  1:707738/107‑1 (MQ=255)
                                                                                                                         cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAttttt                                                                >  2:707738/1‑107 (MQ=255)
                                                                                                                         cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCtttt                      >  2:173159/1‑149 (MQ=255)
                                                                                                                          ttAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTg                     >  3:1198034/1‑149 (MQ=255)
                                                                                                                           taataaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGt                    >  4:294106/1‑149 (MQ=255)
                                                                                                                             atcaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGTAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGttt                  <  2:388279/146‑1 (MQ=255)
                                                                                                                             ataaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGttt                  <  1:1308583/149‑1 (MQ=255)
                                                                                                                                aCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTtttta               <  1:500442/149‑1 (MQ=255)
                                                                                                                                 cAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttat              <  3:794175/149‑1 (MQ=255)
                                                                                                                                  aaTGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttatt             >  3:1183327/1‑149 (MQ=255)
                                                                                                                                   aTGCCGACATAACGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTtttta               <  2:230085/146‑1 (MQ=255)
                                                                                                                                   aTGCCGACATAACGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTtttta               >  1:230085/1‑146 (MQ=255)
                                                                                                                                      ccGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTaa         <  1:1464954/149‑1 (MQ=255)
                                                                                                                                            ttATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAt   >  3:571614/1‑149 (MQ=255)
                                                                                                                                            ttATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAt   <  4:309965/149‑1 (MQ=255)
                                                                                                                                            ttATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAt   <  2:607701/149‑1 (MQ=255)
                                                                                                                                            ttATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAt   >  2:861492/1‑149 (MQ=255)
                                                                                                                                             tATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATg  <  4:1146283/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATG  >  NC_000913/4235440‑4235729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: