Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,488,222 T→C 5.1% intergenic (+35/‑367) PP_2182 → / → PP_2183 hypothetical protein/formate dehydrogenase subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,488,2220TC 5.1% 231.9 / 2.0 79intergenic (+35/‑367)PP_2182/PP_2183hypothetical protein/formate dehydrogenase subunit gamma
Reads supporting (aligned to +/- strand):  ref base T (32/43);  new base C (2/2);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGG  >  NC_002947/2488077‑2488370
                                                                                                                                                  |                                                                                                                                                     
cagGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCa                                                                                                                                                   <  1:905834/150‑1 (MQ=255)
  ggCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGc                                                                                                                                                 >  1:1131764/1‑150 (MQ=255)
        gcgGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGACGGCa                                                                                                                                                   <  2:957047/142‑1 (MQ=255)
        gcgGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGACGGCa                                                                                                                                                   >  1:957047/1‑142 (MQ=255)
         cgGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGc                                                                                                                                          <  1:1331932/150‑1 (MQ=255)
          ggCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGcc                                                                                                                                         >  1:254185/1‑150 (MQ=255)
            cAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCgg                                                                                                                                       <  2:888367/150‑1 (MQ=255)
            cAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCgg                                                                                                                                       <  2:1131764/150‑1 (MQ=255)
              aCGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCgggg                                                                                                                                     <  1:493622/150‑1 (MQ=255)
              aCGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCgggg                                                                                                                                     <  2:254185/150‑1 (MQ=255)
                gCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCa                                                                                                                                              <  1:397885/139‑1 (MQ=255)
                gCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCa                                                                                                                                              >  2:397885/1‑139 (MQ=255)
                   gCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGcc                                                                                                                                >  1:670606/1‑150 (MQ=255)
                    ccAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCg                                                                                                                              >  1:206443/1‑150 (MQ=255)
                        gACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTg                                                                                                                          >  1:1261964/1‑150 (MQ=255)
                         aCCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTgg                                                                                                                         <  1:1538116/150‑1 (MQ=255)
                         aCCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTgg                                                                                                                         <  2:1603938/150‑1 (MQ=255)
                             aCGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAgg                                                                                                                     <  2:1202303/150‑1 (MQ=255)
                                   tCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAg                                                                                                                                                  <  1:1583003/116‑1 (MQ=255)
                                   tCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAg                                                                                                                                                  >  2:1583003/1‑116 (MQ=255)
                                      ggCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGCCTGGTTTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCACATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAAt                                                                                                            <  1:155072/149‑1 (MQ=255)
                                       gCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATg                                                                                                           <  2:1602978/150‑1 (MQ=255)
                                       gCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATg                                                                                                           <  2:1626221/150‑1 (MQ=255)
                                        cGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACAACGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATgg                                                                                                          <  1:579122/150‑1 (MQ=255)
                                         gttgTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGa                                                                                                         >  2:1449701/4‑150 (MQ=255)
                                          aaGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGAc                                                                                                        <  2:1003644/150‑1 (MQ=255)
                                          aaGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGAc                                                                                                        <  2:670606/150‑1 (MQ=255)
                                           aGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGTCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGGCt                                                                                                       <  2:706653/150‑1 (MQ=37)
                                           aGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCa                                                                                                                                              >  1:322406/1‑112 (MQ=255)
                                           aGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCa                                                                                                                                              <  2:322406/112‑1 (MQ=255)
                                           aGTCAAACGCCAAGACCCAGACTCCGGCCCCCGCGCCTTCCCCCCAAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGGCt                                                                                                       <  2:123277/150‑1 (MQ=21)
                                              cAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGc                                                                                                                        <  1:359775/130‑1 (MQ=255)
                                              cAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGc                                                                                                                        >  2:359775/1‑130 (MQ=255)
                                               aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGt                                                                                                   >  1:1424274/1‑150 (MQ=255)
                                               aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGt                                                                                                   >  2:1732953/1‑150 (MQ=255)
                                                         gcccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCga                                                                                         <  2:869129/150‑1 (MQ=255)
                                                         gcccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCga                                                                                         <  1:1732953/150‑1 (MQ=255)
                                                               gcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCa                                                                                                                       >  2:1580908/1‑114 (MQ=255)
                                                               gcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCa                                                                                                                       <  1:1580908/114‑1 (MQ=255)
                                                                cTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTg                                                                                  >  1:1657402/1‑150 (MQ=255)
                                                                 tCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTgg                                                                                 >  1:1703252/1‑150 (MQ=255)
                                                                 tCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGCTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTgg                                                                                 <  1:1650991/150‑1 (MQ=255)
                                                                            ccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACg                                                                      <  1:1347979/150‑1 (MQ=255)
                                                                              gccgACACCTCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGCATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGTTTTCGAGAGCCTGGTGCATGCAACGAt                                                                    <  2:533509/150‑1 (MQ=255)
                                                                              gccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGAt                                                                    <  1:1599061/150‑1 (MQ=255)
                                                                              gccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGAt                                                                    <  2:795283/150‑1 (MQ=255)
                                                                               ccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGATGCCCAGTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATg                                                                   <  1:724608/150‑1 (MQ=255)
                                                                                 gACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGt                                                                 >  2:751176/1‑150 (MQ=255)
                                                                                       gCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGACGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGttt                                                                                            >  2:905102/1‑117 (MQ=255)
                                                                                       gCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGACGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGttt                                                                                            <  1:905102/117‑1 (MQ=255)
                                                                                          aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGAc                                                         <  1:759356/150‑1 (MQ=255)
                                                                                              gTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCAc                                                     <  1:317940/150‑1 (MQ=255)
                                                                                                 cAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcat                                                                            <  2:1852931/124‑1 (MQ=255)
                                                                                                 cAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcat                                                                            >  1:1852931/1‑124 (MQ=255)
                                                                                                   gTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGAt                                                >  1:691999/1‑150 (MQ=255)
                                                                                                    tCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGCCCCGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATc                                               >  1:1087799/1‑150 (MQ=255)
                                                                                                      tACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCgg                                             >  1:62109/1‑150 (MQ=255)
                                                                                                      acgcGACTGAGCTCAGTCGCGTTCTGCCCCGTAGATATGCTGGCTGCATGCTGGTTGCCGGGTCGGCCTTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCgg                                             <  2:1390462/147‑1 (MQ=255)
                                                                                                        ccGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGGGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCa                                           >  2:1398943/1‑150 (MQ=255)
                                                                                                             tGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATgggct                                                                                                       >  1:869535/1‑82 (MQ=17)
                                                                                                             tGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGGCt                                                                                                       <  2:869535/85‑1 (MQ=17)
                                                                                                              gAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGt                                     <  2:1720607/150‑1 (MQ=255)
                                                                                                                gCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCCAGACCCACCTGATCGGCAAGCCGTgg                                   <  1:751176/150‑1 (MQ=255)
                                                                                                                 cTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCAGTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCGGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGc                                  >  1:1699545/1‑150 (MQ=255)
                                                                                                                   cAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAg                                <  2:1261964/150‑1 (MQ=255)
                                                                                                                      tCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCCTTGCAGCGACAGCCGTCGGGTAGGCCCGGCAATGGTCGGGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGCTCGGCAAGCCGTGGCAgatg                             >  1:1208215/1‑150 (MQ=255)
                                                                                                                      tCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAgatg                             >  2:829540/1‑150 (MQ=255)
                                                                                                                      tCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGGTTCAAGACCCACCTGATCGGCAAGCCGTGGCAgatg                             >  1:1273212/1‑150 (MQ=255)
                                                                                                                        gcgTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCAtt                                                                                                                                            >  2:1933772/1‑38 (MQ=255)
                                                                                                                        gcgTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCAtt                                                                                                                                            <  1:1933772/38‑1 (MQ=255)
                                                                                                                        gcgTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAgatgat                           >  2:1434837/1‑150 (MQ=255)
                                                                                                                              tGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCggg                     <  2:206443/150‑1 (MQ=255)
                                                                                                                                  ccGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCCGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGt                 >  1:1867108/1‑150 (MQ=255)
                                                                                                                                        ccATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCAc           >  1:854267/1‑150 (MQ=255)
                                                                                                                                          aTGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGGCCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt         <  1:1791648/150‑1 (MQ=255)
                                                                                                                                           tGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTg        >  2:1777343/1‑150 (MQ=255)
                                                                                                                                             cTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGc                                            <  1:421090/112‑1 (MQ=255)
                                                                                                                                             cTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGc                                            >  2:421090/1‑112 (MQ=255)
                                                                                                                                                 aTGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCgg  <  2:691999/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAGGCTGCGCGGCAACGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGG  >  NC_002947/2488077‑2488370

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: