Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,874 Δ1 bp 15.3% intergenic (‑137/+31) 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,8730C.15.3% 398.1 / 74.3 137intergenic (‑136/+32)ypfJ/purCputative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/52);  new base . (11/10);  total (74/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 141 total aligned reads displayed.

GTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCT  >  NC_000913/2596774‑2596970
                                                                                                   |                                                                                                  
gTAGGTGGTTCAGAGGAAGATTACCTCAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                    >  2:2301654/1‑100 (MQ=255)
gTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                    >  1:2237034/1‑100 (MQ=255)
 tAGGTGGATCAAAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGcc                                                                                                   >  2:2184885/1‑100 (MQ=255)
  aGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCa                                                                                                  <  2:639979/100‑1 (MQ=255)
    gtggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGc                                                                                                >  1:2883487/1‑100 (MQ=255)
      ggtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGCGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                             >  2:483797/1‑100 (MQ=255)
       gtggagagGAAGATTACCAGAGGAAAGCGGGATGAAATGGGGGGGGAATGCGCGGTCAATGAATATATTCAGAGAGCGGTTGGCTTGCCGTGcga                                                                                                  <  2:230498/91‑3 (MQ=255)
       gtTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑Gc                                                                                            >  1:927986/1‑100 (MQ=255)
          cAGAGGAAGATTACCACCGGAAAGCGGGATGACATGGCGGAGGAATGCGCGTTCACTGCATATATTCCGAGCGCGGTTGTCTTGCCGTGACAGCT‑GAACg                                                                                         >  2:850198/1‑100 (MQ=255)
          cAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCAAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                         <  2:3223852/100‑1 (MQ=255)
           agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACgg                                                                                        <  1:4016413/100‑1 (MQ=255)
           agagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACg                                                                                         <  1:1572029/99‑1 (MQ=255)
            gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGc                                                                                       >  2:2359378/1‑100 (MQ=255)
             agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCa                                                                                      <  1:3432255/100‑1 (MQ=255)
              ttccgatctccACAGGAAAGCGGGATGAAATGGCCGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGATGGCACGGCaa                                                                                     <  2:2191190/91‑1 (MQ=255)
              ggAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                     <  1:2135970/100‑1 (MQ=255)
              ggAAGATTAACACAGGAAAGCGGGATGAAATGGCGCAGCAATGAGCGTTCAATGCATATATTCCGAGAGGGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaa                                                                                     <  2:1856338/100‑1 (MQ=255)
               gAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAg                                                                                    <  2:2764994/100‑1 (MQ=39)
                aaGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAga                                                                                   <  1:1642133/100‑1 (MQ=39)
                 aGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGACCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgat                                                                                  >  2:1042850/1‑100 (MQ=255)
                  gATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAgatg                                                                                 <  1:799984/100‑1 (MQ=255)
                  gATTAACACAGAAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTAAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCAAgaca                                                                                 <  2:2483844/100‑3 (MQ=37)
                     tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGatgatg                                                                              <  1:397193/100‑1 (MQ=255)
                     tACCACACGAAAGCGGGATGAAATGGCTGATGAATGCGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCTTGCCAGCT‑GCACGGCAAGatgatg                                                                              <  1:3970458/100‑1 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                             <  1:3522248/100‑1 (MQ=255)
                      aCCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGAtgatga                                                                             <  2:2453928/100‑1 (MQ=255)
                       ccACAGGGAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGCACGGCaaagatcggaa                                                                            >  1:2191204/1‑91 (MQ=255)
                       ccACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGaa                                                                            >  1:23764/1‑100 (MQ=255)
                          caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAg                                                                         <  2:2074587/100‑1 (MQ=255)
                           aGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGa                                                                        >  1:173502/1‑100 (MQ=255)
                            ggAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaa                                                                       >  2:3549044/1‑100 (MQ=255)
                             gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaa                                                                      <  1:1634593/100‑1 (MQ=255)
                             gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaa                                                                      <  1:603128/100‑1 (MQ=255)
                             gAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaa                                                                      <  1:3716500/100‑1 (MQ=255)
                              aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaa                                                                     <  1:3647855/100‑1 (MQ=255)
                              aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaa                                                                     <  1:1175299/100‑1 (MQ=255)
                               aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                    >  2:1091602/1‑100 (MQ=255)
                               aaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGACGTGCCAGCT‑GCACTGCAAGATGATGAACAGaaaaa                                                                    <  2:1207219/100‑1 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGaaaaa                                                                    >  1:722141/1‑98 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                  <  2:1640352/100‑1 (MQ=255)
                                 gCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATc                                                                  <  1:3345616/100‑1 (MQ=255)
                                    ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGt                                                               >  1:1915523/1‑100 (MQ=255)
                                     gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                              <  1:2543311/100‑1 (MQ=255)
                                     gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTATTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTc                                                              <  2:3354076/100‑1 (MQ=255)
                                        gAAATGGCGGAGGAATGAGCGTTCGATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAg                                                           >  1:3151532/1‑100 (MQ=255)
                                         aaaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                          >  2:322435/1‑100 (MQ=255)
                                          aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                         >  2:3039517/1‑100 (MQ=255)
                                           aTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTg                                                        <  2:2759616/100‑1 (MQ=255)
                                              gCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTac                                                     >  2:783244/1‑100 (MQ=255)
                                               cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                    >  2:1451847/1‑100 (MQ=255)
                                               cGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                    <  1:2336294/100‑1 (MQ=255)
                                                  aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                 <  1:2439425/100‑1 (MQ=255)
                                                  aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                 <  1:3728310/100‑1 (MQ=255)
                                                  aggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACAccc                                                 <  2:2515431/100‑1 (MQ=255)
                                                   ggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                <  1:893606/100‑1 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCa                                                >  1:2402636/1‑98 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                              <  1:1864123/100‑1 (MQ=255)
                                                     aaTGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAgg                                              >  2:2115005/1‑100 (MQ=255)
                                                        gAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAggcgg                                           >  1:2295096/1‑100 (MQ=255)
                                                           cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGCTGAACAGAAAAATCACCCCAGCTGTACACCCa                                                >  2:3973062/1‑92 (MQ=255)
                                                           cGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGcggcgg                                        <  1:2116624/100‑1 (MQ=255)
                                                            gTTCAATGCATATATTCCTAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGCCCAGCTGTACACCCAGGcggcgg                                        >  2:1289732/1‑99 (MQ=255)
                                                            gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcggg                                       <  2:1951/100‑1 (MQ=255)
                                                             ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                      <  1:3241792/100‑1 (MQ=255)
                                                              tCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggcg                                     <  1:323521/100‑1 (MQ=255)
                                                                 aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                  <  1:322438/100‑1 (MQ=255)
                                                                 aTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACg                                  >  2:1930094/1‑100 (MQ=255)
                                                                  tGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                 >  1:2644148/1‑100 (MQ=255)
                                                                   gCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGc                                <  2:871933/100‑1 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACgg                                 >  2:3263590/1‑98 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCt                               <  1:1639095/100‑1 (MQ=255)
                                                                    cATATATTCCGAGAGCGGTTGTCTTACCGTGCCAGCT‑GCACAGCAAGATGATGAACAGAAAAATCACTCCAGCTATACACCCAAGCCGCCAGCGa                                    >  2:491740/1‑95 (MQ=255)
                                                                       atatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCa                            <  1:1894068/100‑1 (MQ=255)
                                                                        tatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCAt                           >  2:866690/1‑100 (MQ=255)
                                                                         atTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                          >  2:3539835/1‑100 (MQ=255)
                                                                          ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAg                         <  2:2237613/100‑1 (MQ=255)
                                                                           tCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAgg                        >  1:317610/1‑100 (MQ=255)
                                                                            ccGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGc                       >  1:1210958/1‑100 (MQ=255)
                                                                               agagCGGTTGTCATTCCTTGACAGCT‑GCACTTCAAGTTGATGAACAGAAAAATCAGTCCAGCTGTACTCCCAGGCTGCGGGCGACGGCTTCATAg                         <  2:3931556/95‑1 (MQ=255)
                                                                                gagCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCg                   >  2:1555599/1‑100 (MQ=255)
                                                                                      ttGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg             <  1:1037227/100‑1 (MQ=255)
                                                                                      tcatCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg             >  1:738614/4‑100 (MQ=38)
                                                                                         tCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa          >  1:332206/1‑100 (MQ=255)
                                                                                         tCTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCa          >  1:564599/1‑100 (MQ=39)
                                                                                          cTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  1:2694277/1‑100 (MQ=255)
                                                                                          cTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:1703938/1‑100 (MQ=255)
                                                                                          cTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCCACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:2303729/1‑100 (MQ=255)
                                                                                          cTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:2536396/1‑100 (MQ=39)
                                                                                          cTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  1:3604180/1‑100 (MQ=39)
                                                                                          cTTGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAc         >  2:2470503/1‑100 (MQ=39)
                                                                                           ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc        >  1:2325123/1‑100 (MQ=255)
                                                                                           ttGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGAGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc        >  1:292647/1‑100 (MQ=255)
                                                                                            tGCCGTG‑CAGCTGGCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTGCGATCAGGCCACCg       >  1:3395662/1‑100 (MQ=255)
                                                                                              ccGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTTATAGGCGTCGATCAGGCCACCGAg     <  2:3466510/100‑1 (MQ=255)
                                                                                              ccGTGCCAGCT‑GCACGCCAAGATGATGAACAGAAAAGTCAGTCCAGCTGTACACCCAGGCGGCGGGAGATGGCTTTATAGGCTTCGATGCGGCCACCGAg     <  2:2400798/100‑1 (MQ=255)
                                                                                               cGTGCCAGCT‑GCACGGCAAGATGATGAACAGAACAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg    <  1:2202416/100‑1 (MQ=255)
                                                                                               cGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg    <  1:1930104/100‑1 (MQ=255)
                                                                                               cGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAgg    >  1:2435892/1‑100 (MQ=255)
                                                                                                gTGCCAGAT‑GCACGGCAAGACGATGAGCAGAAAAATCAGTCCAGCTGTACCCCCAGGCGGCGGGCGACGGCTTCATGGGCTTCGATCAGGCCACCGAgg    <  2:1858582/99‑1 (MQ=255)
                                                                                                 tGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGct  >  1:1568577/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GCACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCT  >  NC_000913/2596774‑2596970

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