Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,829,581 C→T 6.3% A17A (GCG→GCA)  PP_3382 ← gluconate 2‑dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,829,5810CT 6.3% 236.5 / 4.7 79A17A (GCG→GCA) PP_3382gluconate 2‑dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base T (2/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

ACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGCCCGGTCGG  >  NC_002947/3829433‑3829730
                                                                                                                                                    |                                                                                                                                                     
ccacCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCg                                                                                                                                                      <  2:1151887/149‑1 (MQ=255)
 cacCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCAGTCGCTACTGCCTGTGCcgccgc                                                                                                                                                     >  1:236372/1‑150 (MQ=255)
  acCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCa                                                                                                                                                    >  2:149812/1‑150 (MQ=255)
       gCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAGgc                                                                                                                                               <  1:461162/150‑1 (MQ=255)
         cAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGACGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGgcgc                                                                                                                                             >  2:802857/1‑150 (MQ=255)
            tCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCa                                                                                                                                          >  2:380822/1‑150 (MQ=255)
            tCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCa                                                                                                                                          >  1:345354/1‑150 (MQ=255)
                 gTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCAcc                                                                                                                                     >  2:871418/1‑150 (MQ=255)
                    tCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAgg                                                                                                                                  <  2:720490/150‑1 (MQ=255)
                    tCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAgg                                                                                                                                  <  1:622491/150‑1 (MQ=255)
                     ccATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAggg                                                                                                                                 >  1:1013172/1‑150 (MQ=255)
                       aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGc                                                                                                                                                     >  2:857608/1‑128 (MQ=255)
                       aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGc                                                                                                                                                     <  1:857608/128‑1 (MQ=255)
                       aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTAAGTGCTCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTg                                                                                                                               >  1:430652/1‑150 (MQ=255)
                        tCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCNTNATTCGCNACGGCCTGTGCANCCGCACCGGCGCCCAGCACCAGGGNgg                                                                                                                              <  1:1042174/150‑1 (MQ=255)
                         cGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGc                                                                                                                             <  2:102361/150‑1 (MQ=255)
                            cAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAt                                                                                                                          >  2:387974/1‑150 (MQ=255)
                               cccGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCaa                                                                                                                       <  1:1065283/150‑1 (MQ=255)
                                      gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACa                                                                                                                     >  1:578477/1‑145 (MQ=255)
                                      gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                                <  1:149812/150‑1 (MQ=255)
                                      gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                                >  2:604303/1‑150 (MQ=255)
                                      gCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGACGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                                <  1:802857/150‑1 (MQ=255)
                                      gCAAAGGGTTTGCCGCCTTTGGCGGGATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACa                                                                                                                     <  2:578477/145‑1 (MQ=255)
                                        aaaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAgg                                                                                                                                                <  2:878403/116‑1 (MQ=255)
                                        aaaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAgg                                                                                                                                                >  1:878403/1‑116 (MQ=255)
                                         aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGTGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                             >  2:1262600/1‑150 (MQ=255)
                                         aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                             >  2:31610/1‑150 (MQ=255)
                                         aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                             <  1:899430/150‑1 (MQ=255)
                                         aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                             >  2:1412365/1‑150 (MQ=255)
                                           gggTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCAt                                                                                                           >  2:969666/1‑150 (MQ=255)
                                                tGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCa                                                                                                      <  1:1407737/150‑1 (MQ=255)
                                                 gccgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCAt                                                                                                     <  1:871418/150‑1 (MQ=255)
                                                   cgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGc                                                                                                                             <  1:412650/124‑1 (MQ=255)
                                                   cgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGc                                                                                                                             >  2:412650/1‑124 (MQ=255)
                                                    gccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGc                                                                                                  <  1:31610/150‑1 (MQ=255)
                                                      cTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGTGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGctc                                                                                                <  1:1262600/150‑1 (MQ=255)
                                                           gCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCaa                                                                                           <  2:285791/150‑1 (MQ=255)
                                                                  ggCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGa                                                                                    <  2:1013172/150‑1 (MQ=255)
                                                                      ggCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGgcgc                                                                                                                                             <  1:1196222/89‑1 (MQ=255)
                                                                      ggCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGgcgc                                                                                                                                             >  2:1196222/1‑89 (MQ=255)
                                                                          aCGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCAcag                                                                                                                                                 <  1:1404559/81‑3 (MQ=255)
                                                                          aCGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCAcag                                                                                                                                                 >  2:1404559/1‑79 (MQ=255)
                                                                             cAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCg                                                                         <  2:345354/150‑1 (MQ=255)
                                                                               gccaacGGCGCGGGCCAGCTACTCGCCTAGGAGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACTGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCggg                                                                       <  2:263130/145‑1 (MQ=255)
                                                                                  aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCaa                                                                                >  2:1525472/1‑138 (MQ=255)
                                                                                  aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCaa                                                                                <  1:1525472/138‑1 (MQ=255)
                                                                                     ggCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGttttt                                                                 <  2:1492271/150‑1 (MQ=255)
                                                                                        gcgGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAg                                                              <  2:228860/150‑1 (MQ=255)
                                                                                              cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTc                                                                <  2:91141/142‑1 (MQ=255)
                                                                                              cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTc                                                                >  1:91141/1‑142 (MQ=255)
                                                                                              cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTg                                                        >  1:67935/1‑150 (MQ=255)
                                                                                                  cacgcGCCCATGCGTACCTTCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCa                                                    <  2:354451/146‑1 (MQ=255)
                                                                                                   aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCAc                                                   <  1:387974/150‑1 (MQ=255)
                                                                                                   aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCAc                                                   >  2:326822/1‑150 (MQ=255)
                                                                                                       gCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAACAACGTTGGGTTCATGCGCATCGGTCCTTAAGGCTGGACGAACGGACCGGGGTTTTTCAGGTACTGCTCACGGAt                                               >  1:1102679/1‑150 (MQ=255)
                                                                                                         ccAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATgg                                             >  1:1232491/1‑150 (MQ=255)
                                                                                                               gTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCggg                                       >  2:1557965/1‑150 (MQ=255)
                                                                                                               gTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCggg                                       <  2:162078/150‑1 (MQ=255)
                                                                                                                 aCCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGcc                                     <  2:966064/150‑1 (MQ=255)
                                                                                                                  ccTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCg                                    >  1:1580083/1‑150 (MQ=255)
                                                                                                                    tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                                  <  2:67935/150‑1 (MQ=255)
                                                                                                                    tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                                  <  2:1232491/150‑1 (MQ=255)
                                                                                                                      cgcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                                <  2:1580083/150‑1 (MQ=255)
                                                                                                                      cgcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                                <  1:1412365/150‑1 (MQ=255)
                                                                                                                        cATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATc                                                                                                                         >  1:1389054/1‑59 (MQ=255)
                                                                                                                        cATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATc                                                                                                                         <  2:1389054/59‑1 (MQ=255)
                                                                                                                              cATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCCTGCTCCTCGCTCCTCACGCCTGGACCATCAATCCGGTGTAATTCAGGAACATCACTCTAATATGGCGGTCCGACCAGTACGTGa                        <  2:430652/150‑1 (MQ=255)
                                                                                                                               aTTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                                   >  1:1151597/1‑90 (MQ=255)
                                                                                                                               aTTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                                   <  2:1151597/90‑1 (MQ=255)
                                                                                                                                   gCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATTCGCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCa                   <  1:604303/150‑1 (MQ=255)
                                                                                                                                    cTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGt                                                                     >  2:21331/1‑99 (MQ=255)
                                                                                                                                    cTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGt                                                                     <  1:21331/99‑1 (MQ=255)
                                                                                                                                     tACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                                >  2:1318446/1‑135 (MQ=255)
                                                                                                                                     tACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGCCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGCTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGc                 >  1:1569630/1‑150 (MQ=255)
                                                                                                                                     tACAGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCa                                <  1:1318446/135‑1 (MQ=255)
                                                                                                                                       cGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGCTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGcca               >  2:165275/1‑150 (MQ=255)
                                                                                                                                        ggCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTg                                                                                      <  1:1473101/78‑1 (MQ=255)
                                                                                                                                        ggCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTg                                                                                      >  2:1473101/1‑78 (MQ=255)
                                                                                                                                            tgtgCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGc          >  1:1302217/1‑150 (MQ=255)
                                                                                                                                            tgtgCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGc          >  2:1198434/1‑150 (MQ=255)
                                                                                                                                                   ccGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGCCCGGTCg   >  2:1292966/1‑150 (MQ=255)
                                                                                                                                                    cGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGCCCGGTCgg  <  1:1557965/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCACCAGCCCGGTCGG  >  NC_002947/3829433‑3829730

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