Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,737,705 A→T 6.5% intergenic (+59/‑40) 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,7050AT 6.5% 299.3 / 6.7 89intergenic (+59/‑40)pheL/pheApheA gene leader peptide/chorismate mutase and prephenate dehydratase, P‑protein
Reads supporting (aligned to +/- strand):  ref base A (39/43);  new base T (3/3);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGG  >  NC_000913/2737557‑2737850
                                                                                                                                                    |                                                                                                                                                 
aaCAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCaaa                                                                                                                                                   >  2:506325/1‑149 (MQ=255)
 aCAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAAt                                                                                                                                                  <  4:351100/149‑1 (MQ=255)
  cAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATc                                                                                                                                                 <  1:139335/149‑1 (MQ=255)
    gTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgg                                                                                                                                               <  4:338844/149‑1 (MQ=255)
     tACGGGTACTGTACTAAAGTCACTTAAGTAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCggg                                                                                                                                              <  3:686122/149‑1 (MQ=255)
      aCGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgggg                                                                                                                                             <  3:434983/149‑1 (MQ=255)
       cGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAAGGCGAAGACGAACAATAAGGCCTCCCAACTCggggg                                                                                                                                            <  2:633189/149‑1 (MQ=255)
        gggTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgggggg                                                                                                                                           >  1:376248/1‑149 (MQ=255)
          gtactgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACTAACAATAAGGCCCCCCGATTTGGGAGGcc                                                                                                                                         <  4:126211/149‑1 (MQ=255)
          gtactgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCCCCCGATTTGGGAGGcc                                                                                                                                         <  1:488577/149‑1 (MQ=255)
            actgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCCCCCGATTTGGGAGGCCtt                                                                                                                                       <  1:166577/149‑1 (MQ=255)
                tactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCtttttt                                                                                                                                   >  3:425967/1‑149 (MQ=255)
                 actAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTa                                                                                                                                  <  2:267233/149‑1 (MQ=255)
                   tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAtt                                                                                                                                >  2:227007/1‑149 (MQ=255)
                   tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAtt                                                                                                                                <  2:584438/149‑1 (MQ=255)
                     aaGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGAAAGCGTTTCGTCGTGTGAAACAGAATGCGAATACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGa                                                                                                                              >  2:309876/1‑149 (MQ=255)
                      aGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGAt                                                                                                                             >  3:206336/1‑149 (MQ=255)
                        tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATaa                                                                                                                           <  3:275933/149‑1 (MQ=255)
                            ttAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGa                                                                                                                              >  1:173900/1‑142 (MQ=255)
                            ttAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGa                                                                                                                              <  2:173900/142‑1 (MQ=255)
                            ttAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaa                                                                                                                       >  4:372255/1‑149 (MQ=255)
                            ttAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaa                                                                                                                       >  2:375127/1‑149 (MQ=255)
                              aaGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                     <  2:534247/149‑1 (MQ=255)
                              aaGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                     <  4:144340/149‑1 (MQ=255)
                              aaGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                     <  1:366039/149‑1 (MQ=255)
                               aGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGTGGGCCTTTTTTATTGATAACAAAAAg                                                                                                                    >  3:305769/1‑149 (MQ=255)
                                  aaacaaacATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCa                                                                                                                 <  3:219141/149‑1 (MQ=255)
                                      aaacaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAacac                                                                                                             >  1:745740/1‑149 (MQ=255)
                                        acaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                 <  4:207944/127‑1 (MQ=255)
                                        acaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                 >  3:207944/1‑127 (MQ=255)
                                        acaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTa                                                                                                           <  2:68251/149‑1 (MQ=255)
                                         caTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTAt                                                                                                          <  2:128491/149‑1 (MQ=255)
                                           tGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgg                                                                                                                                               <  1:470279/110‑1 (MQ=255)
                                           tGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgg                                                                                                                                               >  2:470279/1‑110 (MQ=255)
                                               acacATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATc                                                                                                    <  4:726698/149‑1 (MQ=255)
                                               acacATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATc                                                                                                    <  3:123548/149‑1 (MQ=255)
                                                cacaTACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCg                                                                                                   >  4:318748/1‑149 (MQ=255)
                                                    tACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGTAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGaaa                                                                                               >  1:407427/1‑149 (MQ=255)
                                                      ccGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAAc                                                                                                                          <  2:444320/120‑1 (MQ=255)
                                                      ccGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAAc                                                                                                                          >  1:444320/1‑120 (MQ=255)
                                                      ccGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAAc                                                                                             >  2:775598/1‑149 (MQ=255)
                                                      ccGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAAc                                                                                             <  3:372255/149‑1 (MQ=255)
                                                        gTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAAccc                                                                                           <  1:315802/149‑1 (MQ=255)
                                                         tttttCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCg                                                                                          <  1:506325/149‑1 (MQ=255)
                                                                   cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAaca                                                                                                              <  1:539844/119‑1 (MQ=255)
                                                                   cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAaca                                                                                                              >  2:539844/1‑119 (MQ=255)
                                                                   cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTAt                                                                                                          <  1:8398/123‑1 (MQ=255)
                                                                   cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTAt                                                                                                          >  2:8398/1‑123 (MQ=255)
                                                                      tCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGc                                                                             >  2:760760/1‑149 (MQ=255)
                                                                          ttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCgaga                                                                         >  1:338008/1‑149 (MQ=255)
                                                                              aCCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaa                                                                     >  4:678600/1‑149 (MQ=255)
                                                                               ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaaa                                                                    >  1:107138/1‑149 (MQ=255)
                                                                                cTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt                                                                   >  3:649964/1‑149 (MQ=255)
                                                                                 ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATc                                                                  >  1:73794/1‑149 (MQ=255)
                                                                                    cccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgc                                                               <  4:649964/149‑1 (MQ=255)
                                                                                    cccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgc                                                               <  3:320907/149‑1 (MQ=255)
                                                                                     ccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgcg                                                              >  1:493209/1‑149 (MQ=255)
                                                                                         aaTGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCt                                                          <  2:279121/149‑1 (MQ=255)
                                                                                          aTGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTg                                                         <  3:318748/149‑1 (MQ=255)
                                                                                              gAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCggg                                                                                                                                              <  1:256300/60‑1 (MQ=255)
                                                                                              gAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCggg                                                                                                                                              >  2:256300/1‑60 (MQ=255)
                                                                                                    tttCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCa                                                                 >  1:466656/1‑131 (MQ=255)
                                                                                                    tttCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCa                                                                 <  2:466656/131‑1 (MQ=255)
                                                                                                    tttCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAat                                               >  4:256030/1‑149 (MQ=255)
                                                                                                        gtcgtGTGAAACAAAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAattatt                                           <  2:618623/149‑1 (MQ=255)
                                                                                                           gtgtGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGc                                        <  2:3046/149‑1 (MQ=255)
                                                                                                                aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTAc                                   >  1:286505/1‑149 (MQ=255)
                                                                                                                aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTAc                                   >  1:10931/1‑149 (MQ=255)
                                                                                                                    aGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGc                               <  3:325792/149‑1 (MQ=255)
                                                                                                                    aGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGc                               >  1:445781/1‑149 (MQ=255)
                                                                                                                        tGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGaa                           <  4:206336/149‑1 (MQ=255)
                                                                                                                         gCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAAc                          <  2:338008/149‑1 (MQ=255)
                                                                                                                           gAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACgg                        >  3:304244/1‑149 (MQ=255)
                                                                                                                            aaGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGgc                       <  3:256030/149‑1 (MQ=255)
                                                                                                                             aGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAAAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGgcg                      <  4:23594/149‑1 (MQ=255)
                                                                                                                                cGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGa                   >  4:461846/1‑149 (MQ=255)
                                                                                                                                 gAACAATAAGGCCTCTCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGaa                  >  2:598317/1‑149 (MQ=255)
                                                                                                                                 gAACAATAAGGCCTCCCGAATCGGGGGGCCTTTTTTATTGAGAACAAAAGGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGACAGAAAATCAGCGCGCTGGATGTAAAGTTATTAGCGTTACTGGCAGAACGGCGCGaa                  >  4:468327/1‑149 (MQ=255)
                                                                                                                                  ataaaaaaGGCCCCCCGATTTGGGAGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAAc                 >  2:21836/7‑149 (MQ=255)
                                                                                                                                   aCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAAttatta                                          >  2:493818/1‑123 (MQ=255)
                                                                                                                                   aCAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAAttatta                                          <  1:493818/123‑1 (MQ=255)
                                                                                                                                     aaaaaGGCCCCCCGATTTGGGAGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTgg              >  1:619226/4‑149 (MQ=255)
                                                                                                                                      aaaaGGCCCCCCGATTTGGGAGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGc             >  1:623100/3‑149 (MQ=255)
                                                                                                                                            ccTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGa       >  1:7629/1‑149 (MQ=255)
                                                                                                                                              tCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAgg     <  2:493209/149‑1 (MQ=255)
                                                                                                                                              tCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAgg     >  4:592821/1‑149 (MQ=255)
                                                                                                                                                ccAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTg   <  4:170045/149‑1 (MQ=255)
                                                                                                                                                ccAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTg   >  1:551645/1‑149 (MQ=255)
                                                                                                                                                 cAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTgg  <  3:723384/149‑1 (MQ=255)
                                                                                                                                                 cAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTgg  <  1:598317/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACAGTACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGG  >  NC_000913/2737557‑2737850

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