Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,074 C→T 17.6% 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,0740CT17.6% 393.2 / 66.7 165intergenic (‑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 (68/68);  new base T (8/21);  total (76/89)
Fisher's exact test for biased strand distribution p-value = 3.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 169 total aligned reads displayed.

GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173
                                                                                                   |                                                                                                    
gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAAc                                                                                                      <  2:1605658/100‑1 (MQ=255)
 caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACt                                                                                                     >  1:1494748/1‑100 (MQ=255)
    aaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCcg                                                                                                  <  2:282129/100‑1 (MQ=255)
     aGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCcgc                                                                                                 <  1:2977532/100‑1 (MQ=255)
       cgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGa                                                                                               >  1:394085/1‑100 (MQ=255)
          tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                            >  2:1105150/1‑100 (MQ=255)
           ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                           >  1:2645003/1‑100 (MQ=255)
           ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                           >  1:351155/1‑100 (MQ=255)
            gTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCg                                                                                          <  1:780255/100‑1 (MQ=255)
             ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATcgcgn                                                                                         <  1:1388462/100‑2 (MQ=255)
              tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                         <  2:1409702/99‑1 (MQ=255)
              tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGa                                                                                        >  2:2726477/1‑100 (MQ=255)
                  aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                    <  1:2412668/100‑1 (MQ=255)
                  aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTa                                                                                    <  1:2662697/100‑1 (MQ=255)
                   acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                   >  1:1063650/1‑100 (MQ=255)
                   acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                   >  1:1477262/1‑100 (MQ=255)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                  <  1:3183006/100‑1 (MQ=38)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAcc                                                                                  <  2:683282/100‑1 (MQ=38)
                    caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                  <  1:2333806/100‑1 (MQ=255)
                     aTTAAGCCCAGGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                 <  1:3585513/100‑1 (MQ=38)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                 <  2:1742333/100‑1 (MQ=38)
                     aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                 >  1:1425184/1‑100 (MQ=255)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                <  2:3612646/100‑2 (MQ=38)
                      ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                <  1:2806504/100‑2 (MQ=38)
                      ttAAGCCCACGAGTGTGCCGCTGAAATACTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                <  2:159595/100‑2 (MQ=38)
                       tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttt                                                                               <  2:1548537/100‑1 (MQ=255)
                        aaGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtc                                                                              <  1:816292/100‑4 (MQ=37)
                         cgatctcGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttc                                                                             <  2:1784535/94‑2 (MQ=255)
                         aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcc                                                                             <  1:2260602/100‑5 (MQ=37)
                         aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTACTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                             >  1:2148295/1‑100 (MQ=255)
                          gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtccc                                                                            <  2:1158403/100‑6 (MQ=37)
                          gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           <  2:2977616/100‑1 (MQ=255)
                          gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                           >  2:3001092/1‑100 (MQ=255)
                           cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGTGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                          <  1:1255104/100‑1 (MQ=255)
                           cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcccc                                                                           <  2:426778/100‑7 (MQ=37)
                           cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                          >  2:2910203/1‑100 (MQ=255)
                            ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTcn                                                                          <  2:1306252/99‑2 (MQ=255)
                             cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                        <  2:215281/100‑1 (MQ=255)
                               ctaGTGGGCAGCTTAAATATTTCAGATATTTATTCATACTACTGTTCTCACGGTAAGGTGACAGTTAACTCCGCGATGGCGGAATTCCCTTTTCCTTAtt                                                                      <  2:2621766/98‑1 (MQ=255)
                               cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttcagatcg                                                                       >  1:1784543/1‑93 (MQ=255)
                                 aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                    <  1:3099139/100‑1 (MQ=255)
                                    gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                 <  1:212777/100‑1 (MQ=255)
                                    gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                 >  2:1448866/1‑100 (MQ=255)
                                      gCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgc                                                               >  1:1728521/1‑100 (MQ=255)
                                       ccGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                              >  1:259979/1‑100 (MQ=255)
                                         gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                            <  1:3334996/100‑1 (MQ=255)
                                         gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                            >  2:2201155/1‑100 (MQ=255)
                                            gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                         >  1:23975/1‑100 (MQ=255)
                                               atatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                      <  2:2353514/100‑1 (MQ=255)
                                                 atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTGAGCGCGCCTGGTgaga                                                    >  2:2123781/1‑100 (MQ=255)
                                                 atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                    >  1:2637086/1‑100 (MQ=255)
                                                     cAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgaagaa                                                >  1:1081604/1‑100 (MQ=255)
                                                       gATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                              <  2:2690156/100‑1 (MQ=255)
                                                        atatTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                             <  2:1714146/100‑1 (MQ=255)
                                                         tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                            <  1:1116346/100‑1 (MQ=255)
                                                          atttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                           >  1:1054003/1‑100 (MQ=255)
                                                           tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATCACCTTTTTCTTATTTCAGCGCGCCTGGTCAGAAGAAATcacc                                          <  1:2468186/100‑3 (MQ=255)
                                                            ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                         <  1:3574788/100‑1 (MQ=255)
                                                             tattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                        >  2:460587/1‑100 (MQ=255)
                                                               ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                      >  2:3486503/1‑100 (MQ=255)
                                                                tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTg                                     >  1:913117/1‑100 (MQ=255)
                                                                 cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                    <  1:1168569/100‑1 (MQ=255)
                                                                   taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                  >  2:3508321/1‑100 (MQ=255)
                                                                   taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                  <  2:741051/100‑1 (MQ=255)
                                                                   taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                  <  1:1328372/100‑1 (MQ=255)
                                                                     ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                >  1:2516698/1‑100 (MQ=255)
                                                                     ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                >  1:1000014/1‑100 (MQ=255)
                                                                     ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                >  1:2585356/1‑100 (MQ=255)
                                                                     ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                >  2:1015000/1‑100 (MQ=255)
                                                                      taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               >  1:890150/1‑100 (MQ=255)
                                                                      taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                               <  2:1843512/100‑1 (MQ=255)
                                                                        atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtt                             >  1:2537069/1‑100 (MQ=255)
                                                                          gTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttt                           >  1:2717741/1‑100 (MQ=255)
                                                                          gTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttt                           <  1:3099368/100‑1 (MQ=255)
                                                                             ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  2:3483004/100‑1 (MQ=255)
                                                                             ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                        <  2:3442074/100‑1 (MQ=255)
                                                                              tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                       >  1:455659/1‑100 (MQ=255)
                                                                                     aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                <  2:2110404/100‑1 (MQ=255)
                                                                                      aagaaaaatgTAATTCCGCCATCGCGGAGATACCTTTTTCTTATTTCAGCGCACCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               >  2:1236224/10‑100 (MQ=25)
                                                                                      aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat               >  1:3385239/1‑100 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  2:351446/1‑100 (MQ=255)
                                                                                       ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata              >  2:2491391/1‑100 (MQ=255)
                                                                                            cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         <  1:2251732/100‑1 (MQ=255)
                                                                                            cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         <  1:786426/100‑1 (MQ=255)
                                                                                            cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         <  2:2843678/100‑1 (MQ=255)
                                                                                            aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc         >  2:50346/2‑100 (MQ=38)
                                                                                             aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  1:408331/1‑100 (MQ=38)
                                                                                             aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  2:2349166/1‑100 (MQ=38)
                                                                                             aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        <  1:2195446/100‑1 (MQ=255)
                                                                                             aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        <  1:2882183/100‑1 (MQ=255)
                                                                                             aGGGAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        <  2:3525973/100‑1 (MQ=255)
                                                                                              ggTAATTCCGCCATCGCGGAGTTACCTGTCCCCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt        >  1:3186369/1‑100 (MQ=21)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  1:1523628/1‑100 (MQ=255)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  2:803956/1‑100 (MQ=255)
                                                                                              ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       >  2:3214206/1‑100 (MQ=255)
                                                                                               gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCTCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt       <  2:614078/99‑1 (MQ=255)
                                                                                               gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg      <  1:2998313/100‑1 (MQ=255)
                                                                                                  aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   <  1:3209224/100‑1 (MQ=255)
                                                                                                  aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTc   >  1:1277367/1‑100 (MQ=255)
                                                                                                   cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa  >  1:3025092/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCA  >  NC_000913/4311975‑4312173

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