Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,737,713 (G)6→7 intergenic (+67/‑32) pheL → / → pheA pheA gene leader peptide/chorismate mutase and prephenate dehydratase, P‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,737,7071.G98.4% 226.7 / ‑2.5 63intergenic (+61/‑38)pheL/pheApheA gene leader peptide/chorismate mutase and prephenate dehydratase, P‑protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (34/28);  total (35/28)
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.50e-01

AACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑GG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTG  >  NC_000913/2737557‑2737849
                                                                                                                                                       |                                                                                                                                               
aaCAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATc                                                                                                                                                  <  1:719301/151‑1 (MQ=255)
aaCAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATc                                                                                                                                                  >  2:230271/1‑151 (MQ=255)
 aCAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑g                                                                                                                                                >  2:221520/1‑151 (MQ=255)
  cAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg                                                                                                                                               <  2:212379/151‑1 (MQ=255)
    gTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑gg                                                                                                                                            <  1:526602/151‑1 (MQ=255)
      aCGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑gggg                                                                                                                                          <  2:474/151‑1 (MQ=255)
      aCGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑ggg                                                                                                                                           <  1:549085/150‑1 (MQ=255)
       cGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑ggggg                                                                                                                                         <  1:238899/151‑2 (MQ=255)
          gtactgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑gggggcct                                                                                                                                      <  2:613952/151‑5 (MQ=255)
            actgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑gggggcctt                                                                                                                                     <  2:184099/150‑6 (MQ=255)
            actgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑gg‑gggggcctt                                                                                                                                     <  1:801479/150‑6 (MQ=255)
               gtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGTATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCtttttt                                                                                                                                  <  2:604754/151‑1 (MQ=255)
                   tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTg                                                                                                                              >  2:1035490/1‑151 (MQ=255)
                       gTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATa                                                                                                                           >  1:15089/1‑150 (MQ=255)
                        tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAAc                                                                                                                         >  2:296943/1‑151 (MQ=255)
                        tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAAc                                                                                                                         >  2:298428/1‑151 (MQ=255)
                        tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAAc                                                                                                                         >  1:125178/1‑151 (MQ=255)
                        tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAAc                                                                                                                         >  2:39750/1‑151 (MQ=255)
                          aCTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACa                                                                                                                        >  2:1013871/1‑150 (MQ=255)
                          aCTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACa                                                                                                                        >  1:523105/1‑150 (MQ=255)
                             tAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACaaaa                                                                                                                     >  1:1033775/1‑150 (MQ=255)
                              aaGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCTTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACaaaaa                                                                                                                    <  2:330011/150‑1 (MQ=255)
                               aGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAgg                                                                                                                  <  2:485362/151‑1 (MQ=255)
                                  aaacaaacATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCaa                                                                                                               <  2:1172283/151‑1 (MQ=255)
                                   aacaaacaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAac                                                                                                              <  2:171894/151‑1 (MQ=255)
                                      aaacaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACt                                                                                                           <  2:1018928/151‑1 (MQ=255)
                                                cacaTACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCgg                                                                                                 >  2:1007411/1‑151 (MQ=255)
                                                          ttttCTTCGCATTCTTTTTTACGTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGtt                                                                                       >  1:798612/1‑151 (MQ=255)
                                                           tttCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGtt                                                                                       <  2:773587/150‑1 (MQ=255)
                                                           tttCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGtt                                                                                       >  2:483567/1‑150 (MQ=255)
                                                            ttcttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTa                                                                                      <  1:127292/150‑1 (MQ=255)
                                                            ttcttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTa                                                                                      <  1:392167/150‑1 (MQ=255)
                                                                  gCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGgcgc                                                                               >  1:342215/1‑151 (MQ=255)
                                                                    aTTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTg                                                                             >  1:996454/1‑151 (MQ=255)
                                                                      tCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCg                                                                           >  1:137178/1‑151 (MQ=255)
                                                                         tttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCg                                                                           >  1:255563/1‑148 (MQ=255)
                                                                          ttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑GGAGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgag                                                                         >  2:619923/1‑149 (MQ=255)
                                                                            ttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgagaga                                                                      <  1:777693/150‑1 (MQ=255)
                                                                              aCCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaaa                                                                   <  2:366574/151‑1 (MQ=255)
                                                                              aCCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaa                                                                    >  2:1165152/1‑150 (MQ=255)
                                                                               ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaaa                                                                   <  1:925080/150‑1 (MQ=255)
                                                                               ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt                                                                  <  2:735418/151‑1 (MQ=255)
                                                                                 ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATc                                                                 >  1:220986/1‑150 (MQ=255)
                                                                                  tCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATc                                                                 <  2:209919/149‑1 (MQ=255)
                                                                                  tCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCa                                                                <  1:487558/150‑1 (MQ=255)
                                                                                  tCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGAGAGAGAAAATCa                                                                <  1:1167134/150‑1 (MQ=255)
                                                                                    cccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcg                                                             >  2:277609/1‑151 (MQ=255)
                                                                                     ccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcgc                                                            >  1:1072588/1‑151 (MQ=255)
                                                                                       tGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcgcgc                                                          >  2:280843/1‑151 (MQ=255)
                                                                                         aaTGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTg                                                        <  1:514184/151‑1 (MQ=255)
                                                                                                ggCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGaaa                                                 <  2:91291/151‑1 (MQ=255)
                                                                                                  cGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGaaaaa                                               <  2:714799/151‑1 (MQ=255)
                                                                                                      tcgtcgTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCTAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAattat                                           >  1:891600/1‑151 (MQ=255)
                                                                                                        gtcgtGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAAttatta                                         <  1:221517/151‑1 (MQ=255)
                                                                                                          cgtGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGc                                       >  1:318944/1‑151 (MQ=255)
                                                                                                             gtgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGt                                     <  1:17868/150‑1 (MQ=255)
                                                                                                             gtgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGt                                     <  1:296940/150‑1 (MQ=255)
                                                                                                               gAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTAc                                  <  2:620975/151‑1 (MQ=255)
                                                                                                                aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACt                                 <  2:338825/151‑1 (MQ=255)
                                                                                                                aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACt                                 <  2:703937/151‑1 (MQ=255)
                                                                                                                 aaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTg                                <  1:711266/151‑1 (MQ=255)
                                                                                                                  aCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTgg                               >  2:1057669/1‑151 (MQ=255)
                                                                                                                   cAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGc                              <  2:799662/151‑1 (MQ=255)
                                                                                                                    aGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGc                              >  2:251468/1‑150 (MQ=255)
                                                                                                                     gAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGc                              >  1:87706/1‑149 (MQ=255)
                                                                                                                        tGCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAAc                         >  1:848911/1‑151 (MQ=255)
                                                                                                                         gCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAAc                         <  2:437088/150‑1 (MQ=255)
                                                                                                                         gCGAAGACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACg                        >  2:15605/1‑151 (MQ=255)
                                                                                                                              gACGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGcgcg                   >  2:1187789/1‑151 (MQ=255)
                                                                                                                                cGAACAATAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGa                  >  1:89177/1‑150 (MQ=255)
                                                                                                                                     aaTAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGc            >  1:1185818/1‑151 (MQ=255)
                                                                                                                                       tAAGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCg          >  2:90296/1‑151 (MQ=255)
                                                                                                                                        aaGGCCTCCCAAATCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTAGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGt         >  1:1023262/1‑151 (MQ=255)
                                                                                                                                                  aaaTCGGG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGCAAAATTATTAGCGTTACTGGCAGAACGGCGCGAAATGGCCGTCGAGGTg  >  2:1114849/1‑148 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
AACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATC‑GG‑GGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTG  >  NC_000913/2737557‑2737849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑