Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,760 T→A 6.5% intergenic (+59/+35) guaA → / ← PP_1033 glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,7600TA 6.5% 487.7 / 16.2 151intergenic (+59/+35)guaA/PP_1033glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (75/64);  new base A (3/7);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

AGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGC  >  NC_002947/1178612‑1178909
                                                                                                                                                    |                                                                                                                                                     
aGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTTGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGAtt                                                                                                                                                      <  2:2351521/150‑1 (MQ=255)
  cGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGc                                                                                                                                                    >  1:382798/1‑150 (MQ=255)
    gTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCgg                                                                                                                                                  <  2:401479/150‑1 (MQ=255)
    gTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCgg                                                                                                                                                  <  1:4116811/150‑1 (MQ=255)
       gTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCt                                                                                                                                               <  1:1661833/150‑1 (MQ=255)
          tcatcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTc                                                                                                                                            >  1:4315865/1‑150 (MQ=255)
           catcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                           <  2:772353/150‑1 (MQ=255)
           catcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                           >  2:2479734/1‑150 (MQ=255)
           catcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                           <  1:804204/150‑1 (MQ=255)
            atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                          >  1:3078107/1‑150 (MQ=255)
            atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                          <  1:1053723/150‑1 (MQ=255)
            atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                          >  1:808705/1‑150 (MQ=255)
            atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                          <  1:2857794/150‑1 (MQ=255)
            atcaACGAAATCGACGGTATCTCACGCGTTACCGACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                          <  2:418704/150‑1 (MQ=255)
             tcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTgg                                                                                                                                         >  1:1949961/1‑150 (MQ=255)
               aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGtt                                                                                                                                       >  2:3408138/1‑150 (MQ=255)
               aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGtt                                                                                                                                       <  1:1284619/150‑1 (MQ=255)
                 cGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCt                                                                                                                                     <  1:2165286/150‑1 (MQ=255)
                 cGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTGGtgta                                                                                                                                     <  1:837699/150‑4 (MQ=39)
                   aaaTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTg                                                                                                                                   <  2:1353804/150‑1 (MQ=255)
                   aaaTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTg                                                                                                                                   <  2:1672783/150‑1 (MQ=255)
                      tCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttc                                                                                                                                >  2:437758/1‑150 (MQ=255)
                             tATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGa                                                                                                                         <  2:1960868/150‑1 (MQ=255)
                                 tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                     >  2:2752946/1‑150 (MQ=255)
                                 tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                     >  2:3094211/1‑150 (MQ=255)
                                 tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                     >  1:3807097/1‑150 (MQ=255)
                                 tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                     <  1:1356731/150‑1 (MQ=255)
                                 tcGCGCGGTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                     <  1:2291426/150‑1 (MQ=255)
                                    cgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                                  <  2:4117629/150‑1 (MQ=255)
                                     gcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAt                                                                                                                 >  1:4091281/1‑150 (MQ=255)
                                      cgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAtt                                                                                                                <  2:3905700/150‑1 (MQ=255)
                                          aCCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGGGGGCTTTCTGGTTCTTGTTCGTTGAGATTGGTTAttcttg                                                                                                            >  2:4052150/1‑149 (MQ=255)
                                           ccTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGCCGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                           <  1:3330515/150‑1 (MQ=255)
                                            cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTg                                                                                                          >  2:3235921/1‑150 (MQ=255)
                                             tACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTgg                                                                                                         >  2:1884303/1‑150 (MQ=255)
                                              acgacgTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGc                                                                                                        <  2:1191057/150‑1 (MQ=255)
                                                   gtgtCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                                   >  1:3830298/1‑150 (MQ=255)
                                                         aGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                             <  1:2384752/150‑1 (MQ=255)
                                                             aGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAgc                                                                                         >  1:674534/1‑150 (MQ=255)
                                                               ccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAgcgc                                                                                       >  1:471483/1‑150 (MQ=255)
                                                                  ccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCttt                                                                                    <  2:3508743/150‑1 (MQ=255)
                                                                   cggCCACGATCGAGTGGGAGTGATCGGGTAGCCCGCCCAGGCTGGCACCCGCAGGACGGATACACCCGACAGCCCGCATGATTGCGGGCTTTCTGGTTCTGGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAccgcgtac                                                                                   <  2:2717085/148‑5 (MQ=255)
                                                                     gCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTg                                                                                 <  2:632503/150‑1 (MQ=255)
                                                                        aCTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTa                                                                              >  1:413464/1‑150 (MQ=255)
                                                                          tATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTTGTAgg                                                                            >  2:2531201/1‑150 (MQ=255)
                                                                          tATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAgg                                                                            >  2:3829182/1‑150 (MQ=255)
                                                                          tATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAgg                                                                            <  2:85744/150‑1 (MQ=255)
                                                                              gAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGtt                                                                        >  1:1762485/1‑150 (MQ=255)
                                                                              gAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGtt                                                                        >  1:3265656/1‑150 (MQ=255)
                                                                                gTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtat                                                                      <  1:3715530/150‑1 (MQ=255)
                                                                                 tGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                                     <  2:1883196/150‑1 (MQ=255)
                                                                                 tGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                                     <  2:1811739/150‑1 (MQ=255)
                                                                                 tGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                                     >  1:1078502/1‑150 (MQ=255)
                                                                                  gggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAg                                                                    <  1:4012151/150‑1 (MQ=255)
                                                                                   ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAgg                                                                   <  1:2752946/150‑1 (MQ=255)
                                                                                    gAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTGTCTGGTAGGCGTTATAGGa                                                                  <  1:2780444/150‑1 (MQ=255)
                                                                                      gTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                                <  2:3757760/150‑1 (MQ=255)
                                                                                      gTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                                <  2:1013181/150‑1 (MQ=255)
                                                                                       tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAag                                                               <  2:602970/150‑1 (MQ=255)
                                                                                         aTCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAg                                                                    <  1:926353/143‑1 (MQ=255)
                                                                                          tCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagcag                                                            >  1:849191/1‑150 (MQ=255)
                                                                                            gggTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGcc                                                          >  1:2827442/1‑150 (MQ=255)
                                                                                             ggTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCt                                                         <  2:1627622/150‑1 (MQ=255)
                                                                                             ggTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCt                                                         >  2:3737478/1‑150 (MQ=255)
                                                                                                agccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagcagc                                                           >  2:1958309/1‑145 (MQ=255)
                                                                                                agccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagcagc                                                           <  1:1958309/145‑1 (MQ=255)
                                                                                                 gccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                                     >  1:3364598/1‑150 (MQ=255)
                                                                                                  ccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctcc                                                    >  1:523332/1‑150 (MQ=255)
                                                                                                     gccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCa                                                 <  2:3172487/150‑1 (MQ=255)
                                                                                                     gccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCa                                                 >  1:3878150/1‑150 (MQ=255)
                                                                                                        cAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTg                                              >  2:411252/1‑150 (MQ=255)
                                                                                                         aGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGc                                             >  1:1384504/1‑150 (MQ=255)
                                                                                                             tGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATg                                         >  1:2735754/1‑150 (MQ=255)
                                                                                                                   ccGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCg                                   <  2:4039545/150‑1 (MQ=255)
                                                                                                                      cAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCa                                <  2:668373/150‑1 (MQ=255)
                                                                                                                      cAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTGTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCa                                <  1:2567789/150‑1 (MQ=255)
                                                                                                                                acacCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgc                      <  2:1949961/150‑1 (MQ=255)
                                                                                                                                  acCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGt                    >  1:599003/1‑150 (MQ=255)
                                                                                                                                   ccAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGtt                   >  2:3127940/1‑150 (MQ=255)
                                                                                                                                    cAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTg                  >  1:1428643/1‑150 (MQ=255)
                                                                                                                                    cAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTg                  >  2:971571/1‑150 (MQ=255)
                                                                                                                                    cAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTg                  <  2:1492553/150‑1 (MQ=255)
                                                                                                                                    cAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTg                  <  2:1492435/150‑1 (MQ=255)
                                                                                                                                     aGACAGCCCGCATGAGTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTGTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                 <  1:3656175/150‑1 (MQ=255)
                                                                                                                                     aGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                 <  2:1078502/150‑1 (MQ=255)
                                                                                                                                     aGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                 >  1:1596347/1‑150 (MQ=255)
                                                                                                                                     aGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg                 >  2:3761862/1‑150 (MQ=255)
                                                                                                                                      gAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg                >  1:1344723/1‑150 (MQ=255)
                                                                                                                                         aGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGt             >  2:3403128/1‑150 (MQ=255)
                                                                                                                                           cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc           >  2:369/1‑150 (MQ=255)
                                                                                                                                           cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc           >  2:353/1‑150 (MQ=255)
                                                                                                                                           cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTga           <  2:2202982/150‑2 (MQ=255)
                                                                                                                                           cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTga           <  1:3829182/150‑2 (MQ=255)
                                                                                                                                           cccGCATGAGTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAGTCTTCGGGCTTTCTGGCAAAGAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCCCAACGCGTGGAGGGTga           <  2:3029034/150‑2 (MQ=255)
                                                                                                                                            ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCt          >  1:475906/1‑150 (MQ=255)
                                                                                                                                             cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtt         >  1:4315651/1‑150 (MQ=255)
                                                                                                                                                aTGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTc      <  2:4315865/150‑1 (MQ=255)
                                                                                                                                                aTGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTc      <  1:1291773/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGc  <  2:471483/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGc  >  1:2250577/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTCCTGC  >  NC_002947/1178612‑1178909

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: