Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,879:1 +G 19.9% intergenic (‑142/+26) ypfJ ← / ← purC putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,8781.G19.9% 421.0 / 124.6 146intergenic (‑141/+27)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (51/66);  new base G (12/17);  total (63/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

Only 100 of 152 total aligned reads displayed.

TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977
                                                                                                    |                                                                                                   
tggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                      >  1:2747059/1‑100 (MQ=255)
 ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                    >  1:1504568/1‑100 (MQ=255)
  gtTCAGCGGAAGATTACCACAGGAAACCGGGATGAAATGTCGGAGGAATGACCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTACCGTGCCAGCT‑g                                                                                                    >  2:2408916/1‑99 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                                   >  2:1333348/1‑100 (MQ=255)
  gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                                   >  2:494723/1‑100 (MQ=255)
    tCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCAc                                                                                                 >  1:211375/1‑100 (MQ=255)
     cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                <  2:3273092/100‑1 (MQ=255)
     cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATAGATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                <  1:2356099/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                               <  2:675487/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                               >  2:253686/1‑100 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                               <  1:4060052/100‑1 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                >  1:173569/1‑99 (MQ=255)
      agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCAGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                                <  1:53460/99‑1 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                              >  2:3889558/1‑100 (MQ=255)
       gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGc                                                                                              <  1:3527836/100‑1 (MQ=255)
        agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                             <  2:2059640/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTTAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCa                                                                                             >  1:4368678/1‑99 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                            <  1:2856594/100‑1 (MQ=255)
         ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCaa                                                                                            <  2:2872847/100‑1 (MQ=255)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                           <  2:2954893/100‑1 (MQ=39)
          gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                           <  1:1158992/100‑1 (MQ=39)
           aacaTTACCCCAACAAACCGGTATGAAATGGCGGAGGAATCAGCGTTCAACGCATATATTCCGATAGCGTTTGTCTTGCCTTG‑CAGCTGGCACGGCa                                                                                             <  1:3563226/94‑1 (MQ=255)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAga                                                                                          >  2:3251921/1‑100 (MQ=255)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAga                                                                                          <  2:2273214/100‑1 (MQ=255)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                          <  2:2058543/100‑1 (MQ=39)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                          <  1:2132609/100‑1 (MQ=39)
           aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                          <  1:3652805/100‑1 (MQ=39)
            aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                         >  2:1592512/1‑100 (MQ=255)
              aTTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatga                                                                                       <  2:1404729/100‑1 (MQ=255)
                tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                                     <  1:2278360/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                                    <  2:3922605/100‑1 (MQ=255)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgacaaccg                                                                                    <  1:2759877/100‑7 (MQ=37)
                 aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgacaac                                                                                      <  1:3893976/98‑5 (MQ=38)
                     caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAg                                                                                >  2:2950401/1‑100 (MQ=255)
                       ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaa                                                                              >  2:2893845/1‑100 (MQ=255)
                          aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                           <  2:448966/100‑1 (MQ=255)
                          aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                           <  2:4214741/100‑1 (MQ=255)
                           aGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCTAGCT‑GCACGGCAAGATGCTGAACAGAAAAAt                                                                          >  2:1792833/1‑100 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGGGCGTTCAATTCATATATCCCGAGCGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGACGAACAGAAAAATc                                                                         <  1:26380/100‑1 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                         <  2:1184713/100‑1 (MQ=255)
                            gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                         >  2:642203/1‑100 (MQ=255)
                                    aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                 >  2:1764450/1‑100 (MQ=255)
                                      aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                               <  1:1082956/100‑1 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGtt                                                             <  2:1980990/100‑2 (MQ=255)
                                        ggCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTa                                                             <  2:1495367/100‑1 (MQ=255)
                                         gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                            >  2:2019309/1‑100 (MQ=255)
                                          cGGAGGAATGAGCGTTCAATGCAAATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                           >  1:18069/1‑100 (MQ=255)
                                           ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTacac                                                          <  2:4404960/100‑1 (MQ=255)
                                            gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAcc                                                         <  2:2748974/100‑1 (MQ=255)
                                             aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                        <  2:2735499/100‑1 (MQ=255)
                                              ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                       <  2:2796951/100‑1 (MQ=255)
                                                aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                                     <  2:2122328/100‑1 (MQ=255)
                                                  tGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcg                                                   <  2:3371583/100‑1 (MQ=255)
                                                    aGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcggc                                                 <  1:2228579/100‑1 (MQ=255)
                                                     gCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCGGgcggcg                                                <  2:2315511/100‑1 (MQ=255)
                                                        ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             <  2:3962444/100‑1 (MQ=255)
                                                          cAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                             >  2:4175895/1‑98 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  1:3057666/1‑100 (MQ=255)
                                                           aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGAc                                          >  2:3802741/1‑100 (MQ=255)
                                                            aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                         >  1:3350926/1‑100 (MQ=255)
                                                             ggCATATGTTCCGAGAGGGGATGTCGTGCCGTGACAGCG‑GCACGGCAATATGATGACCAGAAAAATGAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                        <  1:3211163/99‑1 (MQ=255)
                                                              gatcttatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  1:840085/95‑1 (MQ=255)
                                                              gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                       <  2:492286/100‑1 (MQ=255)
                                                               cATATATTCCGAGATCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                                      <  1:679548/100‑1 (MQ=255)
                                                                atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                     >  2:3358550/1‑100 (MQ=255)
                                                                atatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtn                                     <  1:941152/100‑2 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  1:974249/100‑1 (MQ=255)
                                                                 tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                    <  2:3908581/100‑1 (MQ=255)
                                                                  atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAGGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                                   >  1:1614059/1‑100 (MQ=255)
                                                                   tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGcagatc                                  >  2:840088/1‑95 (MQ=255)
                                                                    atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                 >  2:2714938/1‑100 (MQ=255)
                                                                     ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                <  1:4135296/100‑1 (MQ=255)
                                                                      tCCGAGGGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                                <  2:469214/99‑1 (MQ=255)
                                                                      tCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                               >  1:665193/1‑100 (MQ=255)
                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                             >  1:3241114/1‑100 (MQ=255)
                                                                        cGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCt                             <  2:3582357/100‑1 (MQ=255)
                                                                          agagCGGTTGTCTTGCCGTTCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                           >  1:1106317/1‑100 (MQ=255)
                                                                           gagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                          <  2:62488/100‑1 (MQ=255)
                                                                            agCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGa                         <  2:1387086/100‑1 (MQ=255)
                                                                             gCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGAt                        <  1:3906854/100‑1 (MQ=255)
                                                                              tcatcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  2:531289/7‑100 (MQ=37)
                                                                              cGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                       >  1:3613167/1‑100 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  2:691063/100‑1 (MQ=255)
                                                                               ggTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCa                      <  2:525540/100‑1 (MQ=255)
                                                                               ggTTGCCTTGCCGTGCCAGCT‑GCACGTCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCgcac                       <  1:1302912/99‑4 (MQ=255)
                                                                                 ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                    <  2:3040771/100‑1 (MQ=255)
                                                                                   atCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGcc                  >  2:3119400/2‑100 (MQ=39)
                                                                                    tCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:4384593/1‑100 (MQ=255)
                                                                                    tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa                 >  2:2143205/1‑100 (MQ=39)
                                                                                     cTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTGCAGCTTTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc                >  1:2598471/1‑100 (MQ=38)
                                                                                      ttGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  1:742127/1‑100 (MQ=255)
                                                                                      ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc               >  1:1883874/1‑100 (MQ=255)
                                                                                       tGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  2:4133253/1‑100 (MQ=255)
                                                                                       tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCg              >  2:4111783/1‑100 (MQ=255)
                                                                                         ccGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACTGCTTCATAGGCTTCGATCAGGCCACCGAg            <  1:2466999/100‑1 (MQ=255)
                                                                                           gTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGc          <  1:1528005/100‑1 (MQ=255)
                                                                                             gCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc        <  1:2208668/100‑1 (MQ=255)
                                                                                             gCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc        >  2:1148053/1‑100 (MQ=255)
                                                                                                  cT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCg   <  2:3071737/100‑1 (MQ=255)
                                                                                                   t‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCgg  <  1:2262334/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGG  >  NC_000913/2596779‑2596977

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: ‑