Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,083 G→C 15.9% intergenic (‑41/+18) alsB ← / ← alsR D‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,0830GC15.9% 381.7 / 60.5 148intergenic (‑41/+18)alsB/alsRD‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR
Reads supporting (aligned to +/- strand):  ref base G (62/59);  new base C (10/13);  total (74/74)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 155 total aligned reads displayed.

CTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC  >  NC_000913/4311984‑4312182
                                                                                                   |                                                                                                   
cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATg                                                                                                     >  2:482574/1‑100 (MQ=255)
  ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                                    >  1:1005134/1‑99 (MQ=255)
  ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                                   >  1:687281/1‑100 (MQ=255)
    ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATATTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                                 <  2:2230988/100‑1 (MQ=255)
     tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                                 >  2:401442/1‑99 (MQ=255)
     tGACAACATTAAGCCCAAGAGTGTAACAATGAAAGATTTCAGATATTTATTCAAAATAATGTTCTCAAGGTAAGGGGACACGTAACTAAGAGATGGCAGa                                                                                                >  2:1511842/1‑100 (MQ=255)
       acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                              <  2:250090/100‑1 (MQ=255)
        caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                             >  2:2011324/1‑100 (MQ=255)
        caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                             >  2:2136446/1‑100 (MQ=255)
        caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAGCGATGGCGGAAtt                                                                                             >  2:2138019/1‑100 (MQ=255)
        caacaTTAAGCACACGAGTGTGCCGGTGAAATATTTCAGGTATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                               <  1:351196/98‑1 (MQ=255)
          acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                           <  1:1074984/100‑1 (MQ=255)
           caTTATTCCCACTAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                          <  1:1294429/100‑1 (MQ=37)
           caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                          <  2:1786458/100‑1 (MQ=38)
           caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                          <  1:1412274/100‑1 (MQ=38)
           caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                          <  1:1645813/100‑1 (MQ=255)
           caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                          <  2:1511754/100‑1 (MQ=255)
            aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                         <  1:754046/100‑1 (MQ=38)
            aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                         <  2:1524874/100‑1 (MQ=38)
            aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                         <  1:2220007/100‑1 (MQ=255)
            aTTAAGCCCACGAGTGAGCCGATGAAATATTTCAGATATTTATTCATAATAATGTTGTCACGATAAGGGGACAGGAAACTCCACGATGGCGGAATTa                                                                                            >  2:241080/1‑97 (MQ=255)
             ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                        <  1:513527/100‑2 (MQ=38)
               aaccccACGAGTGTGCCACTAAAATATTTCAAATATTTATTCATAATAATGTTCTCACGGTAAGCGGAACGGTAACTCCGCGATAAAACAATTACCttt                                                                                       >  2:1962111/4‑99 (MQ=255)
               aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                                      <  1:1726757/100‑4 (MQ=37)
               aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                                      <  1:2305928/100‑4 (MQ=37)
               aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtttt                                                                                      >  1:818185/1‑100 (MQ=255)
                aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                     >  1:1607487/1‑100 (MQ=255)
                aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                     >  1:1603209/1‑100 (MQ=255)
                 gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGTGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                    <  1:1690549/100‑1 (MQ=255)
                 gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                     >  2:1208607/1‑99 (MQ=255)
                 gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                    >  1:391424/1‑100 (MQ=255)
                 gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                    >  1:2194841/1‑100 (MQ=255)
                  cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcccc                                                                                   <  1:656943/100‑7 (MQ=25)
                   ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                  <  2:1967273/100‑1 (MQ=255)
                    cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAGGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                 >  2:1678443/1‑100 (MQ=255)
                    cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                 <  1:2345856/100‑1 (MQ=255)
                      cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                               <  2:110870/100‑1 (MQ=255)
                      cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                               <  1:457921/100‑1 (MQ=255)
                      cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                               >  2:428403/1‑100 (MQ=255)
                      cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                               <  1:378417/100‑1 (MQ=255)
                        aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                             >  2:580535/1‑100 (MQ=255)
                         gtgtGCAGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGAGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                             <  2:989490/99‑1 (MQ=255)
                          tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                           <  2:973260/100‑1 (MQ=255)
                             gCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                                        >  1:988218/1‑100 (MQ=255)
                              ccGCTGAAATATTTCAGATATTTATTAATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGACTTACCTTTTTCTTATTTCAgcgcg                                                                       >  2:978550/1‑100 (MQ=255)
                               cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                      >  1:1546116/1‑100 (MQ=255)
                                    aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGTTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                                 >  2:81831/1‑100 (MQ=255)
                                     aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                <  1:1427467/100‑1 (MQ=255)
                                     aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                <  2:1244237/100‑1 (MQ=255)
                                     aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                <  1:1342248/100‑1 (MQ=255)
                                      atatTTCAGATATTTATTCATAATAATGTTCTAACGGTAAGGGGACAAGTAACTACGCGATGGCGGAATTACGTTGTTCTTATTTCAGCGCGCCTg                                                                   >  2:2395906/1‑96 (MQ=255)
                                       tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                              <  2:3013/100‑1 (MQ=255)
                                        atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                             >  1:491260/1‑100 (MQ=255)
                                         tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                            <  1:1478820/100‑1 (MQ=255)
                                          ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                           <  2:892406/100‑1 (MQ=255)
                                            cAGATATTTATTCATAATAATGTTCTCAAGGTAAGGGGACAGGTAACTACGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgaagaa                                                         >  2:183456/1‑100 (MQ=255)
                                                tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                                     >  1:1251423/1‑100 (MQ=255)
                                                tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                                     >  2:862648/1‑100 (MQ=255)
                                                  tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                                   <  2:2194815/100‑1 (MQ=255)
                                                  tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                                   <  1:1107622/100‑1 (MQ=255)
                                                    tattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                                 <  1:196768/100‑1 (MQ=255)
                                                     attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                >  1:1114944/1‑100 (MQ=255)
                                                     attCATAANAATGTTCTCACGGTAAGNGGACNGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATNAACAAt                                                >  1:2270362/1‑100 (MQ=255)
                                                       tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGGGGAATTACCTTTTTATTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                              <  1:2179733/100‑1 (MQ=255)
                                                         ataataATGTTCTCACGGTAAGCAGAAATGTAACTCCGCGATGGGGAAATTACCTATGTCTGAAATTAGCGGGCAAGGTGAAAAGAAATCAATAATTgc                                             >  2:1072997/1‑99 (MQ=255)
                                                          taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                           >  2:45858/1‑100 (MQ=255)
                                                          taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                           >  2:73081/1‑100 (MQ=255)
                                                            ataatCTTCACACGGTAAGCGGCTAGGTAACTCCGTGATCGCGGGTTTACGGTTGTCTTATTTCAGCGCGCCTGGTGATAAGGATTCATCAATTCGGCCg                                         >  2:1641036/1‑100 (MQ=255)
                                                             aaatGTTCTCACGGTAAGGGGACAGGAAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                        >  2:1673833/2‑100 (MQ=255)
                                                                    ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                                 <  2:2126015/100‑1 (MQ=255)
                                                                     tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                >  1:1057970/1‑100 (MQ=255)
                                                                      cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                >  2:2069947/1‑99 (MQ=255)
                                                                        cGGTAAGGGGACAGGTAACTCCGCGATAGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGTTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGAata                             >  2:913655/1‑100 (MQ=255)
                                                                          gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                           >  1:2278853/1‑100 (MQ=255)
                                                                            taagaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTACatattaa                         >  1:799331/9‑100 (MQ=25)
                                                                            aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                         <  2:206642/100‑1 (MQ=255)
                                                                             atgggACAGGTAACTACGCGATGGCGGAATTAACTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                         <  1:1848274/97‑1 (MQ=255)
                                                                              ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                       <  2:1160319/100‑1 (MQ=255)
                                                                               gtgACAGGTAACGCCGCGATGGCGGAATTACCTTTTTCTGAGTTCAGCGCGCTTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                       <  2:2127705/97‑1 (MQ=255)
                                                                                ggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                     <  1:959283/100‑1 (MQ=255)
                                                                                  aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg                   >  1:1992869/3‑100 (MQ=38)
                                                                                  aCAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg                   <  2:818176/100‑1 (MQ=255)
                                                                                   aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                  >  1:210396/2‑100 (MQ=38)
                                                                                     ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  1:184293/1‑100 (MQ=38)
                                                                                     ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  2:1822133/1‑100 (MQ=255)
                                                                                     ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                >  1:1013042/1‑100 (MQ=255)
                                                                                      gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg               >  1:295572/1‑100 (MQ=255)
                                                                                      gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg               <  2:865783/100‑1 (MQ=255)
                                                                                        aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGtt             >  1:1626479/1‑100 (MQ=255)
                                                                                         aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc            >  2:651453/1‑100 (MQ=255)
                                                                                         aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc            <  2:1349095/100‑1 (MQ=255)
                                                                                          cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa           <  1:2110902/100‑1 (MQ=255)
                                                                                           tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCaa          <  2:1607467/100‑1 (MQ=255)
                                                                                              gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt       >  1:1390298/1‑100 (MQ=255)
                                                                                              gcgATGGCAGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt       <  1:120087/100‑1 (MQ=255)
                                                                                               cgATGGCCGAATTACCCTTTTCTTATTTCAGCCAGCCCGATGAGCAGAAATCAACAACTGCGCAGGATTTTTGCATATTAATATTACATTGATCAATg        >  2:1880562/1‑98 (MQ=255)
                                                                                                 aTGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGa    >  1:1740567/1‑100 (MQ=255)
                                                                                                 aTGGCGGAATTACATTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGa    <  1:58476/100‑1 (MQ=255)
                                                                                                   ggCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTCATATTAGCTTGTTCAATGTTGAg   >  2:352970/1‑99 (MQ=255)
                                                                                                   ggCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc  >  2:1418870/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGC  >  NC_000913/4311984‑4312182

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: