Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,926,153 Δ1 bp 13.2% coding (382/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,1530A.13.2% 222.4 / 29.0 76coding (382/1629 nt)peaAquinohemoprotein amine dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (28/38);  new base . (5/5);  total (33/43)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTT  >  NC_002947/3926008‑3926297
                                                                                                                                                 |                                                                                                                                                 
ccAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCTCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAgcgc                                                                                                                                               >  1:1109063/1‑150 (MQ=255)
 cAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTATCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCa                                                                                                                                             >  1:564495/1‑150 (MQ=255)
         aaGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                     >  2:5653/1‑150 (MQ=255)
         aaGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCg                                                                                                                                     <  1:381256/150‑1 (MQ=255)
          aGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtgga                                                                                                                                    <  2:657825/150‑4 (MQ=255)
          aGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtgga                                                                                                                                    <  2:1874499/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                                                                                                   <  2:539708/150‑1 (MQ=255)
                                           gcCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCa                                                                                                   >  1:2192360/1‑150 (MQ=255)
                                            ccGTTGCAAGGCGAGACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAg                                                                                                  <  2:2035358/150‑1 (MQ=255)
                                              gTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAt                                                                                                                                       >  1:2262144/1‑111 (MQ=255)
                                              gTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAt                                                                                                                                       <  2:2262144/111‑1 (MQ=255)
                                                 gCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAtt                                                                                             >  1:1581044/1‑150 (MQ=255)
                                                  cAAGGCGCCACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                            <  1:5653/150‑1 (MQ=255)
                                                  cAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                            <  1:1534563/150‑1 (MQ=255)
                                                  cAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                            <  2:1572592/150‑1 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTg                                                                                           >  2:620368/1‑150 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTg                                                                                           >  1:2019870/1‑150 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtg                                                                                                                                      >  1:831089/1‑106 (MQ=255)
                                                   aaGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTGCGCTACGCCAtg                                                                                                                                      <  2:831089/107‑2 (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                                                                                            >  1:1643187/1‑77 (MQ=255)
                                                                                                                           ttGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTAttt                                                                                            <  2:1643187/77‑1 (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)
                                                                                                                                                 |                                                                                                                                                 
CCAAGGTGGAAGTGGACCAGGTGTTCCCACTCCTTGGCCGGGCGCCGTTGCAAGGCGACACGGGCACCCGAGTGGCAACGGCCACAGGTCTCGCTGAGCTGGGTGTCGAAGTGCTCGACAGTATTGAGCCGACGCTCCATGGCGTAGCGCACGCCATCGGTTTCGCTGGGCGCCAGGCCCTGTTTGTCGGCCAGGTATTTGACCAGCGTGCGACGATCGTCGTCGCTGATCTGAAGGCCATGCATCACCTGCATGCGGGCGATGCTCATCAGCCAGCCTTCCGGTGTCTT  >  NC_002947/3926008‑3926297

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: