Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,080 G→C 15.9% intergenic (‑38/+21) 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,0800GC15.9% 555.4 / 94.1 215intergenic (‑38/+21)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 (89/90);  new base C (17/17);  total (107/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

Only 100 of 232 total aligned reads displayed.

AGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTG  >  NC_000913/4311980‑4312179
                                                                                                    |                                                                                                   
aGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCGCTGAAATATTTCAGATAGTTAGTCATAATAATGTTCTCTCGGTAAGGGGACAGGTAACTCcgcg                                                                                                     <  1:181968/100‑1 (MQ=255)
 gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCcgcg                                                                                                     <  2:1865370/100‑1 (MQ=255)
 gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCcgcg                                                                                                     >  1:77290/1‑100 (MQ=255)
  cgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGa                                                                                                    <  1:6833/100‑1 (MQ=255)
  cgcTGGTTGACAACATTAAGCCAAAGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGa                                                                                                    >  1:453823/1‑100 (MQ=255)
   gcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGAt                                                                                                   <  1:1060832/100‑1 (MQ=255)
   gcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGAt                                                                                                   >  1:3389334/1‑100 (MQ=255)
     tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                                 >  2:417037/1‑100 (MQ=255)
     tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATgg                                                                                                 <  2:2481956/100‑1 (MQ=255)
      ngTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGc                                                                                                >  1:2273990/2‑100 (MQ=255)
        ttGACACCACTACGCCCACCGGTGTGCCGCTGAAACATTTCATATATTCATTCATAACACTGTTCTCGCGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                              <  1:4171224/100‑1 (MQ=255)
        ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCgg                                                                                              <  2:3736034/100‑1 (MQ=255)
        ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                              <  1:1074515/100‑1 (MQ=255)
        ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCgg                                                                                              <  2:152667/100‑1 (MQ=255)
            caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                          <  2:4342848/100‑1 (MQ=255)
            caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                          >  1:2238003/1‑100 (MQ=255)
              acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                        >  2:666703/1‑100 (MQ=255)
                aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                      <  2:1048428/100‑1 (MQ=255)
                aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                      <  2:4282058/100‑1 (MQ=255)
                 ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                     <  1:607590/100‑2 (MQ=38)
                 ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                     <  2:529593/100‑2 (MQ=38)
                 ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtt                                                                                     >  1:4353210/1‑100 (MQ=255)
                  tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtg                                                                                     <  2:4387596/99‑2 (MQ=38)
                   aaGCCAACGAGTGTGCAGCTCAACTATTTCAGATATTTATTCATAATAATGTTCACACAGTACGGGGACAGGTACCTCCGCTATGGTGGAATTACCtttt                                                                                   <  1:1461780/100‑1 (MQ=255)
                    aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCtgtcc                                                                                  <  1:3630680/100‑5 (MQ=37)
                     gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                 <  2:3536206/100‑1 (MQ=255)
                      cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCCCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACATTTTTCt                                                                                >  1:352043/1‑100 (MQ=255)
                       ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                               >  1:3157588/1‑100 (MQ=255)
                        cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                              >  2:803636/1‑100 (MQ=255)
                         aCGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAt                                                                             >  2:2941100/1‑100 (MQ=255)
                          cGAGTGTGTCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGTACAGGTAACTCCGCGATTGCGGAATTACCTTTTTCTTAtt                                                                            <  1:4391078/100‑1 (MQ=255)
                          cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                            >  1:975056/1‑100 (MQ=255)
                            aGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTc                                                                          <  1:1734129/100‑1 (MQ=255)
                             gtgtGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCa                                                                         <  1:2021436/100‑1 (MQ=255)
                               gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                       >  2:1408121/1‑100 (MQ=255)
                                   cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                   >  2:4270588/1‑100 (MQ=255)
                                    gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGTAATTAACTTTTTCTTATTTCAGCGCGcc                                                                  >  2:131380/1‑100 (MQ=255)
                                    gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                   <  1:2792622/99‑1 (MQ=255)
                                    gCTGAAATATTTCAGATATTTATTCATAATAATGTTCGCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgcn                                                                  <  1:1423264/100‑2 (MQ=255)
                                     cTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                 >  1:407738/1‑100 (MQ=255)
                                      tGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTg                                                                <  2:1391120/100‑1 (MQ=255)
                                      tGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTg                                                                <  2:2810049/100‑1 (MQ=255)
                                       gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                               <  2:2387371/100‑1 (MQ=255)
                                       gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTgg                                                               >  2:707695/1‑100 (MQ=255)
                                         aaTCTTTCAGATATTTATTCATAATTATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGCATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                             <  2:1298025/100‑1 (MQ=255)
                                           tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                           <  2:790270/100‑1 (MQ=255)
                                            ntttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAACGCGCCTGGTgaga                                                          >  1:1861900/2‑100 (MQ=255)
                                            atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCTGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                          >  2:2631783/1‑100 (MQ=255)
                                            atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                          >  2:1253525/1‑100 (MQ=255)
                                             tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                         <  1:847048/100‑1 (MQ=255)
                                             tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                         <  1:1140432/100‑1 (MQ=255)
                                                cAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgaagaa                                                      >  1:1643186/1‑100 (MQ=255)
                                                  gATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                    >  2:2911236/1‑100 (MQ=255)
                                                    tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                                  >  2:960975/1‑100 (MQ=255)
                                                     atttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                 <  1:1521764/100‑1 (MQ=255)
                                                       ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                               <  2:1192996/100‑1 (MQ=255)
                                                       ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                               <  2:2106934/100‑1 (MQ=255)
                                                         attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                             >  2:3384356/1‑100 (MQ=255)
                                                           tCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTTTTTTTTTTTCGCGCCTGGTGAGAACAAATCAACAAt                                             >  1:4368139/1‑98 (MQ=255)
                                                             ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                        <  1:2536802/100‑1 (MQ=255)
                                                              taataaTGTTCTCACGGTAAGGGGACAGGCAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGGGAAGAAATCAACAATTgcgc                                        >  1:563515/1‑100 (MQ=255)
                                                               aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                       <  2:2274000/100‑1 (MQ=255)
                                                                ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                      <  2:1662360/100‑1 (MQ=255)
                                                                ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                      >  2:2853035/1‑100 (MQ=255)
                                                                  aatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGt                                    >  2:3294257/1‑100 (MQ=255)
                                                                       tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                               >  2:3439012/1‑100 (MQ=255)
                                                                       tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                               >  1:1243095/1‑100 (MQ=255)
                                                                        tctaCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                              <  1:824573/97‑1 (MQ=255)
                                                                        ctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                              >  2:2203240/1‑100 (MQ=255)
                                                                          cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                            <  2:407739/100‑1 (MQ=255)
                                                                            cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCGGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCata                          <  1:4261618/100‑1 (MQ=255)
                                                                            cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCata                          >  2:2744978/1‑100 (MQ=255)
                                                                             ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTgagatc                         >  1:3045441/1‑95 (MQ=255)
                                                                              gTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatatt                        >  1:3322158/1‑100 (MQ=255)
                                                                                taagtataaGGTAATTCCGCCTTCGCGTAGTCACTTTTTTCTTATTTCAGCGCGACTGGTGAGCAGAAATCTACAATTGCGTCGGTTTTTTTCatat                         >  1:2952869/9‑97 (MQ=255)
                                                                                 aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat                     <  2:1776762/100‑1 (MQ=255)
                                                                                 aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCAGTTTTTTGCagattaa                      <  1:1324527/99‑1 (MQ=255)
                                                                                 aGGGGACAGGTAACTCCGCGATGGCGGAATTACATTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat                     >  1:3843499/1‑100 (MQ=255)
                                                                                  ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaata                    <  1:2580612/100‑1 (MQ=255)
                                                                                   gggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat                   >  1:1594241/1‑100 (MQ=255)
                                                                                     gACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatta                 <  2:174220/100‑1 (MQ=255)
                                                                                     gACAGGTAACTCCGCGATGGCGGAATTACCTTTTCCTTACTTCAGCGCGCCTGGTGAGAAGAAAGCAAAAAGTGCGCCGGTCTTTTGCataataatatta                 >  1:1696022/1‑100 (MQ=255)
                                                                                      aaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg                >  1:1000544/3‑100 (MQ=38)
                                                                                       cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc               >  1:3632311/1‑100 (MQ=255)
                                                                                       aaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc               >  1:1611678/2‑100 (MQ=38)
                                                                                        aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTTCATATTAATATTAg                >  1:483373/1‑98 (MQ=38)
                                                                                        aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              >  2:1608299/1‑100 (MQ=38)
                                                                                        aGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGAGCCGGTTTTTTGCATATTAATATTAGCt              >  2:1546467/1‑100 (MQ=38)
                                                                                        aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCTCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              <  1:399112/100‑1 (MQ=255)
                                                                                        aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt              <  2:3037993/100‑1 (MQ=255)
                                                                                        aGGTAACTCCGCGATGGCGGAATTACCTTGTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACGATTGCGCCGGTTTTTTGCATATTAATATTAGCt              <  1:2640075/100‑1 (MQ=255)
                                                                                         ggTAATTTCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAATTCAACAATTGCGCAGGTTTTTTTCATATTAATCTTAACtt             >  1:3696822/1‑100 (MQ=255)
                                                                                         ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt             >  1:139554/1‑100 (MQ=38)
                                                                                         ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt             >  2:3628607/1‑100 (MQ=255)
                                                                                              cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa        >  1:1281769/1‑100 (MQ=255)
                                                                                               tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCaa       <  1:231479/100‑1 (MQ=255)
                                                                                                ccGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt      >  2:2469371/1‑100 (MQ=255)
                                                                                                 cgcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATg     <  1:525100/100‑1 (MQ=255)
                                                                                                  gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt      <  1:515103/98‑1 (MQ=255)
                                                                                                  gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt    <  1:2881400/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
AGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTG  >  NC_000913/4311980‑4312179

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