Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,844,168 A→G 6.2% D245G (GAC→GGC)  PP_5122 → GMC family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,844,1680AG 6.2% 279.4 / 6.2 97D245G (GAC→GGC) PP_5122GMC family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base G (3/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

CCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCTTCCTG  >  NC_002947/5844025‑5844316
                                                                                                                                               |                                                                                                                                                    
cccGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGNCCa                                                                                                                                                <  2:1196791/150‑1 (MQ=255)
 ccGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAg                                                                                                                                               <  1:252396/150‑1 (MQ=255)
  cGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAgg                                                                                                                                              <  1:1357956/150‑1 (MQ=255)
  cGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGAATTCAGAGCCTGACATGCCAGGCCCTGGGCGTCCAgg                                                                                                                                              <  2:1667890/150‑1 (MQ=255)
   ggTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAggg                                                                                                                                             <  2:172113/150‑1 (MQ=255)
    gTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGc                                                                                                                                            >  2:1886269/1‑150 (MQ=255)
    gTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGc                                                                                                                                            >  2:139514/1‑150 (MQ=255)
      gAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGGCGTCCAGGGcct                                                                                                                                          <  1:1228618/150‑3 (MQ=255)
         tGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATaca                                                                                                                                       <  1:535189/150‑1 (MQ=255)
         tGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATaca                                                                                                                                       >  1:1863868/1‑150 (MQ=255)
           ccAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACAcc                                                                                                                                     >  2:1529420/1‑150 (MQ=255)
            cAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCACGCCCTGGACGCCCAGGGCATACAccc                                                                                                                                    >  1:486951/1‑150 (MQ=255)
              aaCAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTa                                                                                                                                  <  1:1228338/150‑1 (MQ=255)
              aaCAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTa                                                                                                                                  >  2:761772/1‑150 (MQ=255)
               aCAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTAc                                                                                                                                 <  2:700301/150‑1 (MQ=255)
                  aTCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCGGGACGCCCAGGGCATACACCCTACCgg                                                                                                                              >  2:400718/1‑150 (MQ=255)
                      aTGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggcagg                                                                                                                          <  1:421120/150‑1 (MQ=255)
                                cccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCg                                                                                                                <  2:394731/150‑1 (MQ=255)
                                 ccGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGt                                                                                                               >  2:478747/1‑150 (MQ=255)
                                  cGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGTATTCAGAGCCTGACATGCCAGCCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTg                                                                                                              <  2:979417/150‑1 (MQ=255)
                                  cGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTg                                                                                                              >  2:323997/1‑150 (MQ=255)
                                    tataCCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCa                                                                                                            >  1:916078/1‑150 (MQ=255)
                                             cACCCTCGGGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCCTGCACGctg                                                                                                   <  2:1041359/150‑3 (MQ=255)
                                                ccTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTa                                                                                                >  2:1216184/1‑150 (MQ=255)
                                                  tCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAg                                                                                                                                               <  1:1851883/101‑1 (MQ=255)
                                                  tCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAg                                                                                                                                               >  2:1851883/1‑101 (MQ=255)
                                                   cGAGCAAGGCGGCAAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCAACACGGCGAACGCATTCAGAGCCTGTCATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCATGTCAGCGTGCAGGCACGCCATTAcct                                                                                             <  2:578574/150‑1 (MQ=255)
                                                    gAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                                            >  2:813962/1‑150 (MQ=255)
                                                    gAGCAACGCCGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTAcctc                                                                                            >  2:417954/1‑150 (MQ=255)
                                                        aaGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTgg                                                                                        >  2:35613/1‑150 (MQ=255)
                                                          ggcggcGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACAcc                                                                                                                                     <  1:226735/103‑1 (MQ=255)
                                                          ggcggcGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACAcc                                                                                                                                     >  2:226735/1‑103 (MQ=255)
                                                           gcggcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCg                                                                                     >  1:517340/1‑150 (MQ=255)
                                                           gcggcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCg                                                                                     >  1:1044077/1‑150 (MQ=255)
                                                            cggcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCgg                                                                                    <  2:517340/150‑1 (MQ=255)
                                                              gcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGgcaggcag                                                                                                                       <  2:772270/113‑1 (MQ=255)
                                                              gcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGgcaggcag                                                                                                                       <  2:772269/113‑1 (MQ=255)
                                                              gcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGgcaggcag                                                                                                                       >  1:772270/1‑113 (MQ=255)
                                                              gcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGgcaggcag                                                                                                                       >  1:772269/1‑113 (MQ=255)
                                                              gcgAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTg                                                                                  <  2:539303/150‑1 (MQ=255)
                                                               cgAACTGCGGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTgg                                                                                 <  1:861464/150‑1 (MQ=255)
                                                                gAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGc                                                                                <  1:139514/150‑1 (MQ=255)
                                                                gAACTGCTGTATCTGGCGCGGGCCGAGAGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCAAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTggg                                                                                >  1:647933/1‑149 (MQ=255)
                                                                 aaCTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                                               <  1:794548/150‑1 (MQ=255)
                                                                 aaCTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                                               <  1:11437/150‑1 (MQ=255)
                                                                 aaCTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                                               <  2:199712/150‑1 (MQ=255)
                                                                 aaCTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                                               <  1:1886269/150‑1 (MQ=255)
                                                                          aTCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACa                                                                      <  1:1529420/150‑1 (MQ=255)
                                                                           tCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggca                                                                                                                            <  2:407528/95‑1 (MQ=255)
                                                                           tCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCggca                                                                                                                            >  1:407528/1‑95 (MQ=255)
                                                                                cgcgGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCgg                                                                <  1:761772/150‑1 (MQ=255)
                                                                                 gcgGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGc                                                               >  2:1761763/1‑150 (MQ=255)
                                                                                 gcgGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGc                                                               >  2:860106/1‑150 (MQ=255)
                                                                                 gcgGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGc                                                               >  2:457483/1‑150 (MQ=255)
                                                                                    ggCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCct                                                            <  1:655668/150‑1 (MQ=255)
                                                                                       cGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAg                                                                                                                                               >  1:530137/1‑64 (MQ=255)
                                                                                       cGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAg                                                                                                                                               <  2:530137/64‑1 (MQ=255)
                                                                                         aGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctgctgc                                                       <  1:860106/150‑1 (MQ=255)
                                                                                              ttCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTgcgc                                                  >  1:998822/1‑150 (MQ=255)
                                                                                                 gCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGctct                                               >  1:40325/1‑150 (MQ=255)
                                                                                                   tCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGc                                             <  2:1044077/150‑1 (MQ=255)
                                                                                                     acaACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCcg                                           >  1:306547/1‑150 (MQ=255)
                                                                                                       aaCGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCcgcg                                         >  1:1331422/1‑150 (MQ=255)
                                                                                                             gAACGCATTCAGAGCCTGACATGCCAGGCCCTGAACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAc                                   >  2:138991/1‑150 (MQ=255)
                                                                                                              aaCGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTGTCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAcc                                  <  2:1266779/150‑1 (MQ=255)
                                                                                                               aCGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGAccc                                 <  1:1226876/150‑1 (MQ=255)
                                                                                                                  cATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctg                                                           >  2:651981/1‑121 (MQ=255)
                                                                                                                  cATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctg                                                           <  1:651981/121‑1 (MQ=255)
                                                                                                                   aTTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCAt                             >  2:591780/1‑150 (MQ=255)
                                                                                                                    ttCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCCCTGCCGCGCCCGACCCGCATg                            >  1:1608106/1‑150 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCTGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGCCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTc                         <  2:998822/150‑1 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCAc                                                                                                       <  2:415156/72‑1 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCAc                                                                                                       >  1:415156/1‑72 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctgctgct                                                      >  1:1280617/1‑121 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCctgctgct                                                      <  2:1280617/121‑1 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTc                         >  1:1747131/1‑150 (MQ=255)
                                                                                                                       agagCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTc                         <  2:595872/150‑1 (MQ=255)
                                                                                                                         ggCCTGGCATGTCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGc                       >  1:1146327/2‑150 (MQ=255)
                                                                                                                           ccTGACGTGCCAGGCCCTGGACTCCCAGGGCATACACTCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCCTTGCCTGCTGGCCGGTGGCGGTATCCACAGCCCGGCCTTGCTGCGGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCt                     <  2:40325/150‑1 (MQ=255)
                                                                                                                            cTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTg                    >  1:23325/1‑150 (MQ=255)
                                                                                                                             tGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTgg                   >  1:1765980/1‑150 (MQ=255)
                                                                                                                             tGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGACTgg                   >  1:1012036/1‑150 (MQ=255)
                                                                                                                                   tgtatgCCCTGGGCGTCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTc                         <  1:1455802/133‑1 (MQ=255)
                                                                                                                                   tgtatgCCCTGGGCGTCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTc                         >  2:1455802/6‑138 (MQ=255)
                                                                                                                                    ccAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCTGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCg            >  2:1036457/1‑150 (MQ=255)
                                                                                                                                     tatgCCCTGGGCGTCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTTCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGa           >  1:1154488/4‑150 (MQ=255)
                                                                                                                                     tatgCCCTGGGCGTCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTAGCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGGCCCGCATGGGCGGCTGGGCAAGCGa           >  1:392125/4‑150 (MQ=255)
                                                                                                                                     cAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGa           >  2:1381665/1‑150 (MQ=255)
                                                                                                                                        gCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAAcc        <  2:486951/150‑1 (MQ=255)
                                                                                                                                           cTGGACGCCCAGGGCATACACCCTACTACCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGTCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCTTc     <  1:478747/150‑1 (MQ=255)
                                                                                                                                           cTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCTTc     <  2:1863868/150‑1 (MQ=255)
                                                                                                                                            tGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                                               >  1:868089/1‑75 (MQ=255)
                                                                                                                                            tGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCg                                                                               <  2:868089/75‑1 (MQ=255)
                                                                                                                                            tGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGc             >  1:45760/1‑141 (MQ=255)
                                                                                                                                            tGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGc             <  2:45760/141‑1 (MQ=255)
                                                                                                                                             ggACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCTTCCt   <  1:1448503/150‑1 (MQ=255)
                                                                                                                                              gACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCAGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCTTCCTg  >  2:243194/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCCGGTGAATGCCAAACAATCGATGCTGGTGACCCGTATACCTGCCACCCTCGAGCAAGGCGGCGAACTGCTGTATCTGGCGCGGGCCGAGCGCTTCGCTCACAACGGCGAACGCATTCAGAGCCTGACATGCCAGGCCCTGGACGCCCAGGGCATACACCCTACCGGCAGGCAGGTCAGCGTGCAGGCACGCCATTACCTCCTGGCCGGTGGCGGTATCAACAGCCCGGCCCTGCTGCTGCGCTCTGCCGCGCCCGACCCGCATGGTCGCCTGGGCAAGCGAACCTTCCTG  >  NC_002947/5844025‑5844316

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: