Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,446,579 T→C T299A (ACC→GCC)  aroC ← chorismate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,446,5790TC98.4% 206.3 / ‑5.9 61T299A (ACC→GCC) aroCchorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (36/24);  minor base A (0/1);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAG  >  NC_000913/2446442‑2446714
                                                                                                                                         |                                                                                                                                       
cAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCaa                                                                                                                                       <  1:552555/140‑1 (MQ=255)
     gATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTg                                                                                                                                  >  1:1348273/1‑140 (MQ=255)
       tCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGa                                                                                                                                >  1:450030/1‑140 (MQ=255)
                aCGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGttt                                                                                                                       <  2:450030/140‑1 (MQ=255)
                 cGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTc                                                                                                                      <  2:736249/140‑1 (MQ=255)
                   aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                    >  1:77118/1‑140 (MQ=255)
                   aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAg                                                                                                                    <  2:1295157/140‑1 (MQ=255)
                     cGCCAGCATCGCTTCTGCGATCTGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCa                                                                                                                     >  1:357601/1‑137 (MQ=255)
                     cGCCAGCATCGCTTCTGCGATCTGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCa                                                                                                                     <  2:357601/137‑1 (MQ=255)
                         aGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCat                                                                                                              <  2:1278490/140‑1 (MQ=255)
                              cGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATggg                                                                                                         <  2:1477881/140‑1 (MQ=255)
                                     gCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGAt                                                                                                  <  2:1348273/140‑1 (MQ=255)
                                     gCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGAt                                                                                                  >  1:1240498/1‑140 (MQ=255)
                                        aTCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtg                                                                                               >  1:63083/1‑140 (MQ=255)
                                           ggCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgctg                                                                                            >  1:953071/1‑140 (MQ=255)
                                           ggCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgctg                                                                                            >  1:436783/1‑140 (MQ=255)
                                           ggCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgctg                                                                                            >  1:450263/1‑140 (MQ=255)
                                            gCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTtgctgc                                                                                           <  2:1240498/140‑1 (MQ=255)
                                             cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                          <  2:928289/140‑1 (MQ=255)
                                             cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                          <  2:600626/140‑1 (MQ=255)
                                             cACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGcc                                                                                          <  2:63083/140‑1 (MQ=255)
                                                tGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCgc                                                                                       >  2:638584/1‑140 (MQ=255)
                                                  cgcgGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCgctg                                                                                     >  1:1424770/1‑140 (MQ=255)
                                                       aTCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGAt                                                                                >  2:802739/1‑140 (MQ=255)
                                                               acaGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCCGTTAATGGTCCGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCga                                                                        <  2:1357763/140‑1 (MQ=255)
                                                                 aGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCgaga                                                                      >  2:1119601/1‑140 (MQ=255)
                                                                 aGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCgaga                                                                      >  2:1241033/1‑140 (MQ=255)
                                                                   ggATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                    >  2:1290698/1‑140 (MQ=255)
                                                                   ggATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                    >  1:524457/1‑140 (MQ=255)
                                                                         tGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgccgcc                                                              >  2:814797/1‑140 (MQ=255)
                                                                           aCGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGccgcccg                                                            >  2:418914/1‑140 (MQ=255)
                                                                             ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgc                                                                  >  1:967562/1‑132 (MQ=255)
                                                                             ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCa                                                          >  2:1155185/1‑140 (MQ=255)
                                                                             ggCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGAAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgc                                                                  <  2:967562/132‑1 (MQ=255)
                                                                               ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                        >  1:1285845/1‑140 (MQ=255)
                                                                               ccTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATg                                                        >  2:777310/1‑140 (MQ=255)
                                                                                cTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATgg                                                       <  1:1155185/140‑1 (MQ=255)
                                                                                  ttGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGtt                                                     <  1:638584/140‑1 (MQ=255)
                                                                                    ggTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGc                                                   >  1:370761/1‑140 (MQ=255)
                                                                                           atcTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGaa                                            <  2:1424770/140‑1 (MQ=255)
                                                                                             ctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAAc                                          >  1:391770/1‑140 (MQ=255)
                                                                                                     cttcGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTg                                  >  1:502755/1‑140 (MQ=255)
                                                                                                       tcGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGt                                >  2:271212/1‑140 (MQ=255)
                                                                                                       tcGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGt                                >  1:346120/1‑140 (MQ=255)
                                                                                                       tcGCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGt                                >  2:657847/1‑140 (MQ=255)
                                                                                                         gCCAAAGCGGTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGa                              <  1:418914/140‑1 (MQ=255)
                                                                                                                 ggTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATcg                      <  1:271212/140‑1 (MQ=255)
                                                                                                                 ggTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATcg                      <  2:953071/140‑1 (MQ=255)
                                                                                                                 ggTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATcg                      <  2:436783/140‑1 (MQ=255)
                                                                                                                 ggTTAATGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATcg                      <  2:450263/140‑1 (MQ=255)
                                                                                                                      ttGGTACGACCCGGCAGGGCAAGGCTGGGGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGtt                 <  2:1091156/139‑1 (MQ=255)
                                                                                                                       tGGTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTc                >  2:914060/1‑140 (MQ=255)
                                                                                                                         gTACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTg              >  2:875036/1‑140 (MQ=255)
                                                                                                                          tACGACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTgg             >  1:22395/1‑140 (MQ=255)
                                                                                                                             gACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTg          >  2:655728/1‑140 (MQ=255)
                                                                                                                             gACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTg          >  2:506815/1‑140 (MQ=255)
                                                                                                                             gACCCGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTAGGTGATTTCATCGCGGTTCTGGCTg          >  1:429868/1‑140 (MQ=255)
                                                                                                                               cccGGCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGcc        >  1:893115/1‑140 (MQ=255)
                                                                                                                                 cggcacggcaATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgc      <  2:1285845/140‑1 (MQ=255)
                                                                                                                                 cggcacggcaATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgc      >  2:625159/1‑140 (MQ=255)
                                                                                                                                 cggcacggcaATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgc      >  2:433206/1‑140 (MQ=255)
                                                                                                                                  ggCACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgcg     <  1:625159/140‑1 (MQ=255)
                                                                                                                                    cACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCa   <  1:655728/140‑1 (MQ=255)
                                                                                                                                    cACGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCa   >  1:1309041/1‑140 (MQ=255)
                                                                                                                                     aCGGCAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAg  <  2:1090600/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATCCCGACACAGGGATCGTGACGGCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAG  >  NC_000913/2446442‑2446714

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