Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,145 G→C 16.5% intergenic (‑196/+33) ibsC ← / ← serA toxic membrane protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1450GC16.5% 576.4 / 99.2 218intergenic (‑196/+33)ibsC/serAtoxic membrane protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (95/86);  new base C (19/17);  total (114/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 219 total aligned reads displayed.

CATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGC  >  NC_000913/3057046‑3057243
                                                                                                   |                                                                                                  
    naTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                                >  1:1297265/2‑100 (MQ=255)
    gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                                <  1:3886925/100‑1 (MQ=255)
     aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGGATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAACAAACGGGCAAGACAGTGaac                                                                                               >  2:185104/1‑98 (MQ=255)
     aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                               >  2:1701045/1‑100 (MQ=255)
        tGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                            >  1:2094177/1‑100 (MQ=255)
        tGACTTTGATTTATTTTCTGGAGCAGACTAGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGAct                                                                                               <  2:4082683/97‑2 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGccc                                                                                          <  1:1969978/99‑1 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCg                                                                                         <  2:1158024/100‑1 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                         <  1:1178715/100‑1 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                         >  1:2668134/1‑100 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                         >  1:475021/1‑100 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                         >  1:1587049/1‑100 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                         >  1:4101606/1‑100 (MQ=255)
           cTTTGATTTATTTTCTGAAGCAGCCTCGCACAGTAGAATGCGCAACGCGGCACCGGCGCGGAGAAGGGATAAAAAAACGAGAAAGTAAGTTAcc                                                                                               >  2:3289759/1‑94 (MQ=255)
               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAATGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGCGCAAGTAAGTGACATGCCCGTTg                                                                                      >  2:3391872/1‑99 (MQ=255)
               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                     <  2:2063121/100‑3 (MQ=38)
               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                     <  1:2919876/100‑3 (MQ=38)
               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                     <  1:1451280/100‑3 (MQ=38)
                atttgtttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCAGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                                    <  2:2652597/100‑4 (MQ=37)
                 tttatttaCTGGAGCAGGCTCGCAAAGTAGAATGCGCAACGAGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                                    <  2:1558739/99‑4 (MQ=37)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                                   <  1:1699774/100‑5 (MQ=38)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCGGGTGAt                                                                                    >  2:2097806/1‑99 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                   <  1:3343347/100‑1 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                   >  2:3648395/1‑100 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                   >  1:1317971/1‑100 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                   >  1:552842/1‑100 (MQ=255)
                   tatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                                 >  1:2185658/1‑100 (MQ=255)
                   tatgtTGTGGAGCAGGCTGGCATAGTAGAATGCGCAACGCGGCCACGGTGTGGAGAATGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttt                                                                                      <  2:3319940/95‑2 (MQ=21)
                    atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                >  1:776363/1‑100 (MQ=255)
                       ttCTGGAGCAGACTCGCAAAGTAGAAGGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                             <  2:413967/100‑1 (MQ=255)
                        tCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag                                                                            >  2:1880606/1‑100 (MQ=255)
                          tGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAgaga                                                                           <  1:631605/99‑1 (MQ=255)
                           ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                         >  1:2721629/1‑100 (MQ=255)
                             aGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAggg                                                                       <  1:2988228/100‑1 (MQ=255)
                             aGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAggg                                                                       <  2:2185645/100‑1 (MQ=255)
                              gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                      >  1:371098/1‑100 (MQ=255)
                               cAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGa                                                                     <  1:2569988/100‑1 (MQ=255)
                                aGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaa                                                                    >  2:2959598/1‑100 (MQ=255)
                                  aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                  <  1:1518538/100‑1 (MQ=255)
                                   cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                 >  1:2561424/1‑100 (MQ=255)
                                   cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                 >  1:3625029/1‑100 (MQ=255)
                                   cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                 >  2:2044845/1‑100 (MQ=255)
                                    tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAg                                                                <  1:3320942/100‑1 (MQ=255)
                                    tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAg                                                                <  2:371095/100‑1 (MQ=255)
                                      gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                              >  1:2640264/1‑100 (MQ=255)
                                       cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                             >  1:2826630/1‑100 (MQ=255)
                                        aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTATTAca                                                            >  1:1848175/1‑100 (MQ=255)
                                        aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAca                                                            <  2:320157/100‑1 (MQ=255)
                                         aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                           <  1:3904351/100‑1 (MQ=255)
                                          aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                          <  1:1032566/100‑1 (MQ=255)
                                          aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                          <  2:3057059/100‑1 (MQ=255)
                                          aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                          <  2:3017430/100‑1 (MQ=255)
                                             aGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGa                                                       >  1:1726719/1‑100 (MQ=255)
                                                aTGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACgg                                                    >  1:2702197/1‑100 (MQ=255)
                                                 tGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCTGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACggg                                                   <  2:1413262/100‑1 (MQ=255)
                                                  gcgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                                  >  1:925233/1‑100 (MQ=255)
                                                   cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                 >  2:19894/1‑100 (MQ=255)
                                                     cAACGCGGCAACGGTGTGGAGAAGGGATAAAACAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                               >  2:1292373/1‑100 (MQ=255)
                                                        cgcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt                                            <  1:2379151/100‑1 (MQ=255)
                                                         gcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATg                                           <  1:1604742/100‑1 (MQ=255)
                                                            gCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                        >  2:4021492/1‑100 (MQ=255)
                                                             cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                       >  2:149324/1‑100 (MQ=255)
                                                             cAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                       >  2:2889631/1‑100 (MQ=255)
                                                              aaCGTTGTGGAGAAGGGATAAAAAAACGGGAACGTCGGTGACCTGCCCGGTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAcc                                      <  2:658719/100‑1 (MQ=255)
                                                                 ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCg                                    >  2:223409/1‑99 (MQ=255)
                                                                 ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGCGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                                   <  2:2207851/100‑1 (MQ=255)
                                                                    gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                >  2:24765/1‑100 (MQ=255)
                                                                     tgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAcccag                                   >  1:2946412/1‑94 (MQ=255)
                                                                      ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                              <  2:600911/100‑1 (MQ=255)
                                                                      ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                              >  1:2067769/1‑100 (MQ=255)
                                                                       gagaAGGGATAAAAAAGCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAAGGGAATTAGTACAGCAGACGTGCGCGAATGGTACCTGGAATAGc                             >  2:887249/1‑100 (MQ=255)
                                                                       gagaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGGAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAgn                             <  1:531322/100‑2 (MQ=255)
                                                                        agaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCt                            <  1:1158172/100‑1 (MQ=255)
                                                                         gaAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGGATGGTACCCGGAATAGCtt                           >  2:2596925/1‑100 (MQ=255)
                                                                          aaGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCttt                          <  2:546527/100‑1 (MQ=255)
                                                                              gATAAAAAAACGGGTAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                      >  1:1686857/1‑100 (MQ=255)
                                                                              gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                         <  2:2915552/97‑1 (MQ=255)
                                                                                tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                    <  2:2417638/100‑1 (MQ=255)
                                                                                tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                    <  1:1243773/100‑1 (MQ=255)
                                                                                tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                    >  1:3430834/1‑100 (MQ=255)
                                                                                 aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                   <  1:1006550/100‑1 (MQ=255)
                                                                                  aatcaaCGGGCAGGTCACTGACTTGCCCGTTTATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                  >  1:232392/5‑100 (MQ=37)
                                                                                   atcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                 >  2:483642/4‑100 (MQ=38)
                                                                                    tcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                >  2:3941729/3‑100 (MQ=38)
                                                                                     caaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa               >  1:214139/2‑100 (MQ=38)
                                                                                      aaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAg              >  1:2502680/1‑100 (MQ=39)
                                                                                      aaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTGTCATTGCCTGccg              >  2:3894050/1‑98 (MQ=38)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGTAGTTAGTACAGCAGACGGGCGCGAATGGTACCGGGAATAGCGTTCATTTCCTGCAgc             >  2:343584/1‑100 (MQ=37)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGGACCAGGAATAGCCTTCATTGCCTGCAg              >  2:1303753/1‑99 (MQ=38)
                                                                                        cGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg            <  2:2163714/100‑1 (MQ=255)
                                                                                         gggCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc           >  1:3475002/1‑100 (MQ=255)
                                                                                            naaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGCCGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttt        >  1:3489966/2‑100 (MQ=255)
                                                                                             aaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGCATAGCTTTCATTGCCTGCAGCGCtttt       >  2:1649543/1‑100 (MQ=255)
                                                                                             aaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt       <  1:3150463/100‑1 (MQ=255)
                                                                                             aaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt       >  2:2822692/1‑100 (MQ=255)
                                                                                             aaGTCAGTGACCTGCCCGTTGATTGTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt       >  2:2308417/1‑100 (MQ=255)
                                                                                              aatCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt      >  1:2789858/3‑100 (MQ=255)
                                                                                                tCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGTTGGGCGCGAATGGTGCCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg    <  2:34023/100‑1 (MQ=255)
                                                                                                 cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg   <  1:2836974/100‑1 (MQ=255)
                                                                                                  aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGc  <  1:3990479/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGC  >  NC_000913/3057046‑3057243

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: