Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,926,157:1 +T 13.5% coding (378/1629 nt) peaA ← quinohemoprotein amine dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,926,1571.T13.5% 234.2 / 38.3 74coding (378/1629 nt)peaAquinohemoprotein amine dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base . (26/38);  new base T (5/5);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGT  >  NC_002947/3926008‑3926305
                                                                                                                                                      |                                                                                                                                                    
ccAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCTCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAgcgc                                                                                                                                                       >  1:1109063/1‑150 (MQ=255)
 cAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTATCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCa                                                                                                                                                     >  1:564495/1‑150 (MQ=255)
         aaGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                             <  1:381256/150‑1 (MQ=255)
         aaGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                             >  2:5653/1‑150 (MQ=255)
          aGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtgga                                                                                                                                            <  2:1874499/150‑4 (MQ=255)
          aGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtgga                                                                                                                                            <  2:657825/150‑4 (MQ=255)
          aGTGGAACAGCGGTTCCCACTCCTTGGCCGGGCGCGGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCgg                                                                                                                                            >  2:691876/1‑150 (MQ=255)
            tGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtggagc                                                                                                                                          <  2:1980266/150‑6 (MQ=255)
               aCCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCg                                                                                                                                       >  2:2345452/1‑150 (MQ=255)
                 cAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCt                                                                                                                                     >  1:238653/1‑150 (MQ=255)
                  aGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTg                                                                                                                                    <  1:2345452/150‑1 (MQ=255)
                   ggTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTgg                                                                                                                                   <  1:2146796/150‑1 (MQ=255)
                    gtgtTCCCGCTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTggg                                                                                                                                  >  1:1024676/1‑150 (MQ=255)
                    gtgtCCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTggg                                                                                                                                  <  1:903442/150‑1 (MQ=255)
                     tgtTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGCCGCGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgc                                                                                                                                 <  2:564495/150‑1 (MQ=255)
                     tgtTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgc                                                                                                                                 <  2:2387167/150‑1 (MQ=255)
                       ttCCCACTCCTTGGCCGGCCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGgcgc                                                                                                                               >  2:1818283/1‑150 (MQ=255)
                          ccACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAg                                                                                                                            >  2:1487793/1‑150 (MQ=255)
                          ccACTCCTTGGCCGGCCGCCGTTGCAAGACGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAg                                                                                                                            <  1:1818283/150‑1 (MQ=255)
                           cACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAgg                                                                                                                           <  1:1529186/150‑1 (MQ=255)
                                cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGCGACGGCACGCGATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                      >  2:1880411/1‑150 (MQ=255)
                                cTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCGCGCCATCGGTTTCGCTGGGCGCCAGGCCCTg                                                                                                                      >  1:1053620/1‑150 (MQ=255)
                                     ccGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGt                                                                                                                 <  2:2374101/150‑1 (MQ=255)
                                      cGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTc                                                                                                                <  1:320621/150‑1 (MQ=255)
                                       gggCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCg                                                                                                               <  1:2372628/150‑1 (MQ=255)
                                           gcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCa                                                                                                           >  1:2192360/1‑150 (MQ=255)
                                           gcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCa                                                                                                           <  2:539708/150‑1 (MQ=255)
                                            ccGTTGCAAGGCGAGACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAg                                                                                                          <  2:2035358/150‑1 (MQ=255)
                                              gTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAt                                                                                                                                               <  2:2262144/111‑1 (MQ=255)
                                              gTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAt                                                                                                                                               >  1:2262144/1‑111 (MQ=255)
                                                 gCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAtt                                                                                                     >  1:1581044/1‑150 (MQ=255)
                                                  cAAGGCGCCACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                                    <  1:5653/150‑1 (MQ=255)
                                                  cAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                                    <  2:1572592/150‑1 (MQ=255)
                                                  cAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                                    <  1:1534563/150‑1 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTg                                                                                                   >  1:2019870/1‑150 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTg                                                                                                   >  2:620368/1‑150 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtg                                                                                                                                              <  2:831089/107‑2 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtg                                                                                                                                              >  1:831089/1‑106 (MQ=255)
                                                    aGGCGACACGGGCACCCGGGTGGCCACGGCTGCTGTTCCCGCTTAGCTGGGTGTCGAAGTAATCACCAATCCAGTTCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGa                                                                                                  <  1:1224564/150‑1 (MQ=255)
                                                    aGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGa                                                                                                  <  2:1143139/150‑1 (MQ=255)
                                                           acGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCGCCTTCTGTCGGCCAGGTATTTGACCAGCGt                                                                                           <  1:1880411/150‑1 (MQ=255)
                                                                       gTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtcgtc                                                                               <  1:620368/150‑1 (MQ=255)
                                                                       gTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtcgtc                                                                               <  1:1478447/150‑1 (MQ=255)
                                                                            aaCGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGc                                                                          >  2:2370305/1‑150 (MQ=255)
                                                                               ggCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGa                                                                       <  2:1581044/150‑1 (MQ=255)
                                                                                     aGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACACCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGaa                                                                 >  1:1767391/1‑150 (MQ=255)
                                                                                       gTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAgg                                                               >  1:1867740/1‑150 (MQ=255)
                                                                                           cGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATTGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCACGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGACCTGAAGGCCAt                                                           <  2:1230015/150‑1 (MQ=255)
                                                                                            gCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGACAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATg                                                          >  2:1771789/1‑150 (MQ=255)
                                                                                                  cTGGGTGTCGACGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCAc                                                    <  1:1228804/150‑1 (MQ=255)
                                                                                                  cTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCAc                                                    <  2:1867740/150‑1 (MQ=255)
                                                                                                       tgtCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCa                                               >  1:165710/1‑150 (MQ=255)
                                                                                                       tgtCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCAGTTTGTCGGCCAGTTCCTTTTCCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCa                                               <  1:1487793/150‑1 (MQ=255)
                                                                                                           gAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCg                                           >  2:1935427/1‑150 (MQ=255)
                                                                                                                  tCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTTTCGGCCAGGTATTTGACCAGCGTGCGACGAtcgtcgtc                                                                            >  1:1749653/1‑110 (MQ=255)
                                                                                                                  tCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtcgtcgtc                                                                            <  2:1749653/110‑1 (MQ=255)
                                                                                                                  tCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATg                                    <  1:2262853/150‑1 (MQ=255)
                                                                                                                  tCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTAGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATg                                    <  1:691876/150‑1 (MQ=255)
                                                                                                                           ttGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                                    <  2:1643187/77‑1 (MQ=255)
                                                                                                                           ttGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                                    >  1:1643187/1‑77 (MQ=255)
                                                                                                                           ttGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATcagc                           >  1:1524843/1‑150 (MQ=255)
                                                                                                                           ttGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCAGGCATCACCTGCATGCGGGCGATGCTCATcagc                           <  2:1036568/150‑1 (MQ=255)
                                                                                                                               gCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtcgtcgtc                                                                            >  2:1937895/1‑97 (MQ=255)
                                                                                                                               gCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGAtcgtcgtc                                                                            <  1:1937895/97‑1 (MQ=255)
                                                                                                                               gCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCgg                                          <  2:886655/131‑1 (MQ=255)
                                                                                                                               gCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCgg                                          >  1:886655/1‑131 (MQ=255)
                                                                                                                                ccGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATg                                    >  2:2443998/1‑136 (MQ=255)
                                                                                                                                ccGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATg                                    <  1:2443998/136‑1 (MQ=255)
                                                                                                                                     gCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCg                 <  2:179486/150‑1 (MQ=255)
                                                                                                                                     gCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCg                 <  2:644581/150‑1 (MQ=255)
                                                                                                                                       tCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgt               >  2:489226/1‑150 (MQ=255)
                                                                                                                                        cgaTGGCGTGCGCTACGCCATGGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtg              >  2:1927403/3‑150 (MQ=255)
                                                                                                                                         gaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtgt             >  2:135910/2‑150 (MQ=255)
                                                                                                                                         gaTGGCGTGCGCTACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGgtgt             >  1:459507/2‑150 (MQ=255)
                                                                                                                                          aTGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGGGTGCGTCTTTCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTc            <  2:2019870/150‑1 (MQ=255)
                                                                                                                                            ggCGTAGCGCACGCCATAGGGTCCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGCCCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCtt          <  1:1305704/150‑1 (MQ=255)
                                                                                                                                                    gcACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCtggt  <  1:1096077/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTTGCGCTGGT  >  NC_002947/3926008‑3926305

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