Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,559 A→G 19.4% intergenic (‑382/+15) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,5590AG19.4% 203.8 / 47.8 103intergenic (‑382/+15)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (38/45);  new base G (9/11);  total (47/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

Only 100 of 104 total aligned reads displayed.

AAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAA  >  NC_000913/699460‑699655
                                                                                                   |                                                                                                
aaaGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCa                                                                                                  <  2:1692431/100‑1 (MQ=255)
  aGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGa                                                                                                <  1:1948557/100‑1 (MQ=255)
      tGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTg                                                                                            >  1:1259798/1‑100 (MQ=255)
      tGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTg                                                                                            >  1:425123/1‑100 (MQ=255)
        aaaaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGt                                                                                          <  2:761309/100‑1 (MQ=255)
         aaaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGtt                                                                                         <  2:1042667/100‑1 (MQ=255)
          aaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTg                                                                                        <  2:371841/100‑1 (MQ=255)
          aaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTg                                                                                        <  2:112460/100‑1 (MQ=255)
          aaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTg                                                                                        <  2:279253/100‑1 (MQ=255)
          aaTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTg                                                                                        <  2:2487174/100‑1 (MQ=255)
           aTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTgg                                                                                       <  2:1836730/100‑2 (MQ=255)
           aTTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTgg                                                                                       <  2:1983158/100‑2 (MQ=255)
            ttCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTggc                                                                                      <  2:1619716/100‑3 (MQ=255)
              cAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTAt                                                                                    >  1:665750/1‑100 (MQ=255)
              cAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTAt                                                                                    <  2:1246075/100‑1 (MQ=255)
               aaTATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTAtt                                                                                   >  2:2419924/1‑100 (MQ=255)
                 tatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa                                                                                 <  2:53094/100‑1 (MQ=255)
                  atagTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa                                                                                 <  1:454876/99‑1 (MQ=255)
                   tatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGTTggctctgatg                                                                               <  2:1147668/100‑10 (MQ=25)
                      tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATaa                                                                            <  1:2357655/100‑1 (MQ=255)
                      tGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                             >  1:1065684/1‑99 (MQ=255)
                          ggAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTGTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGGGGTTGCTATTCAGATAACGTc                                                                        >  1:763162/1‑100 (MQ=255)
                           gAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                       <  1:132405/100‑1 (MQ=255)
                           gAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                       <  1:2186084/100‑1 (MQ=255)
                             gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGAt                                                                     >  2:1819480/1‑100 (MQ=255)
                               tGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGAttt                                                                   >  1:181719/1‑100 (MQ=255)
                                gCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                  <  2:2446366/100‑1 (MQ=255)
                                  gTAGGCCTGATGAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                <  2:1134184/100‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                <  1:2568622/100‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                <  2:1013871/100‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                <  1:639269/100‑1 (MQ=255)
                                    aGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCg                                                              <  2:2044587/100‑1 (MQ=255)
                                     ggCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGa                                                             >  1:154480/1‑100 (MQ=255)
                                     ggCCTGATAAGCGTAGCGCATCAGGCAGTTGGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCg                                                              <  1:1531249/99‑1 (MQ=255)
                                      gCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGAc                                                            >  2:2319586/1‑100 (MQ=255)
                                       ccTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACt                                                           <  2:1558323/100‑1 (MQ=255)
                                         tGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGa                                                         <  1:927323/100‑1 (MQ=255)
                                          gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTTCTATTCAGATAACGTCGATTTCAGCGACTGAc                                                        <  2:2463305/100‑1 (MQ=255)
                                          gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGAc                                                        >  2:37319/1‑100 (MQ=255)
                                            tAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACgg                                                      <  1:1105219/100‑1 (MQ=255)
                                            tAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACgg                                                      >  1:1115199/1‑100 (MQ=255)
                                             tctcGTAGCGCATCAGGCAGTTTGTCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACggg                                                     <  2:477815/97‑1 (MQ=255)
                                             aaGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCTGCGACTGACggg                                                     <  2:2082889/100‑1 (MQ=255)
                                             aaGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACggg                                                     >  2:2056510/1‑100 (MQ=255)
                                              aGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGt                                                    >  1:869819/1‑100 (MQ=255)
                                               gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGACCGCAGACGTGGTTGATATGCAGATAAAGTCGCTTTCAGCGACCAACGGGTa                                                   >  1:2367971/1‑100 (MQ=37)
                                                cGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACgggaga                                                  >  1:477813/1‑97 (MQ=255)
                                                cGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTaa                                                  >  2:1022376/1‑100 (MQ=255)
                                                 gTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTaaa                                                 >  1:164735/1‑100 (MQ=255)
                                                  tAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAAt                                                <  2:1909573/100‑1 (MQ=255)
                                                   aGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATc                                               >  2:781450/1‑100 (MQ=255)
                                                     cgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCa                                             >  2:2422294/1‑100 (MQ=255)
                                                     cgcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGCCGATTTCAGCGACTGACGGGTAAATCCa                                             >  1:1589910/1‑100 (MQ=255)
                                                      gcATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAg                                            >  2:451544/1‑100 (MQ=255)
                                                         tCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTg                                         <  2:1968029/100‑1 (MQ=255)
                                                           aGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTggg                                       <  2:1726214/100‑1 (MQ=255)
                                                            ggCAGTTTGGCGTTTGTCATCAGAGCCATCCACGTCCGCATACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTgggg                                      <  1:530586/100‑1 (MQ=255)
                                                            ggCAGTTTGGCGTTTGTCATCAGAGCCAACCAGGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTgggg                                      >  2:2459380/1‑100 (MQ=255)
                                                            ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTgggg                                      >  1:786571/1‑100 (MQ=255)
                                                            ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTgggg                                      >  1:595392/1‑100 (MQ=255)
                                                                gTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGa                                  >  2:216130/1‑100 (MQ=255)
                                                                 tttGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGaa                                 >  2:707923/1‑100 (MQ=255)
                                                                     gCGTTTGTCATCAGAGCCAACCACGTTCGCAGACGTGGGTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGc                             <  1:2589021/100‑1 (MQ=255)
                                                                      cGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                            <  1:216128/100‑1 (MQ=255)
                                                                      cGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTCCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                            <  1:1819470/100‑1 (MQ=255)
                                                                            tcatcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGt                      >  2:1624610/1‑100 (MQ=255)
                                                                             catcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTc                     <  2:1683857/100‑1 (MQ=255)
                                                                             catcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTc                     >  2:372041/1‑100 (MQ=255)
                                                                             cagcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTc                     <  1:1148085/97‑1 (MQ=255)
                                                                              atcaGAGCCAACCACGTCCGCAGACGTGGTTGGTATTCAGATAAGGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTc                     <  1:2070159/99‑1 (MQ=255)
                                                                              atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCg                    <  1:1318682/100‑1 (MQ=255)
                                                                              atcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCg                    <  1:1940363/100‑1 (MQ=255)
                                                                               tgaatagcAACCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGa                   >  2:1058261/8‑100 (MQ=38)
                                                                               tgaatagcAACCACGTCTGCGGACGTGGTTGCGATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGa                   >  2:2458893/8‑100 (MQ=38)
                                                                               tcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGa                   <  1:2111537/100‑1 (MQ=255)
                                                                               tcaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGa                   <  2:1927877/100‑1 (MQ=255)
                                                                               ncaGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGa                   >  1:2120585/2‑100 (MQ=255)
                                                                                caGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGat                  <  2:1231787/100‑1 (MQ=255)
                                                                                  gagCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                <  2:665753/100‑1 (MQ=255)
                                                                                  atagcAACCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                >  1:1270543/5‑100 (MQ=39)
                                                                                  atagcAACCACGTCTGCGGACGAGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATACAGCTGGGGCGGAAAGGCAGACTGTCGata                 >  1:1916131/5‑99 (MQ=38)
                                                                                   agCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtc               >  2:1334769/1‑100 (MQ=255)
                                                                                    gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATATCGTCGATTTCAGCGACTGACGGGTAAAACCAGCTCGGGCGGAAAGGCATACTGTCGATAtct              >  1:1365605/1‑99 (MQ=255)
                                                                                    gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg              >  2:1198699/1‑100 (MQ=255)
                                                                                    gCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcg              >  2:393751/1‑100 (MQ=255)
                                                                                     gcAACCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCCGGGGCGGAAAGGCCTACTGTCGATAtcgt             >  1:785616/2‑100 (MQ=255)
                                                                                     gcAACCACGTCTACGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgt             >  2:1701001/2‑100 (MQ=255)
                                                                                     ccAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgt             <  1:463930/100‑1 (MQ=255)
                                                                                      cAACCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtc            >  2:2455830/1‑100 (MQ=255)
                                                                                      cAACCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtc            >  2:2401548/1‑100 (MQ=255)
                                                                                      cAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtc            >  1:1876624/1‑100 (MQ=255)
                                                                                       aaCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg           >  2:2003802/1‑100 (MQ=255)
                                                                                       aaCCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg           <  2:2164837/100‑1 (MQ=255)
                                                                                         ccACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAg         <  1:1345865/100‑1 (MQ=255)
                                                                                            cGTCCGCAGACGTGTTTGCTATTCATATAACGTCGATTTCAGCTACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcg       <  1:2403035/99‑1 (MQ=255)
                                                                                            cGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcga      >  2:1806981/1‑100 (MQ=255)
                                                                                             gTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgac     <  1:968958/100‑1 (MQ=255)
                                                                                              tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGcgacg    >  1:906547/1‑100 (MQ=255)
                                                                                                cGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGaa  <  1:271864/100‑1 (MQ=255)
                                                                                                cGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGaa  >  2:629829/1‑100 (MQ=255)
                                                                                                   |                                                                                                
AAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAA  >  NC_000913/699460‑699655

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