Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,839,321 C→G 6.9% intergenic (‑57/‑203) ascG ← / → ascF DNA‑binding transcriptional repressor AscG/beta‑glucoside specific PTS enzyme IIBC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,839,3210CG 6.9% 367.0 / 12.3 111intergenic (‑57/‑203)ascG/ascFDNA‑binding transcriptional repressor AscG/beta‑glucoside specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base C (51/52);  new base G (4/4);  total (55/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 112 total aligned reads displayed.

CTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCC  >  NC_000913/2839172‑2839468
                                                                                                                                                     |                                                                                                                                                   
ctcCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTc                                                                                                                                                     >  2:679238/1‑150 (MQ=255)
    tGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTg                                                                                                                                                 <  2:336721/150‑1 (MQ=255)
     ggCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCCTCCCTTTACCCGCAATCAACGCAGTGTACTGCACCGTTTGCCGCTTGTCCTTGc                                                                                                                                                >  2:1167106/1‑150 (MQ=255)
          aCGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAAt                                                                                                                                           <  1:625245/150‑1 (MQ=255)
           cGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATc                                                                                                                                          <  1:1420368/150‑1 (MQ=255)
           cGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATc                                                                                                                                          <  1:1304102/150‑1 (MQ=255)
            gTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACCCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAAAcg                                                                                                                                         >  1:158918/1‑150 (MQ=255)
              aGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                       <  1:1408592/150‑1 (MQ=255)
              aGCCATTACCTGAAAGCACGCGGGAAACGGGCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCAAGGACAATcggc                                                                                                                                       <  2:1034283/150‑1 (MQ=255)
               gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTTCCGATTGTCCTTGCACAATaggcg                                                                                                                                      >  1:171115/1‑150 (MQ=255)
               gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggcg                                                                                                                                      >  1:1006987/1‑150 (MQ=255)
               gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggcg                                                                                                                                      >  1:873269/1‑150 (MQ=255)
               gCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggcg                                                                                                                                      >  1:450361/1‑150 (MQ=255)
                ccATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATcggcgg                                                                                                                                     <  1:679238/150‑1 (MQ=255)
                 cATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCggg                                                                                                                                    <  2:386833/150‑1 (MQ=255)
                 cATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTAATCCCTTTACACGCAATCAACGCAGTGTACTGCCCCGTTTGCCGATTGTCCTTGCACAATCGGCggg                                                                                                                                    >  1:379433/1‑150 (MQ=255)
                  aTTACCTGAAAGCACGCGTGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGa                                                                                                                                   <  1:214348/150‑1 (MQ=255)
                      ccTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACCCGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaaaa                                                                                                                               >  2:156608/1‑150 (MQ=255)
                      ccTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGaaaaa                                                                                                                               >  2:1113415/1‑150 (MQ=255)
                        tGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAata                                                                                                                             >  1:1364364/1‑150 (MQ=255)
                         gAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAatat                                                                                                                            >  2:1024497/1‑150 (MQ=255)
                           aaGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTc                                                                                                                          <  2:1582251/150‑1 (MQ=255)
                           aaGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTc                                                                                                                          <  2:1416949/150‑1 (MQ=255)
                             gCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAg                                                                                                                        >  1:1176770/1‑150 (MQ=255)
                              cACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAgg                                                                                                                       <  2:690140/150‑1 (MQ=255)
                                cgcgGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTg                                                                                                                     >  1:1596847/1‑150 (MQ=255)
                                   gggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGAcc                                                                                                                  <  2:1006987/150‑1 (MQ=255)
                                    ggAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCg                                                                                                                 >  2:334208/1‑150 (MQ=255)
                                         cGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTc                                                                                                            <  2:1529258/150‑1 (MQ=255)
                                         cGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTc                                                                                                            >  2:503222/1‑150 (MQ=255)
                                          ggTCGCTTTTGAAACCCCGGCGCGATTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTGCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTca                                                                                                           <  2:444510/150‑1 (MQ=255)
                                           ggcGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGGCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcac                                                                                                          <  2:917128/148‑1 (MQ=255)
                                           gTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcac                                                                                                          <  2:207002/150‑1 (MQ=255)
                                           gTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcac                                                                                                          >  2:541846/1‑150 (MQ=255)
                                           gTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTcac                                                                                                          >  2:1545093/1‑150 (MQ=255)
                                                    gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTAGTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaa                                                                                                 <  2:595506/150‑1 (MQ=255)
                                                    gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaa                                                                                                 >  1:1254341/1‑150 (MQ=255)
                                                    gAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaa                                                                                                 >  2:112010/1‑150 (MQ=255)
                                                     aaaCCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATaaaa                                                                                                <  1:1009893/150‑1 (MQ=255)
                                                            ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                                         >  2:237672/1‑150 (MQ=255)
                                                            ggCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAAc                                                                                         >  2:373965/1‑150 (MQ=255)
                                                                     cGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActc                                                                                >  1:40751/1‑150 (MQ=255)
                                                                      gCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCActct                                                                               >  2:526336/1‑150 (MQ=255)
                                                                         aCTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCtt                                                                            <  2:1207293/150‑1 (MQ=255)
                                                                          cTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTg                                                                           <  2:1176770/150‑1 (MQ=255)
                                                                           ttCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGc                                                                          <  1:1036275/150‑1 (MQ=255)
                                                                           ttCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGc                                                                          <  2:624770/150‑1 (MQ=255)
                                                                            tCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCt                                                                         >  1:1521471/1‑150 (MQ=255)
                                                                             ccAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCtt                                                                        >  2:1281030/1‑150 (MQ=255)
                                                                                gCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTAtt                                                                     >  1:223235/1‑150 (MQ=255)
                                                                                 cATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTAttt                                                                    >  1:1557436/1‑150 (MQ=255)
                                                                                    cgtcgtCATCATTTTCATCCCTTTACACGCAATGAACGCAGTGTAGTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGt                                                                 <  2:158918/150‑1 (MQ=255)
                                                                                     gtcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTg                                                                >  2:1048241/1‑150 (MQ=255)
                                                                                     gtcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTg                                                                <  2:873269/150‑1 (MQ=255)
                                                                                      tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGa                                                               <  1:1256294/150‑1 (MQ=255)
                                                                                       cgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGAc                                                              <  1:106474/150‑1 (MQ=255)
                                                                                          catcatTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAAc                                                           <  1:1545093/150‑1 (MQ=255)
                                                                                              atTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTAtt                                                       <  2:1167671/150‑1 (MQ=255)
                                                                                               ttttCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTc                                                      >  2:502631/1‑150 (MQ=255)
                                                                                                tttCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCa                                                     <  2:441881/150‑1 (MQ=255)
                                                                                                 ttCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCAt                                                    <  1:531494/150‑1 (MQ=255)
                                                                                                 ttCATCCCTTTACACGCAATCAACGCACTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCAt                                                    <  1:323654/150‑1 (MQ=255)
                                                                                                  tCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATg                                                   <  1:886413/150‑1 (MQ=255)
                                                                                                  tCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATg                                                   <  2:1557436/150‑1 (MQ=255)
                                                                                                     tCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtt                                                >  2:1050769/1‑150 (MQ=255)
                                                                                                      cccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAAGATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAttt                                               >  2:1434106/1‑150 (MQ=255)
                                                                                                       ccTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtttt                                              >  1:418476/1‑150 (MQ=255)
                                                                                                          ttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                           <  1:522515/150‑1 (MQ=255)
                                                                                                          ttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                           <  2:40751/150‑1 (MQ=255)
                                                                                                          ttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTtgtg                                           >  1:556573/1‑150 (MQ=255)
                                                                                                                 cAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                    >  1:37972/1‑150 (MQ=255)
                                                                                                                 cAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                    >  1:92651/1‑150 (MQ=255)
                                                                                                                     cAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCCGGTGACCGGTTTCACAAATATAAAAAATGGACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGtgtg                                >  2:198834/1‑147 (MQ=255)
                                                                                                                     cAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                <  2:1521471/150‑1 (MQ=255)
                                                                                                                      aaCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCt                               <  1:1048241/150‑1 (MQ=255)
                                                                                                                      aaCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCt                               <  1:1050769/150‑1 (MQ=255)
                                                                                                                      aaCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCt                               >  1:490022/1‑150 (MQ=255)
                                                                                                                       aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                              >  2:60966/1‑150 (MQ=255)
                                                                                                                       aCGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                              <  1:1113415/150‑1 (MQ=255)
                                                                                                                          cAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGTGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAAt                           <  1:1608490/150‑1 (MQ=255)
                                                                                                                          cAGTGTACTGCACCGTTTGCCGATTGTCCTTGAACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAAt                           <  2:1505421/150‑1 (MQ=255)
                                                                                                                            gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTc                         <  2:556573/150‑1 (MQ=255)
                                                                                                                            gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTc                         <  2:389381/150‑1 (MQ=255)
                                                                                                                              gtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                >  1:792777/1‑141 (MQ=255)
                                                                                                                              gtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTc                                <  2:792777/141‑1 (MQ=255)
                                                                                                                               tACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGt                      <  1:373965/150‑1 (MQ=255)
                                                                                                                                aCTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGtt                     <  1:865441/150‑1 (MQ=255)
                                                                                                                                   gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCa                  <  1:1130958/150‑1 (MQ=255)
                                                                                                                                   gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCa                  >  1:157142/1‑150 (MQ=255)
                                                                                                                                   gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCa                  >  1:256074/1‑150 (MQ=255)
                                                                                                                                   gCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCa                  >  1:967364/1‑150 (MQ=255)
                                                                                                                                     aCCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGc                >  2:1373878/1‑150 (MQ=255)
                                                                                                                                         tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCg            <  1:502631/150‑1 (MQ=255)
                                                                                                                                         cccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCg            >  1:880375/4‑150 (MQ=255)
                                                                                                                                         cccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCg            >  1:153853/4‑150 (MQ=255)
                                                                                                                                          ccgCCGATTGTGCAAGGACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTACGTTTCAGCATCgg           >  1:1564037/3‑150 (MQ=255)
                                                                                                                                            gCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCa         <  1:1281030/150‑1 (MQ=255)
                                                                                                                                               gATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCAttt      >  2:257775/1‑150 (MQ=255)
                                                                                                                                                   gTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc  >  1:797598/1‑150 (MQ=255)
                                                                                                                                                   gTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTcc  <  1:526336/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTCCTGGCTGACGTAGCCATTACCTGAAAGCACGCGGGAAACGGTCGCTTTTGAAACCCCGGCGCGCTTCGCCACTTCCAGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCC  >  NC_000913/2839172‑2839468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: