Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,151,972 C→T 7.6% Q231* (CAG→TAG) ‡ edd → phosphogluconate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,151,9720CT 7.6% 294.2 / 11.1 105Q231* (CAG→TAG) ‡eddphosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (41/55);  new base T (4/4);  total (45/60)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 115 total aligned reads displayed.

TGGTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGA  >  NC_002947/1151823‑1152120
                                                                                                                                                      |                                                                                                                                                    
 ggTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACaccaacca                                                                                                                                                     <  1:1612582/150‑1 (MQ=255)
  gTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACAAg                                                                                                                                                    <  1:1664532/150‑1 (MQ=255)
       cGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACAACAGCCTCATAACCTGCACCTTCTACGTCACCGCCAACACCAACCAGTTGGt                                                                                                                                               >  2:145271/1‑150 (MQ=255)
         gCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGa                                                                                                                                             <  1:2397559/150‑1 (MQ=255)
           aTCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATg                                                                                                                                           <  1:2328787/150‑1 (MQ=255)
              tcCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGaa                                                                                                                                        <  2:2184828/150‑1 (MQ=255)
              tcCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCAAACACCAACCAGTTGGTGATGGaa                                                                                                                                        <  1:852214/150‑1 (MQ=255)
               ccAACCAGCAGAAGGCCGATGTGCCCCAGCGCTATGCCGACCGCGAGCACCAGCCAGAAAATGGAATGGGGTGCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAg                                                                                                                                       <  1:280554/150‑1 (MQ=255)
                 aacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGCGCTGTTGGAGTCCGAGAAGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTg                                                                                                                                     >  1:290615/1‑150 (MQ=255)
                 aacaaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACGGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTg                                                                                                                                     >  1:2439087/1‑150 (MQ=255)
                   caaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGAt                                                                                                                                   >  2:888336/1‑150 (MQ=255)
                   caaGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCACGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGAt                                                                                                                                   <  1:1787782/150‑1 (MQ=255)
                       cAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATgggg                                                                                                                               <  1:202558/150‑1 (MQ=255)
                        aGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACCCCACCCAGTTCGTGATGGAAGTGATGGGGc                                                                                                                              <  2:1971626/150‑1 (MQ=255)
                         gAAGGCTGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTg                                                                                                                                           <  1:722789/136‑1 (MQ=255)
                         gAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTgttg                                                                                                                                           >  2:722789/1‑133 (MQ=255)
                             gCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCAc                                                                                                                         <  1:1244468/150‑1 (MQ=255)
                               cGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTgg                                                                                                                                          >  2:618574/1‑131 (MQ=255)
                               cGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCATCACCAACTGGTTGGTGTTgg                                                                                                                                          <  1:618574/131‑1 (MQ=255)
                                gATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTg                                                                                                                      >  2:2216361/1‑150 (MQ=255)
                                 aTGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGc                                                                                                                     <  1:865131/150‑1 (MQ=255)
                                 aTGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGATGCAccaaa                                                                                                                     <  2:1174517/150‑4 (MQ=255)
                                   gtgCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCg                                                                                                                   <  2:1825059/150‑1 (MQ=255)
                                     gcgcCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCggg                                                                                                                 >  1:1398989/1‑150 (MQ=255)
                                           gcgcTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGTCACCGCCAACACCAGCCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtc                                                                                                           <  1:1227080/150‑1 (MQ=255)
                                           gcgcTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGTCACCGCCAACACCAGCCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtc                                                                                                           <  1:1227038/150‑1 (MQ=255)
                                            cgcTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTGCGAGTCCGAGATGTATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcg                                                                                                          <  2:665881/150‑1 (MQ=255)
                                             gcTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTg                                                                                                                                                 >  1:929357/1‑110 (MQ=255)
                                             gcTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTg                                                                                                                                                 <  2:929357/110‑1 (MQ=255)
                                              cTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGACTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCtcgtt                                                                                                        <  1:725284/150‑1 (MQ=255)
                                                  gCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTc                                                                                                    <  2:2430834/150‑1 (MQ=255)
                                                   ccGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCACCGACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTAGTCa                                                                                                   >  1:158315/1‑150 (MQ=255)
                                                    cGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCaa                                                                                                  >  2:909897/1‑150 (MQ=255)
                                                     gAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACAAAGTCTTCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAAc                                                                                                 >  1:1332763/1‑150 (MQ=255)
                                                       aggcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAAccc                                                                                               >  2:1057970/1‑150 (MQ=255)
                                                         gcaaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGt                                                                                             >  1:454106/1‑150 (MQ=255)
                                                          caaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCTCCTGCACCTTCTACGGCACCGCCAACACCAAACCGATGATGCTGCAGTGGGTGGGGGTGTACCTGCCGGGCGCCTCGTTCGTCAACCCGTa                                                                                            >  2:1241976/1‑150 (MQ=255)
                                                          caaggcCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTa                                                                                            <  2:1896209/150‑1 (MQ=255)
                                                                  aGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACAcgccg                                                                                     >  1:2283433/1‑150 (MQ=255)
                                                                    ccGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCt                                                                                   >  1:1164339/1‑150 (MQ=255)
                                                                        gAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACCCATACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGt                                                                               <  1:759125/150‑1 (MQ=255)
                                                                        gAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGt                                                                               <  1:2419432/150‑1 (MQ=255)
                                                                           gagCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGAc                                                                            <  1:991664/150‑1 (MQ=255)
                                                                           gagCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGAc                                                                            <  1:888336/150‑1 (MQ=255)
                                                                           gagCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCAAATCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGAc                                                                            <  1:1308987/150‑1 (MQ=255)
                                                                            agCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCAGCGGCGTgtac                                                                           <  1:2098834/150‑4 (MQ=255)
                                                                             gCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCAGCGGCGTgtacg                                                                          <  2:1277331/150‑5 (MQ=255)
                                                                              cTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCa                                                                         >  2:1585697/1‑150 (MQ=255)
                                                                                gTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCAGCGGCGTgtacgggt                                                                       <  1:458352/150‑8 (MQ=255)
                                                                                 ttGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                  <  2:1545787/90‑1 (MQ=255)
                                                                                 ttGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATg                                                                                                                                  >  1:1545787/1‑90 (MQ=255)
                                                                                   ggAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAc                                                                    >  1:857799/1‑150 (MQ=255)
                                                                                   ggAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAc                                                                    >  1:1633083/1‑150 (MQ=255)
                                                                                       tCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCGCTCAccgcc                                                                >  1:236014/1‑150 (MQ=255)
                                                                                        ccGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTGCACGCCGCTGCGTGACGCACTCAccgccg                                                               <  1:1241976/150‑1 (MQ=255)
                                                                                         cGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGa                                                              <  2:290615/150‑1 (MQ=255)
                                                                                         cGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGa                                                              <  2:1560932/150‑1 (MQ=255)
                                                                                         cGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGa                                                              <  1:1192627/150‑1 (MQ=255)
                                                                                         cGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGa                                                              >  2:2286563/1‑150 (MQ=255)
                                                                                          gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                             >  1:1053270/1‑150 (MQ=255)
                                                                                          gagaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                             >  1:84935/1‑150 (MQ=255)
                                                                                          gagaTGAATTCCTACCACAGCCCAGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGGTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAg                                                             >  2:25615/1‑150 (MQ=255)
                                                                                           agaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAc                                                                    >  1:2161668/1‑142 (MQ=255)
                                                                                           agaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAc                                                                    <  2:2161668/142‑1 (MQ=255)
                                                                                           agaTGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAgg                                                            <  1:575026/150‑1 (MQ=255)
                                                                                               gAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTgc                                                        >  2:540131/1‑150 (MQ=255)
                                                                                                  ttCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTGTTggcgg                                                                                                                                       <  1:1678368/68‑4 (MQ=37)
                                                                                                  ttCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACTGGTTGGTGTTggcgg                                                                                                                                       >  2:1678368/1‑65 (MQ=37)
                                                                                                     cTACCACGGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcagca                                                  <  2:2439087/150‑1 (MQ=255)
                                                                                                     cTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcagca                                                  <  2:2016635/150‑1 (MQ=255)
                                                                                                     cTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcagca                                                  >  1:2016635/1‑150 (MQ=255)
                                                                                                      tACCACAGCCCCGGCACCTGCACCTTCTACGTCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGcagcag                                                 >  1:1883335/1‑150 (MQ=255)
                                                                                                           caGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGAc                                            >  2:611385/1‑150 (MQ=255)
                                                                                                                     acctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGa                                             <  2:807804/139‑1 (MQ=255)
                                                                                                                     acctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGa                                             >  1:807804/1‑139 (MQ=255)
                                                                                                                       ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGCGGCTGCGCAGCAGGTGACGCGCATGACCaa                                <  1:2286563/150‑1 (MQ=255)
                                                                                                                       ctgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCaa                                <  2:1398989/150‑1 (MQ=255)
                                                                                                                        tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                               <  1:2216361/150‑1 (MQ=255)
                                                                                                                        tgcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGCCTTAGTGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                               <  1:145271/150‑1 (MQ=255)
                                                                                                                        tgcacctTCTACGGCACCGACAACCCCAACCTGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAg                               <  2:857799/150‑1 (MQ=255)
                                                                                                                         gcacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAgg                              >  2:1752045/1‑150 (MQ=255)
                                                                                                                          cacctTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGc                             <  1:1585697/150‑1 (MQ=255)
                                                                                                                              ttCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCGGGTGACGCGCATGa                                    <  2:2259005/139‑1 (MQ=255)
                                                                                                                              ttCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCGGGTGACGCGCATGa                                    >  1:2259005/1‑139 (MQ=255)
                                                                                                                               tCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCTGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCg                        <  1:25615/150‑1 (MQ=255)
                                                                                                                                 tACGGCACCGCCAACACCAACCAAATGACGACGACAGTGTTCAGTAAGCACGGGGCGGGCTGCTCTGTGGTCGGCGCGTACTCGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGc                      >  2:1466572/1‑150 (MQ=255)
                                                                                                                                  aCGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCAccgccg                                                               >  2:1234523/1‑108 (MQ=255)
                                                                                                                                         cGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcagca                                                  <  2:455664/114‑1 (MQ=255)
                                                                                                                                         cGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCgcagca                                                  >  1:455664/1‑114 (MQ=255)
                                                                                                                                           ccATCACCAACTGGTTGGTGTTGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGc            <  1:3830/150‑1 (MQ=255)
                                                                                                                                           ccAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGc            <  1:1057970/150‑1 (MQ=255)
                                                                                                                                            cAACACCAAACAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGa                                                                             <  2:1207878/84‑1 (MQ=255)
                                                                                                                                            cAACACCAAACAGTTGGTGAGGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGa                                                                             >  1:1207878/1‑84 (MQ=255)
                                                                                                                                                accaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTg       >  2:815510/1‑150 (MQ=255)
                                                                                                                                                 ccaaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTgg      <  2:1633083/150‑1 (MQ=255)
                                                                                                                                                  caaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTTGTAAAGGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTggg     >  2:1901270/1‑150 (MQ=255)
                                                                                                                                                  caaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGAc                                   >  2:2430175/1‑120 (MQ=255)
                                                                                                                                                  caaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTggg     >  2:681888/1‑150 (MQ=255)
                                                                                                                                                  caaccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAACCTCAGCGGCAGCTTCATGCCGCTggg     <  1:323776/150‑1 (MQ=255)
                                                                                                                                                     ccaGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGa  >  2:2070766/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGGTTTCCGGCATCTCCAACAAGCAGAAGGCCGATGTGCGCCAGCGCTATGCCGAAGGCAAGGCCAGCCGCGAAGAGCTGTTGGAGTCCGAGATGAATTCCTACCACAGCCCCGGCACCTGCACCTTCTACGGCACCGCCAACACCAACCAGTTGGTGATGGAAGTGATGGGGCTGCACCTGCCGGGCGCCTCGTTCGTCAACCCGTACACGCCGCTGCGTGACGCACTCACCGCCGAGGCTGCGCAGCAGGTGACGCGCATGACCAAGGCCAGCGGCAGCTTCATGCCGCTGGGCGA  >  NC_002947/1151823‑1152120

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