Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,629 G→A 10.7% 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,6290GA10.7% 727.4 / 58.6 261Q96* (CAG→TAG) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (134/99);  new base A (14/14);  total (148/113)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 271 total aligned reads displayed.

CACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATG  >  NC_002947/6005480‑6005777
                                                                                                                                                     |                                                                                                                                                    
 aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     >  1:2575689/1‑150 (MQ=255)
 aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  2:4284931/150‑1 (MQ=255)
 aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  2:3307230/150‑1 (MQ=255)
 aCGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTgg                                                                                                                                                     <  1:177874/150‑1 (MQ=255)
  cGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTggg                                                                                                                                                    <  2:4716350/150‑1 (MQ=255)
   gCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGt                                                                                                                                                   >  2:4242978/1‑150 (MQ=255)
      gACCACACCCAGCGGGTGCCCGGTGTCCCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGAAGGCCTGGGTGaa                                                                                                                                                <  2:3748271/150‑2 (MQ=255)
         cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                                                                                                                                             <  1:973969/150‑1 (MQ=255)
         cacaCCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCtt                                                                                                                                             >  2:4182793/1‑150 (MQ=255)
          acacCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTg                                                                                                                                            >  1:2980561/1‑150 (MQ=255)
           cacCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGc                                                                                                                                           >  2:2619688/1‑150 (MQ=255)
            acCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGcc                                                                                                                                          >  1:3652105/1‑150 (MQ=255)
              ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCgg                                                                                                                                        >  2:2602892/1‑150 (MQ=255)
              ccAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCgg                                                                                                                                        >  2:1120968/1‑150 (MQ=255)
                aGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCt                                                                                                                                      <  1:4142522/150‑1 (MQ=255)
                aGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCt                                                                                                                                      >  2:1394771/1‑150 (MQ=255)
                   gggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcg                                                                                                                                   >  2:3464955/1‑150 (MQ=255)
                    ggTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcga                                                                                                                                  <  1:2636964/150‑1 (MQ=255)
                    gatGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcga                                                                                                                                  <  2:2918516/148‑1 (MQ=255)
                    gatGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcga                                                                                                                                  >  1:705349/3‑150 (MQ=255)
                    gatGCCAGGTTTCCATCATACGATGGCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTtcga                                                                                                                                  <  2:360364/148‑1 (MQ=255)
                      tGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatc                                                                                                                                <  2:619250/150‑1 (MQ=255)
                        ccAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTtcgatcga                                                                                                                              <  1:806358/150‑1 (MQ=255)
                         cAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCgag                                                                                                                             >  2:4787118/1‑150 (MQ=255)
                           ggTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGAc                                                                                                                           <  1:318377/150‑1 (MQ=255)
                             tttCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCa                                                                                                                         >  2:1729238/1‑150 (MQ=255)
                                 aCGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTc                                                                                                                     >  1:4093920/1‑150 (MQ=255)
                                       acgatggCCAGCGCGTTCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAg                                                                                                               <  2:151784/146‑1 (MQ=255)
                                            ggCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCgg                                                                                                          >  2:2147924/1‑150 (MQ=255)
                                               ccGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCATCTCGCCGAGGTCGGCCt                                                                                                       >  2:5062981/1‑150 (MQ=255)
                                                   gcgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgtg                                                                                                   <  1:1380509/150‑1 (MQ=255)
                                                    cgttCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTtggaa                                                                                                  <  2:4522445/147‑4 (MQ=255)
                                                    cgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGc                                                                                                  <  2:2724596/150‑1 (MQ=255)
                                                      ttCGGTAGCGATTGTGAGGCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTtggaagc                                                                                                >  1:571423/2‑145 (MQ=255)
                                                       tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa                                                                                               >  1:3875566/1‑150 (MQ=255)
                                                       tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa                                                                                               >  2:4983420/1‑150 (MQ=255)
                                                         ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg                                                                                             <  1:359043/150‑1 (MQ=255)
                                                         ggAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACg                                                                                             >  2:3925326/1‑150 (MQ=255)
                                                               gATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCAc                                                                                       >  1:4802546/1‑150 (MQ=255)
                                                                     cAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGcc                                                                                 <  1:1587609/150‑1 (MQ=255)
                                                                        gCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGaa                                                                              >  2:771225/1‑150 (MQ=255)
                                                                        gCCGTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGag                                                                              <  2:1668557/150‑2 (MQ=255)
                                                                         ccGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                             >  2:684559/1‑150 (MQ=255)
                                                                           gTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCATGTACTTCGCCGATGCCTTCCGAAGGAATCTTGCAAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAg                                                                           >  1:4596882/1‑150 (MQ=255)
                                                                             aCAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGgc                                                                         <  1:1452615/150‑1 (MQ=255)
                                                                                 cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                     <  1:1078374/150‑1 (MQ=255)
                                                                                 cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                     >  2:909694/1‑150 (MQ=255)
                                                                                 cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                     >  1:3269267/1‑150 (MQ=255)
                                                                                   ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                   <  1:4058171/150‑1 (MQ=255)
                                                                                      gcgACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGtt                                                                >  1:1199593/1‑150 (MQ=255)
                                                                                           agaccAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCa                                                           >  2:3440843/4‑150 (MQ=255)
                                                                                           aGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCa                                                           <  1:3189900/150‑1 (MQ=255)
                                                                                             gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACg                                                         >  2:2414534/1‑150 (MQ=255)
                                                                                              ccAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACgg                                                        <  1:198348/150‑1 (MQ=255)
                                                                                              ccAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACgg                                                        >  2:2224208/1‑150 (MQ=255)
                                                                                               cAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgct                                                       >  1:2424698/1‑148 (MQ=255)
                                                                                                aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcg                                                      <  1:503025/150‑1 (MQ=255)
                                                                                                aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGGTCGCCACGgcg                                                      >  2:4417800/1‑150 (MQ=255)
                                                                                                aaCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCGCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcg                                                      >  2:2966808/1‑150 (MQ=255)
                                                                                                 actgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgcttt                                                     >  2:571976/4‑146 (MQ=255)
                                                                                                 aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcgc                                                     <  2:1866340/150‑1 (MQ=255)
                                                                                                 aCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGgcgc                                                     <  1:2225373/150‑1 (MQ=255)
                                                                                                   tgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttgg                                                   >  1:3933014/2‑144 (MQ=255)
                                                                                                   tgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttgg                                                   >  1:2809887/2‑144 (MQ=255)
                                                                                                    gCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGa                                                  <  2:3373986/150‑1 (MQ=255)
                                                                                                    gCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttggt                                                  <  2:557770/150‑8 (MQ=255)
                                                                                                      gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                <  2:3146112/150‑1 (MQ=255)
                                                                                                      gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                >  1:817235/1‑150 (MQ=255)
                                                                                                      gAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGc                                                >  2:1927937/1‑150 (MQ=255)
                                                                                                         gTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCgg                                             <  1:4529400/150‑1 (MQ=255)
                                                                                                             cAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCGCGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGcaca                                         >  2:3960384/1‑150 (MQ=255)
                                                                                                                     gATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                                 <  1:4478533/150‑1 (MQ=255)
                                                                                                                     gATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                                 >  2:2052969/1‑150 (MQ=255)
                                                                                                                     gATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                                 >  1:3574958/1‑150 (MQ=255)
                                                                                                                     gATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                                 >  1:1651546/1‑150 (MQ=255)
                                                                                                                      aTCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAg                                >  2:1779962/1‑150 (MQ=255)
                                                                                                                       tCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAgg                               >  1:5078813/1‑150 (MQ=255)
                                                                                                                       tCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAgg                               >  1:3492196/1‑150 (MQ=255)
                                                                                                                          tttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGcc                                                                                                        >  1:228806/1‑74 (MQ=255)
                                                                                                                           ttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTc                           <  1:4174258/150‑1 (MQ=255)
                                                                                                                           ttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTc                           >  2:1596875/1‑150 (MQ=255)
                                                                                                                             ccTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGa                         <  1:4221521/150‑1 (MQ=255)
                                                                                                                                   gatCACCCAGGAAGGCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcgc                   >  1:863474/3‑150 (MQ=255)
                                                                                                                                   cTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcgc                   >  2:509284/1‑150 (MQ=255)
                                                                                                                                   cTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcgc                   >  2:975630/1‑150 (MQ=255)
                                                                                                                                   cTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcgc                   >  2:3246535/1‑150 (MQ=255)
                                                                                                                                    ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                  >  1:1634613/1‑150 (MQ=255)
                                                                                                                                    ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                  >  2:1285256/1‑150 (MQ=255)
                                                                                                                                     tCACCCAGGAAGGCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGctc                 >  2:3328541/1‑150 (MQ=255)
                                                                                                                                       aCCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTg               >  2:3721828/1‑150 (MQ=255)
                                                                                                                                       aCCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTg               >  2:2005166/1‑150 (MQ=255)
                                                                                                                                        cccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTgg              <  1:5086088/150‑1 (MQ=255)
                                                                                                                                           aGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGcc           >  1:205193/1‑150 (MQ=255)
                                                                                                                                           aGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGcc           <  1:4865569/150‑1 (MQ=255)
                                                                                                                                              ccTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGa        <  1:1886046/150‑1 (MQ=255)
                                                                                                                                                  ccTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGa    <  2:2402465/150‑1 (MQ=255)
                                                                                                                                                   cTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGAt   <  1:2602892/150‑1 (MQ=255)
                                                                                                                                                   cTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGAt   >  2:2866299/1‑150 (MQ=255)
                                                                                                                                                    tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg  >  2:3501053/1‑150 (MQ=255)
                                                                                                                                                    tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg  <  1:1063328/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACGCCGACCACACCCAGCGGGTGCCAGGTTTCACGCATGTGGTGGCCCGGGCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATG  >  NC_002947/6005480‑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: