Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,629 G→A 8.6% Q96* (CAG→TAG) ‡ amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,6290GA 8.6% 681.8 / 38.8 233Q96* (CAG→TAG) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (126/87);  new base A (9/11);  total (135/98)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 242 total aligned reads displayed.

ACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGC  >  NC_002947/6005481‑6005778
                                                                                                                                                    |                                                                                                                                                     
aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                      <  2:4869976/150‑1 (MQ=255)
      aCCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATc                                                                                                                                                <  2:1264043/150‑1 (MQ=255)
       ccACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt                                                                                                                                               >  2:2761846/1‑150 (MQ=255)
        cacaCCCAGCGGGTGCCAGGGTTCACGCAGGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACCGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTGCCGGGGTGATCtt                                                                                                                                              >  2:5219754/1‑150 (MQ=255)
           acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                                                                                                                                           >  1:4498888/1‑150 (MQ=255)
            cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                                          <  2:701939/150‑1 (MQ=255)
            cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                                          <  2:5073871/150‑1 (MQ=255)
              cAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGc                                                                                                                                        >  1:3285845/1‑150 (MQ=255)
              cAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGc                                                                                                                                        >  1:3233702/1‑150 (MQ=255)
                 cGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtc                                                                                                                                     <  2:2490343/150‑1 (MQ=255)
                  gggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcg                                                                                                                                    <  1:5053576/150‑1 (MQ=255)
                   ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                                                                                                                                   <  2:2954615/150‑1 (MQ=255)
                   ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                                                                                                                                   <  1:2483978/150‑1 (MQ=255)
                     tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                                                                                                                                 <  2:964454/150‑1 (MQ=255)
                       ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                               <  1:2261960/150‑1 (MQ=255)
                       ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                               >  1:192613/1‑150 (MQ=255)
                        cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCCACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                              >  1:4803847/1‑150 (MQ=255)
                        cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                              >  2:5415669/1‑150 (MQ=255)
                        cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                              >  2:4060492/1‑150 (MQ=255)
                        cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                              >  2:186214/1‑150 (MQ=255)
                          ggTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGAc                                                                                                                            <  2:5309993/150‑1 (MQ=255)
                            tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCa                                                                                                                          >  2:1766675/1‑150 (MQ=255)
                            tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCa                                                                                                                          >  2:2045363/1‑150 (MQ=255)
                                  gCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGc                                                                                                                    >  2:4308316/1‑150 (MQ=255)
                                        ggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                                                                                                                                  <  2:4098225/130‑1 (MQ=255)
                                        ggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                                                                                                                                  >  1:4098225/1‑130 (MQ=255)
                                         gtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTc                                                                                                             <  2:4220529/150‑1 (MQ=255)
                                          tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGa                                                                                                                                                  >  1:3987897/1‑112 (MQ=255)
                                          tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGa                                                                                                                                                  <  2:3987897/112‑1 (MQ=255)
                                          tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCg                                                                                                            >  1:1043364/1‑150 (MQ=255)
                                               cGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCtt                                                                                                       <  2:16739/150‑1 (MQ=255)
                                               agcgcgtTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCtt                                                                                                       <  1:2447977/147‑1 (MQ=255)
                                                gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTGCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtg                                                                                                      <  2:4091418/150‑1 (MQ=255)
                                                 ggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgt                                                                                                     <  2:4392960/150‑1 (MQ=255)
                                                   cgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGc                                                                                                   >  2:2082903/1‑150 (MQ=255)
                                                         gAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgc                                                                                             >  1:5359240/1‑150 (MQ=255)
                                                                  gTCAGGACGTACAGCGGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCg                                                                                    <  1:1570652/150‑1 (MQ=255)
                                                                   tCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                                   >  2:2532719/1‑150 (MQ=255)
                                                                   tCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                                   >  1:4902594/1‑150 (MQ=255)
                                                                   tCAGGCCGTACAGCTGGCGCGACAGGCAAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCGTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                                   <  1:89542/150‑1 (MQ=255)
                                                                       gCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGc                                                                                                   <  2:1522837/130‑1 (MQ=255)
                                                                        ccGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGAGGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                              <  2:2363999/150‑1 (MQ=255)
                                                                         cGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACa                                                                             <  1:3720556/150‑1 (MQ=255)
                                                                          gTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAg                                                                            <  1:3030600/150‑1 (MQ=255)
                                                                              aGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCgg                                                                                                                                         <  1:4213279/85‑1 (MQ=255)
                                                                              aGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGgcgc                                                                        >  2:3171207/1‑150 (MQ=255)
                                                                               gCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGgcgga                                                                       <  2:1925798/150‑3 (MQ=255)
                                                                                cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                              <  1:2058323/142‑1 (MQ=255)
                                                                                cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                      >  2:4581772/1‑150 (MQ=255)
                                                                                  ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                    <  2:4254694/150‑1 (MQ=255)
                                                                                  ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                    <  2:4830637/150‑1 (MQ=255)
                                                                                     gcgACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGtt                                                                 >  1:5267235/1‑150 (MQ=255)
                                                                                     gcgACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGtt                                                                 <  1:4264560/150‑1 (MQ=255)
                                                                                     gcgACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGtt                                                                 >  1:4094927/1‑150 (MQ=255)
                                                                                         cAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGcc                                                             <  2:1373024/150‑1 (MQ=255)
                                                                                         cAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGcc                                                             >  1:1459831/1‑150 (MQ=255)
                                                                                          aGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGt                                                                                                              <  2:3873735/100‑1 (MQ=255)
                                                                                          aGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGt                                                                                                              >  1:3873735/1‑100 (MQ=255)
                                                                                          aGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCa                                                            >  2:3470227/1‑150 (MQ=255)
                                                                                            gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACg                                                          >  2:1635546/1‑150 (MQ=255)
                                                                                            gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACg                                                          >  1:2636669/1‑150 (MQ=255)
                                                                                             ccAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgc                                                         <  1:1413702/150‑2 (MQ=255)
                                                                                              cAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgc                                                        >  1:3624850/1‑150 (MQ=255)
                                                                                              cAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgc                                                        >  2:1033455/1‑150 (MQ=255)
                                                                                              cAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgc                                                        >  2:4681210/1‑150 (MQ=255)
                                                                                               aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcg                                                       <  2:89871/150‑1 (MQ=255)
                                                                                                aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcgc                                                      <  1:622624/150‑1 (MQ=255)
                                                                                                aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcgc                                                      <  2:3641628/150‑1 (MQ=255)
                                                                                                 cGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGcgcg                                                     <  1:1734226/150‑1 (MQ=255)
                                                                                                     gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                 >  2:4518905/1‑150 (MQ=255)
                                                                                                     gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                 <  1:4382082/150‑1 (MQ=255)
                                                                                                     gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                 >  1:2552975/1‑150 (MQ=255)
                                                                                                         tCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGc                                             >  1:2813146/1‑150 (MQ=255)
                                                                                                             aGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCacac                                         >  2:644819/1‑150 (MQ=255)
                                                                                                              gATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACt                                        >  1:3010747/1‑150 (MQ=255)
                                                                                                               aTGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGc                                                              >  2:179533/1‑127 (MQ=255)
                                                                                                               aTGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGc                                                              <  1:179533/127‑1 (MQ=255)
                                                                                                               aTGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTg                                       >  1:652297/1‑150 (MQ=255)
                                                                                                                         tttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCtt                             >  1:3804364/1‑150 (MQ=255)
                                                                                                                         tttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCtt                             >  1:3447892/1‑150 (MQ=255)
                                                                                                                          ttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTc                            >  2:1868899/1‑150 (MQ=255)
                                                                                                                          ttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTc                            >  2:2282716/1‑150 (MQ=255)
                                                                                                                            ccTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGa                          >  1:31152/1‑150 (MQ=255)
                                                                                                                             cTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGaa                         >  1:3898642/1‑150 (MQ=255)
                                                                                                                             cTGCACGTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGaa                         <  1:4955194/150‑1 (MQ=255)
                                                                                                                              tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg                        >  2:4971390/1‑150 (MQ=255)
                                                                                                                              tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg                        >  2:2792212/1‑150 (MQ=255)
                                                                                                                               gCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAgg                       >  1:381275/1‑150 (MQ=255)
                                                                                                                               gCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAgg                       <  2:2461401/150‑1 (MQ=255)
                                                                                                                                  cTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcgc                    >  1:4144819/1‑150 (MQ=255)
                                                                                                                                   ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                   >  2:4061377/1‑150 (MQ=255)
                                                                                                                                   ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                   >  2:3334061/1‑150 (MQ=255)
                                                                                                                                   ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                   >  2:586901/1‑150 (MQ=255)
                                                                                                                                       cccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTgg               <  1:4060492/150‑1 (MQ=255)
                                                                                                                                         cAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGc             <  2:1043364/150‑1 (MQ=255)
                                                                                                                                            gCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTg          >  2:3110413/1‑150 (MQ=255)
                                                                                                                                                   tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGTACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg   >  2:3822892/1‑150 (MQ=255)
                                                                                                                                                   tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg   >  2:3657603/1‑150 (MQ=255)
                                                                                                                                                   tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg   >  1:4356216/1‑150 (MQ=255)
                                                                                                                                                    gggTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGc  <  1:2761846/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGC  >  NC_002947/6005481‑6005778

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: