Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,074 C→T 15.3% intergenic (‑32/+27) 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,0740CT15.3% 346.6 / 47.4 144intergenic (‑32/+27)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 C (64/58);  new base T (9/13);  total (73/71)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-02

Only 100 of 150 total aligned reads displayed.

GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173
                                                                                                   |                                                                                                   
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     >  1:2382349/1‑100 (MQ=255)
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                     <  1:802764/100‑1 (MQ=255)
       cgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGa                                                                                              >  2:296389/1‑100 (MQ=255)
         cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATg                                                                                            >  2:482574/1‑100 (MQ=255)
         cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATg                                                                                            >  2:1097062/1‑100 (MQ=255)
          tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                           <  1:27540/100‑1 (MQ=255)
           ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                           >  1:1005134/1‑99 (MQ=255)
           ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                          <  2:1970138/100‑1 (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:2136446/1‑100 (MQ=255)
                 caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                    >  2:2011324/1‑100 (MQ=255)
                 caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                    >  1:1701425/1‑100 (MQ=255)
                 caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCAGCGATGGCGGAAtt                                                                                    >  2:2138019/1‑100 (MQ=255)
                   acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                  <  2:1827482/100‑1 (MQ=255)
                   acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                  <  1:1074984/100‑1 (MQ=255)
                    caTTATTCCCACTAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                 <  1:1294429/100‑1 (MQ=37)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                 <  2:1786458/100‑1 (MQ=38)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                 <  1:1645813/100‑1 (MQ=255)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                <  2:1524874/100‑1 (MQ=38)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                <  2:965168/100‑1 (MQ=38)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                <  1:754046/100‑1 (MQ=38)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                <  1:432147/100‑1 (MQ=38)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                <  1:2220007/100‑1 (MQ=255)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                               <  1:513527/100‑2 (MQ=38)
                        aaccccACGAGTGTGCCACTAAAATATTTCAAATATTTATTCATAATAATGTTCTCACGGTAAGCGGAACGGTAACTCCGCGATAAAACAATTACCttt                                                                              >  2:1962111/4‑99 (MQ=255)
                        aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                             <  1:1726757/100‑4 (MQ=37)
                        aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtttt                                                                             >  1:818185/1‑100 (MQ=255)
                         aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                            >  1:1607487/1‑100 (MQ=255)
                          gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           >  1:2194841/1‑100 (MQ=255)
                           cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcccc                                                                          <  1:656943/100‑7 (MQ=25)
                             cACGAGTTTGCTGCTGAAATATTTCAAATATTTATTCATAATAATGTTCTCACGGTAAGGGTACAGGTAACTCCGCGATTTCGGAATTACCTTTTTCTTa                                                                        <  2:1884737/100‑1 (MQ=255)
                             cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                        <  1:2345856/100‑1 (MQ=255)
                               cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      >  2:428403/1‑100 (MQ=255)
                               cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      <  1:378417/100‑1 (MQ=255)
                               cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                      <  1:457921/100‑1 (MQ=255)
                                  gtgtGCAGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGAGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                    <  2:989490/99‑1 (MQ=255)
                                   tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                  <  2:973260/100‑1 (MQ=255)
                                    gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGCATTACCTTTTTCTTATTTCAgc                                                                 >  2:853315/1‑100 (MQ=255)
                                      gCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                               >  1:988218/1‑100 (MQ=255)
                                        cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                             >  1:1546116/1‑100 (MQ=255)
                                         gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                              >  2:372762/1‑98 (MQ=255)
                                         gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                            <  1:356070/100‑1 (MQ=255)
                                          cTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                           >  1:1555443/1‑100 (MQ=255)
                                            gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                         >  1:1614579/1‑100 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:1734675/100‑1 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  2:1244237/100‑1 (MQ=255)
                                              aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                       <  1:1427467/100‑1 (MQ=255)
                                               atatTTCAGATATTTATTCATAATAATGTTCTAACGGTAAGGGGACAAGTAACTACGCGATGGCGGAATTACGTTGTTCTTATTTCAGCGCGCCTg                                                          >  2:2395906/1‑96 (MQ=255)
                                                 atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                    >  1:491260/1‑100 (MQ=255)
                                                  tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                   <  1:1364536/100‑1 (MQ=255)
                                                   ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                  <  2:988207/100‑1 (MQ=255)
                                                   ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                  <  2:892406/100‑1 (MQ=255)
                                                    tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaaga                                                 >  1:1451095/1‑100 (MQ=255)
                                                       gATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                              <  1:2233523/100‑1 (MQ=255)
                                                         tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            >  2:862648/1‑100 (MQ=255)
                                                           tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                          <  2:2194815/100‑1 (MQ=255)
                                                           tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                          <  1:1107622/100‑1 (MQ=255)
                                                           tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCAGAATTACCTTTTGCTTATTTCAGCGCGCCTGTTGAGAAGAAATcaac                                          >  2:681785/1‑100 (MQ=255)
                                                             tattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                        <  1:196768/100‑1 (MQ=255)
                                                              attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                       >  1:1114944/1‑100 (MQ=255)
                                                                tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGGGGAATTACCTTTTTATTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                     <  1:2179733/100‑1 (MQ=255)
                                                                 cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                    >  2:889674/1‑100 (MQ=255)
                                                                  ataataATGTTCTCACGGTAAGCAGAAATGTAACTCCGCGATGGGGAAATTACCTATGTCTGAAATTAGCGGGCAAGGTGAAAAGAAATCAATAATTgc                                    >  2:1072997/1‑99 (MQ=255)
                                                                      aaatGTTCTCACGGTAAGGGGACAGGAAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               >  2:1673833/2‑100 (MQ=255)
                                                                            tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                         >  1:2287476/1‑100 (MQ=255)
                                                                            tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                         <  1:662319/100‑1 (MQ=255)
                                                                             ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  2:2126015/100‑1 (MQ=255)
                                                                              tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                       >  1:1057970/1‑100 (MQ=255)
                                                                              tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                       <  2:2399249/100‑1 (MQ=255)
                                                                               cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                       >  2:2069947/1‑99 (MQ=255)
                                                                               cACGGTAAGGGGACAGGTAACTACGCGATGGCGGAATGACATTTTTCTTATATCAGCGCCCCTGGTGAGAAGAAATCAACAATTACGCCGGtttttt                         >  2:256469/1‑97 (MQ=255)
                                                                                   gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                  >  1:2278853/1‑100 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  2:206642/100‑1 (MQ=255)
                                                                                      atgggACAGGTAACTACGCGATGGCGGAATTAACTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  1:1848274/97‑1 (MQ=255)
                                                                                      aagaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               >  2:523957/8‑100 (MQ=37)
                                                                                      aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               <  1:1265694/100‑1 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              <  2:1160319/100‑1 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGAAGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  2:1460449/1‑100 (MQ=255)
                                                                                       agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  1:786035/7‑100 (MQ=25)
                                                                                        gtgACAGGTAACGCCGCGATGGCGGAATTACCTTTTTCTGAGTTCAGCGCGCTTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              <  2:2127705/97‑1 (MQ=255)
                                                                                        gaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat             >  2:746040/6‑100 (MQ=37)
                                                                                         ggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt            <  1:959283/100‑1 (MQ=255)
                                                                                          gACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatta           <  2:1662094/100‑1 (MQ=255)
                                                                                           aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg          >  1:1992869/3‑100 (MQ=38)
                                                                                           aCAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg          <  2:818176/100‑1 (MQ=255)
                                                                                            cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         <  2:1190006/100‑1 (MQ=255)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  1:210396/2‑100 (MQ=38)
                                                                                            aaGGTAATTCCGCCATCGCGGAGATACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:126580/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)
                                                                                                  aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   >  2:651453/1‑100 (MQ=255)
                                                                                                  aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   <  2:1349095/100‑1 (MQ=255)
                                                                                                   cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa  >  1:1681346/1‑100 (MQ=255)
                                                                                                   cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa  >  1:1398384/1‑100 (MQ=255)
                                                                                                   cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa  <  1:2110902/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173

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: