Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,796,968 G→A 7.9% H970H (CAC→CAT gltB ← L‑glutamate synthase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,796,9680GA 7.9% 450.1 / 20.9 152H970H (CAC→CATgltBL‑glutamate synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (69/71);  new base A (6/6);  total (75/77)
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.60e-01

Only 100 of 158 total aligned reads displayed.

TCAGGTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTGACCAG  >  NC_002947/5796820‑5797117
                                                                                                                                                    |                                                                                                                                                      
tCAGGTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGgt                                                                                                                                                       <  2:4064599/150‑1 (MQ=255)
  aGGTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGgtgt                                                                                                                                                     >  2:4028395/1‑150 (MQ=255)
   ggTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGtgtg                                                                                                                                                    <  1:75979/150‑1 (MQ=255)
    gTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTgtgg                                                                                                                                                   >  2:3109369/1‑150 (MQ=255)
    gTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTgtgg                                                                                                                                                   <  1:3369255/150‑1 (MQ=255)
    gTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTgtgg                                                                                                                                                   <  1:4393031/150‑1 (MQ=255)
    gTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTgtgg                                                                                                                                                   <  1:160496/150‑1 (MQ=255)
      cgaagcgggCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTgtggcg                                                                                                                                                 >  1:4604460/8‑150 (MQ=255)
         cgcgtgCCTTGGCAACACCGGCAGCAACCGTACCAACGCCCACTTCCGACACCAGCTTGACCGAAATATCAGCCGCCGGATTGACGTTCTTCAGATCGTAGATCAGCTGTGCCAGATCTTCGATCGAGTAGATGTCATGATgtggcggag                                                                                                                                              <  1:5288063/145‑3 (MQ=255)
           tAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTgtggcggtggc                                                                                                                                            <  1:5202112/150‑1 (MQ=255)
             ggCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCga                                                                                                                                          <  1:1489253/150‑1 (MQ=255)
                cTTGGCAACACCGGCAGCAACCGTACCAACGCCCACTTCCGACACCAGCTTGACCGAAATATCAGCCGCCGGATTGACGTTCTTCAGATCGTAGATCAGCTGTGCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAAt                                                                                                                                       >  2:633961/1‑150 (MQ=255)
                 ttGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATc                                                                                                                                      >  1:2519307/1‑150 (MQ=255)
                  tGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCa                                                                                                                                     <  1:151347/150‑1 (MQ=255)
                  tGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCa                                                                                                                                     <  2:1340796/150‑1 (MQ=255)
                    gCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAgg                                                                                                                                   <  2:3878265/150‑1 (MQ=255)
                    gCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAgg                                                                                                                                   <  2:5193173/150‑1 (MQ=255)
                     ccACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCCGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCgaga                                                                                                                                        <  1:1157217/144‑1 (MQ=255)
                     ccACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAggg                                                                                                                                  <  1:2428121/150‑1 (MQ=255)
                      cacaCCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGt                                                                                                                                 >  2:5433442/1‑150 (MQ=255)
                      cacaCCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGt                                                                                                                                 >  2:3572625/1‑150 (MQ=255)
                       acacCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTc                                                                                                                                <  1:705382/150‑1 (MQ=255)
                       acacCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTc                                                                                                                                >  2:618053/1‑150 (MQ=255)
                         acCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCAc                                                                                                                              >  2:3550137/1‑150 (MQ=255)
                         acCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCAc                                                                                                                              >  2:5381923/1‑150 (MQ=255)
                           cTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGc                                                                                                                            <  1:3339406/150‑1 (MQ=255)
                                      gtcgatgcGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCa                                                                                                                 >  2:3601521/8‑150 (MQ=255)
                                      gTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCa                                                                                                                 >  1:1757941/1‑150 (MQ=255)
                                        gCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATa                                                                                                               >  2:4516221/1‑150 (MQ=255)
                                              gCCCACTTCCGACACCAGCTTGACCGAAATATCAGCCGCCGGATTGACGTTCTTCAGATCGTAGATCAGCTGTGCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCgaaa                                                                                                                                        <  1:4792164/119‑3 (MQ=255)
                                                 cGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCtt                                                                                                      <  1:3728950/150‑1 (MQ=255)
                                                 cGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCtt                                                                                                      <  1:3728974/150‑1 (MQ=255)
                                                  gCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTg                                                                                                     <  2:3465324/150‑1 (MQ=255)
                                                         cTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCa                                                                                              <  1:5433442/150‑1 (MQ=255)
                                                         cTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCa                                                                                              <  1:218372/150‑1 (MQ=255)
                                                          tACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAg                                                                                             >  2:4644416/1‑150 (MQ=255)
                                                           aCCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAgg                                                                                            >  1:1952749/1‑150 (MQ=255)
                                                              aGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCg                                                                                         >  1:2059020/1‑150 (MQ=255)
                                                              aGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCg                                                                                         >  1:2686780/1‑150 (MQ=255)
                                                                  tGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGa                                                                                     <  1:4995614/150‑1 (MQ=255)
                                                                   gACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGAc                                                                                    >  1:4008145/1‑150 (MQ=255)
                                                                   gACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGCGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGAc                                                                                    <  1:2610047/150‑1 (MQ=255)
                                                                    aCCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGAcc                                                                                   >  2:2823789/1‑150 (MQ=255)
                                                                     ccGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAgg                                                                                                                                   >  1:1804737/1‑101 (MQ=255)
                                                                     ccGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCt                                                                                 >  1:1293791/1‑150 (MQ=255)
                                                                      cGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCtt                                                                                >  1:686457/1‑150 (MQ=255)
                                                                       gATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCttt                                                                               <  2:1193326/150‑1 (MQ=255)
                                                                            cAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGcc                                                                          >  1:3617418/1‑150 (MQ=255)
                                                                            cAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGcc                                                                          >  2:4982750/1‑150 (MQ=255)
                                                                            cAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGcc                                                                          >  1:3512293/1‑150 (MQ=255)
                                                                             aGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGccc                                                                         >  1:4182468/1‑150 (MQ=255)
                                                                              gcgcCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGACCGAGTAAATGTCATGgtg                                                                                                                                                      >  1:2155997/1‑73 (MQ=255)
                                                                              gcgcCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGACCGAGTAAATGTCATGgtg                                                                                                                                                      <  2:2155997/73‑1 (MQ=255)
                                                                                  cTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCaa                                                                    >  1:2368601/1‑150 (MQ=255)
                                                                                       ggTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGgc                                                               >  1:2636256/1‑150 (MQ=255)
                                                                                        gTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGgcc                                                              >  2:5243330/1‑150 (MQ=255)
                                                                                            aCCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTTCCGCCCGGCAACTGgccgcc                                                           <  2:3869031/150‑1 (MQ=255)
                                                                                            aCCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGccc                                                          >  1:3129880/1‑150 (MQ=255)
                                                                                             ccTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTgg                                                                <  2:4505897/143‑1 (MQ=255)
                                                                                             ccTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTgg                                                                >  1:4505897/1‑143 (MQ=255)
                                                                                               tGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCa                                                       <  2:2686780/150‑1 (MQ=255)
                                                                                                gCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCAc                                                      >  2:5019105/1‑150 (MQ=255)
                                                                                                gCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCAc                                                      >  2:331521/1‑150 (MQ=255)
                                                                                                gCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCAc                                                      >  2:331492/1‑150 (MQ=255)
                                                                                                 cTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCAcc                                                     <  2:2275097/150‑1 (MQ=255)
                                                                                                  ttGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCAccc                                                    <  1:1649881/150‑1 (MQ=255)
                                                                                                    gAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCgg                                                  >  2:4157685/1‑150 (MQ=255)
                                                                                                    cagatCGTAGATCAGCTGTGCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAATGAGGCCAACGCCCGGCGTCGAGTGACGGGTCTTGGCAATCGTCGCATCAACCTTATGGCCCGGCAGCTGACCGCCTTCACCCgg                                                  <  1:5443385/146‑1 (MQ=255)
                                                                                                       gTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCtt                                               >  2:5234990/1‑150 (MQ=255)
                                                                                                       atCGTAGATCAGCTGTGCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAATGAGGCCAACGCCCGGCGTCGAGTGACGGGTCTTGGCAATCGTCGCATCAACCTTATGGCCCGGCAGCTGACCGCCTTCACCCGGCtt                                               >  2:3557207/2‑150 (MQ=255)
                                                                                                        tCGTAGATCAGCTGTGCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAATGAGGCCAACGCCCGGCGTCGAGTGACGGGTCTTGGCAATCGTCGCATCAACCTTATGGCCCGGCAGCTGACCGCCTTCACCCGGCttt                                              >  2:1587054/1‑149 (MQ=255)
                                                                                                         cGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTgg                                             <  1:1987592/150‑1 (MQ=255)
                                                                                                          gTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGgc                                            >  2:5076516/1‑150 (MQ=255)
                                                                                                               tCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCt                                       <  1:3579875/150‑1 (MQ=255)
                                                                                                               tCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCt                                       <  1:3601521/150‑1 (MQ=255)
                                                                                                               tCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCt                                       >  1:517017/1‑150 (MQ=255)
                                                                                                               tCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCt                                       <  1:1628098/150‑1 (MQ=255)
                                                                                                                 aGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTgg                                     <  2:1746117/150‑1 (MQ=255)
                                                                                                                 aGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTgg                                     <  2:1006554/150‑1 (MQ=255)
                                                                                                                 aGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTgg                                     <  2:1792418/150‑1 (MQ=255)
                                                                                                                   cTGTGCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAATGAGGCCAACGCCCGGCGTCGAGTGACGGGTCTTGGCAATCGTCGCATCAACCTTATGGCCCGGCAGCTGACCGCCTTCACCCGGCTTTGCACCCTGAGc                                   >  1:1095749/1‑150 (MQ=255)
                                                                                                                    tgtgCCAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAATGAGGCCAACGCCCGGCGTCGAGTGACGGGTCTTGGCAATCGTCGCATCAACCTTATGGCCCGGCAGCTGACCGCCTTCACCCGGCTTTGCACCCTGAGCa                                  >  2:1805805/4‑150 (MQ=255)
                                                                                                                        ccAGATCTTCGATCGAGTAGATGTCATGATGTGGCGGAGGCGAAATGAGGCCAACGCCCGGCGTCGAGTGACGGGTCTTGGCAATCGTCGCATCAACCTTATGGCCCGGCAGCTGACCGCCTTCACCCGGCTTTGCACCCTGAGCAACCt                              >  2:1976558/1‑150 (MQ=255)
                                                                                                                           gATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGa                           <  1:3772774/150‑1 (MQ=255)
                                                                                                                            aTCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGAt                          <  2:5283507/150‑1 (MQ=255)
                                                                                                                               ttCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTg                       >  2:2615415/1‑150 (MQ=255)
                                                                                                                               ttCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTg                       >  2:2075813/1‑150 (MQ=255)
                                                                                                                               ttCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTg                       >  1:4525609/1‑150 (MQ=255)
                                                                                                                               ttCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTg                       <  1:2886802/150‑1 (MQ=255)
                                                                                                                                 cgatcgaGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTgca                     <  1:108379/150‑1 (MQ=255)
                                                                                                                                 cgatcgaGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTgca                     <  1:5383121/150‑1 (MQ=255)
                                                                                                                                 cgatcgaGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTgca                     <  1:5383066/150‑1 (MQ=255)
                                                                                                                                 cgatcgaGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTgca                     <  1:2804223/150‑1 (MQ=255)
                                                                                                                                     cgaGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCAc                 <  1:1017910/150‑1 (MQ=255)
                                                                                                                                       aGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACtt               >  1:3506285/1‑150 (MQ=255)
                                                                                                                                           aaTGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGc           >  1:1329813/1‑150 (MQ=255)
                                                                                                                                             tGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGt         >  1:4036520/1‑150 (MQ=255)
                                                                                                                                               tCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTg       >  2:101160/1‑150 (MQ=255)
                                                                                                                                                 aTGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTGAc     <  1:2550934/150‑1 (MQ=255)
                                                                                                                                                    gtgtGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCACATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTGACCAg  >  2:713790/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TCAGGTCGGCATAGGCCTTGGCCACACCTGCGGCGATGGTGCCCACGCCCGCTTCGGCTACCAGCTTGACCGATACCAGCGCCTGCGGGTTGACCTGCTTGAGGTCGTAGATCAGCTGGGCCAGATCTTCGATCGAGTAAATGTCATGGTGTGGCGGTGGCGAGATCAGGGTCACGCCCGGTACCGCATAACGCAGCTTGGCGATCAGGCCGTTGACCTTTCCGCCCGGCAACTGGCCGCCCTCACCCGGCTTGGCGCCCTGGGCAACCTTGATCTGCAGCACTTCGGCGTTGACCAG  >  NC_002947/5796820‑5797117

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