Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,628 T→G 8.7% Q96P (CAG→CCG) ‡ 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,6280TG 8.7% 740.1 / 45.5 237Q96P (CAG→CCG) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (124/92);  new base G (10/11);  total (134/103)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 242 total aligned reads displayed.

TCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATG  >  NC_002947/6005479‑6005777
                                                                                                                                                     |                                                                                                                                                     
  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  1:528525/150‑1 (MQ=255)
  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  2:2208579/150‑1 (MQ=255)
  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     >  1:4866343/1‑150 (MQ=255)
  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     >  1:4722022/1‑150 (MQ=255)
  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  1:2796757/150‑1 (MQ=255)
  aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  2:4869976/150‑1 (MQ=255)
     ccGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTg                                                                                                                                                  <  1:3866624/150‑1 (MQ=255)
       gACCACACCCAGCGGGTGCCAGGTTTCACGCAAGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACGTCACCCAGGAAGGCCTGGGTGaa                                                                                                                                                <  2:125621/150‑2 (MQ=255)
         ccACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCt                                                                                                                                              >  2:2761846/1‑150 (MQ=255)
          cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                                                                                                                                             <  2:1387548/150‑1 (MQ=255)
          cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                                                                                                                                             <  1:418855/150‑1 (MQ=255)
          cacaCCCAGCGGGTGCCAGGGTTCACGCAGGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACCGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTGCCGGGGTGATCtt                                                                                                                                             >  2:5219754/1‑150 (MQ=255)
              cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                                         <  2:5073871/150‑1 (MQ=255)
              cccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCg                                                                                                                                         <  2:701939/150‑1 (MQ=255)
                   cGGATGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtc                                                                                                                                    <  2:2490343/150‑1 (MQ=255)
                     ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                                                                                                                                  <  1:2483978/150‑1 (MQ=255)
                     ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                                                                                                                                  <  2:2954615/150‑1 (MQ=255)
                       tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                                                                                                                                <  2:964454/150‑1 (MQ=255)
                         ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                              >  1:192613/1‑150 (MQ=255)
                         ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                              >  1:3770085/1‑150 (MQ=255)
                         ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                              <  1:2261960/150‑1 (MQ=255)
                          cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                              >  1:3160155/1‑149 (MQ=255)
                          cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                             >  2:5415669/1‑150 (MQ=255)
                          cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                             >  2:4060492/1‑150 (MQ=255)
                            ggTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGAc                                                                                                                           <  2:5309993/150‑1 (MQ=255)
                              tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCa                                                                                                                         >  2:1766730/1‑150 (MQ=255)
                                    gCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGc                                                                                                                   >  2:4308316/1‑150 (MQ=255)
                                        gtggtgGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAg                                                                                                               >  1:3926149/1‑150 (MQ=255)
                                          ggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                                                                                                                                 <  2:4098225/130‑1 (MQ=255)
                                          ggtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTAGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgat                                                                                                                                 >  1:4098194/1‑130 (MQ=255)
                                           gtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTc                                                                                                            <  1:631784/150‑1 (MQ=255)
                                           gtggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTc                                                                                                            <  2:4220529/150‑1 (MQ=255)
                                            tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGa                                                                                                                                                 <  2:3987897/112‑1 (MQ=255)
                                            tggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCg                                                                                                           >  1:1043364/1‑150 (MQ=255)
                                              gCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGc                                                                                                         <  1:4423869/150‑1 (MQ=255)
                                              gCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGAGATCTTGCCGGCTtcgatc                                                                                                                                >  1:5203614/1‑127 (MQ=255)
                                                ccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCt                                                                                                       >  1:2941850/1‑150 (MQ=255)
                                                 agcgcgtTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCtt                                                                                                      <  1:2447977/147‑1 (MQ=255)
                                                  gggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTGCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtg                                                                                                     <  2:4091418/150‑1 (MQ=255)
                                                   ggCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgt                                                                                                    <  2:4392960/150‑1 (MQ=255)
                                                          ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg                                                                                             >  2:1293873/1‑150 (MQ=255)
                                                          ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg                                                                                             <  1:5369190/150‑1 (MQ=255)
                                                            aaGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgca                                                                                           >  1:4099313/1‑150 (MQ=255)
                                                                   ggTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTc                                                                                    >  1:4907226/1‑150 (MQ=255)
                                                                   ggTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTc                                                                                    <  1:3290011/150‑1 (MQ=255)
                                                                    gTCAGGACGTACAGCGGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCg                                                                                   <  1:1570652/150‑1 (MQ=255)
                                                                     tCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                                  >  2:2532719/1‑150 (MQ=255)
                                                                     tCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                                  <  1:89624/150‑1 (MQ=255)
                                                                          ccGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGAGGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                             <  2:2363999/150‑1 (MQ=255)
                                                                              aaaGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGgc                                                                         >  2:1525212/3‑150 (MQ=255)
                                                                                aGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGgcgc                                                                       >  2:3171207/1‑150 (MQ=255)
                                                                                 gCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGgcgga                                                                      <  2:1925637/150‑3 (MQ=255)
                                                                                  cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                             <  1:2058323/142‑1 (MQ=255)
                                                                                  cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                             >  2:2058323/1‑142 (MQ=255)
                                                                                  cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                     >  2:4581772/1‑150 (MQ=255)
                                                                                   tGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAcc                                                                    <  1:4041062/150‑1 (MQ=255)
                                                                                    ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                   >  2:1651153/1‑150 (MQ=255)
                                                                                    ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                   <  2:4830637/150‑1 (MQ=255)
                                                                                    ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                   <  2:4254694/150‑1 (MQ=255)
                                                                                           cAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGcc                                                            >  1:1459831/1‑150 (MQ=255)
                                                                                            agaccAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCa                                                           >  2:1034135/4‑150 (MQ=255)
                                                                                            aGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGt                                                                                                             <  2:3873735/100‑1 (MQ=255)
                                                                                            aGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCa                                                           >  2:3470227/1‑150 (MQ=255)
                                                                                              gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACg                                                         >  1:2636669/1‑150 (MQ=255)
                                                                                               ccAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgc                                                        <  1:1413702/150‑2 (MQ=255)
                                                                                                cAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgc                                                       >  2:4681210/1‑150 (MQ=255)
                                                                                                 aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcg                                                      <  2:89871/150‑1 (MQ=255)
                                                                                                  aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcgc                                                     <  2:3641628/150‑1 (MQ=255)
                                                                                                   ctgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttg                                                    <  1:3677564/148‑6 (MQ=255)
                                                                                                   cGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGcgcg                                                    <  1:1734226/150‑1 (MQ=255)
                                                                                                     gCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGa                                                  <  1:3548848/150‑1 (MQ=255)
                                                                                                       gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                >  1:2552975/1‑150 (MQ=255)
                                                                                                       gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                >  2:4518905/1‑150 (MQ=255)
                                                                                                           tCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGc                                            >  1:3267160/1‑150 (MQ=255)
                                                                                                               aGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCacac                                        >  2:644819/1‑150 (MQ=255)
                                                                                                               aGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCacac                                        >  1:3048747/1‑150 (MQ=255)
                                                                                                                 aTGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGc                                                             <  1:179533/127‑1 (MQ=255)
                                                                                                                  tGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgc                                     <  2:2293182/150‑1 (MQ=255)
                                                                                                                   gTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCtt                                                                                                                                     <  2:1856735/53‑1 (MQ=255)
                                                                                                                    tCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgcgc                                   >  1:3016469/1‑150 (MQ=255)
                                                                                                                        tCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAgg                               >  1:3882690/1‑150 (MQ=255)
                                                                                                                        tCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAgg                               >  1:4213946/1‑150 (MQ=255)
                                                                                                                            ttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTc                           >  2:1868688/1‑150 (MQ=255)
                                                                                                                               cTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGaa                        <  1:4955195/150‑1 (MQ=255)
                                                                                                                               cTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGaa                        <  2:192613/150‑1 (MQ=255)
                                                                                                                                tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCAc                                                          <  1:4317109/115‑1 (MQ=255)
                                                                                                                                tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg                       >  2:4971390/1‑150 (MQ=255)
                                                                                                                                tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg                       >  2:2792212/1‑150 (MQ=255)
                                                                                                                                 gCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAgg                      <  2:2461401/150‑1 (MQ=255)
                                                                                                                                 gCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAgg                      >  2:1699958/1‑150 (MQ=255)
                                                                                                                                     ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                  >  2:3334061/1‑150 (MQ=255)
                                                                                                                                     ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                  >  2:586901/1‑150 (MQ=255)
                                                                                                                                     ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                  >  2:4061377/1‑150 (MQ=255)
                                                                                                                                         cccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACGAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACTCTGCGCCAGGCTTCGAAGGCGCTCTgg              <  1:4060600/150‑1 (MQ=255)
                                                                                                                                          ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTggg             <  1:728906/150‑1 (MQ=255)
                                                                                                                                              gCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTg         >  2:3110413/1‑150 (MQ=255)
                                                                                                                                                cTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGAt       >  1:3805063/1‑150 (MQ=255)
                                                                                                                                                     tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGTACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg  >  2:3822892/1‑150 (MQ=255)
                                                                                                                                                     tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg  >  1:4356336/1‑150 (MQ=255)
                                                                                                                                                     tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg  >  2:3657603/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATG  >  NC_002947/6005479‑6005777

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: