Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,145 G→C 16.4% 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.4% 639.2 / 112.9 239intergenic (‑196/+33)ibsC/serAtoxic membrane protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (91/108);  new base C (16/23);  total (107/132)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-01

Only 100 of 245 total aligned reads displayed.

CATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCA  >  NC_000913/3057046‑3057244
                                                                                                   |                                                                                                   
  ttGATCTGACTTTGATTTATTTTCTGGAGCATACTCGCAAAGTAGAATGCGCAACGCGGCACCGGTGTGGAGAAGGGATAAAAAAATGGGCAAGTCAGTg                                                                                                   <  1:2723459/100‑1 (MQ=255)
    gATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGAGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAc                                                                                                 >  1:3527664/1‑100 (MQ=255)
      tCTGACTTTGATTAATTTTCTGGAACATACTCAAAAAATAGAATGCACAACGAGACAATGGTGTGGAGAATTGATAAAAAAACGGGCAAGTCAGTg                                                                                                   <  1:4175288/96‑1 (MQ=255)
       cTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTg                                                                                              >  1:1615027/1‑100 (MQ=255)
        tGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                             >  2:1345506/1‑100 (MQ=255)
        tGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                             <  2:2344385/100‑1 (MQ=255)
         gACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGcc                                                                                            >  2:752752/1‑100 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCg                                                                                          <  2:356444/100‑1 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                          <  2:4085263/100‑1 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                          >  1:4094004/1‑100 (MQ=255)
           cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGGAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                          >  1:2913313/1‑100 (MQ=255)
            tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                          <  1:778474/99‑1 (MQ=255)
            tttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGt                                                                                         >  2:156849/1‑100 (MQ=255)
               gATTTATTTTTTGGAGCAGACTCGCAAAGTAAAATGCGTAACGCGTCAATGGTGTGGAAAAAAAAAAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                      <  1:1876681/100‑3 (MQ=16)
               gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                      >  2:172039/1‑100 (MQ=255)
                atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                                     <  2:2184743/100‑4 (MQ=38)
                atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                      <  1:3956689/99‑3 (MQ=38)
                atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt                                                                                     >  2:1519130/1‑100 (MQ=255)
                atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAt                                                                                     >  1:1269776/1‑100 (MQ=255)
                 tttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtt                                                                                    >  2:3060843/1‑100 (MQ=255)
                 tttatgtGCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttttt                                                                                    <  2:4141055/100‑5 (MQ=37)
                  ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAttt                                                                                   >  2:3405173/1‑100 (MQ=255)
                  ttatttTCTGGAGCAGACTCGCAAAGTAGAATGAGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttttt                                                                                   <  2:2044280/100‑6 (MQ=37)
                    atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                 <  1:4382071/100‑1 (MQ=255)
                    atttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTc                                                                                 <  2:2837749/100‑1 (MQ=255)
                     ttttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                <  1:2334479/100‑1 (MQ=255)
                      tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                               <  2:3545960/100‑1 (MQ=255)
                       ttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                              >  2:4063626/1‑100 (MQ=255)
                       ntCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGGGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                               >  1:1385792/2‑99 (MQ=255)
                        tCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGTGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCagag                                                                             <  1:3707297/100‑1 (MQ=255)
                           ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTGAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                          >  1:3480551/1‑100 (MQ=255)
                             aGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAggg                                                                        <  2:1276452/100‑1 (MQ=255)
                              gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                       <  1:1884382/100‑1 (MQ=255)
                                 gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                                                                    >  2:3975292/1‑100 (MQ=255)
                                   cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTa                                                                  >  2:224785/1‑100 (MQ=255)
                                   cTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAATTTCAGAGAAGGGGAATTa                                                                  >  1:4388455/1‑100 (MQ=255)
                                    tCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGaaatca                                                                 <  1:3816716/100‑5 (MQ=255)
                                     cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                >  2:2574053/1‑100 (MQ=255)
                                      gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGCGGAATTAg                                                                 >  1:2458531/1‑98 (MQ=255)
                                      gCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTAAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTa                                                               >  2:3879299/1‑100 (MQ=255)
                                       cAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAc                                                              >  2:2787090/1‑100 (MQ=255)
                                        aaaGTAGAATGCGCAACGCGGCAACTGTGTGGAGAAGGGATAAAAAAACGGGCAAGGCAATGACCTGCCCGTTGATTTTCAGAGAAGTGGAATTAGTAca                                                             <  1:4189417/100‑1 (MQ=255)
                                         aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                            <  2:2913330/100‑1 (MQ=255)
                                          aGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcagc                                                           <  2:2051368/100‑1 (MQ=255)
                                                cgatctaaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACgg                                                     <  1:2311427/94‑1 (MQ=255)
                                                   cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                  <  2:373430/100‑1 (MQ=255)
                                                     cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAACCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                                <  2:3090089/100‑1 (MQ=255)
                                                     cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                                <  1:1594730/100‑1 (MQ=255)
                                                     cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                                >  1:4334317/1‑100 (MQ=255)
                                                      aaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGa                                               >  2:1629760/1‑100 (MQ=255)
                                                      aaCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGa                                               <  1:508645/100‑1 (MQ=255)
                                                        cgcgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAAt                                             <  2:1269780/100‑1 (MQ=255)
                                                          cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATgg                                           <  1:3732943/100‑1 (MQ=255)
                                                          cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATgg                                           <  1:223394/100‑1 (MQ=255)
                                                            ggaagggTGTGGAGAAGGGATAAATAAACGGGGCAGTCAGTGTCGTGCCCGTTGGTTTTCAGAGGAGGGGAATTAGTACGGCAGACGAGCGCGAATgg                                           <  1:3565279/93‑1 (MQ=255)
                                                              aaCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAcc                                       <  1:928189/100‑1 (MQ=255)
                                                               aCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAccc                                      <  1:2000739/100‑1 (MQ=255)
                                                                 ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                                    <  2:1780845/100‑1 (MQ=255)
                                                                 ggTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCgg                                    >  1:3441074/1‑100 (MQ=255)
                                                                  gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATggagatcgga                                   >  2:4169966/1‑92 (MQ=255)
                                                                  gtgtGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGa                                   <  1:3041283/100‑1 (MQ=255)
                                                                   tgtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGaa                                  <  2:1489386/100‑1 (MQ=255)
                                                                    gtgGAGAAGGGATAAAAAAGCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGCGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                 <  1:1058724/100‑1 (MQ=255)
                                                                    gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                 <  1:984675/100‑1 (MQ=255)
                                                                    gtgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                 <  2:83853/100‑1 (MQ=255)
                                                                     tgGAGAAGGGAGAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATa                                <  2:812837/100‑1 (MQ=255)
                                                                      ggAGAAGGGATAATAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                               >  2:2961727/1‑100 (MQ=255)
                                                                           aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                          >  2:3190789/1‑100 (MQ=255)
                                                                            gggATAAAAAAATGGGCAAGTGAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGTTACCCGGAATAGCttgc                          <  1:837404/99‑3 (MQ=255)
                                                                            gggATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCa                         >  2:460427/1‑100 (MQ=255)
                                                                            gggATAAAAAAACGGACGAGTCAGTGACCTGCCCGTTGATTTTCAGATAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCtttn                          <  2:565180/99‑2 (MQ=255)
                                                                             aaaaaaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCa                         <  1:3502680/95‑1 (MQ=255)
                                                                              tgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGGGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAt                        >  1:2615214/9‑99 (MQ=25)
                                                                              gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGCATGGTACCCGGACTAGCTTTCAt                        >  1:2174689/1‑99 (MQ=255)
                                                                              gATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                       <  1:1749171/100‑1 (MQ=255)
                                                                               aTAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTg                      <  2:320346/100‑1 (MQ=255)
                                                                                 naaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGACTTAGTACAGCAGAGGGGTGCGAATGGTAGCCGGAATAGCTTTCATTGc                     >  1:3147688/2‑99 (MQ=255)
                                                                                 aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                    <  1:2944009/100‑1 (MQ=255)
                                                                                  aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                   >  2:4174593/1‑100 (MQ=255)
                                                                                  aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                   <  2:4240064/100‑1 (MQ=255)
                                                                                  aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                   >  2:4268818/1‑100 (MQ=255)
                                                                                   atcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                  >  1:1670126/4‑100 (MQ=38)
                                                                                   aaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                  >  1:403533/1‑100 (MQ=255)
                                                                                     caaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                >  1:1010992/2‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  1:2058462/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  2:3661409/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  1:1855478/1‑100 (MQ=39)
                                                                                       aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  2:2664008/1‑100 (MQ=39)
                                                                                       aCGGGCACGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc              >  2:2311545/1‑100 (MQ=255)
                                                                                         gggCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGCATAGCTTTCATTGCCTGCAgcgc            >  1:1998034/1‑100 (MQ=255)
                                                                                         gggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcgc            <  1:3365330/100‑1 (MQ=255)
                                                                                          ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt           <  1:1974405/100‑1 (MQ=255)
                                                                                          ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt           >  2:1119679/1‑100 (MQ=255)
                                                                                            cAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttt         >  1:959748/1‑100 (MQ=255)
                                                                                              aGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt       >  1:2479509/1‑100 (MQ=255)
                                                                                                 cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg    <  2:261512/100‑1 (MQ=255)
                                                                                                  aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGc   >  2:1148563/1‑100 (MQ=255)
                                                                                                  aGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGc   <  1:106610/100‑1 (MQ=255)
                                                                                                   gTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTGGTGGAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg    <  1:720324/98‑1 (MQ=255)
                                                                                                   gTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCa  >  2:908867/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CATTGATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCA  >  NC_000913/3057046‑3057244

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: