Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,307 A→C 6.0% Q152P (CAG→CCG)  thiE → thiamine‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,451,3070AC 6.0% 914.9 / 8.2 292Q152P (CAG→CCG) thiEthiamine‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (150/117);  new base C (22/2);  total (173/119)
Fisher's exact test for biased strand distribution p-value = 3.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

CGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGCCGCGCCCGCGCCTTC  >  NC_002947/5451158‑5451455
                                                                                                                                                     |                                                                                                                                                    
 ggCTCCCCCTGCCACGCCAGCCTCGACCGGGCCGCCCTGGGTGCCAGCGACGGCGCCAGCTTCGTGGCGTTTGGCCGCTTGTTCACTTCGGTCCCCCCGCCGGGGGCGCCCGGTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAg                                                                                                                                                     <  2:6576700/150‑1 (MQ=255)
     ccACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGa                                                                                                                                                 >  2:6892868/1‑150 (MQ=255)
        ccTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAAc                                                                                                                                              <  2:4893753/150‑1 (MQ=255)
         cTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCCCGCCCCCGGGAAAca                                                                                                                                             >  1:4613768/1‑149 (MQ=255)
           accACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCAGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGc                                                                                                                                           <  2:1792979/149‑1 (MQ=255)
            ccACGCCAGCCTCGAACGGGCCGCCCTGGCGGCCAGCGACGGCGCCAGCTAGGTGGCCTTGGGCCGCTTCTTCACTTGGGTCCCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGcc                                                                                                                                          <  2:2236324/150‑1 (MQ=255)
                  cAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCGCGAGCAGGCCCGCGCCCCCGGGAAACGGCCGAgcgc                                                                                                                                    >  1:353508/1‑150 (MQ=255)
                    gCCTCGAACTGGCCGCCCAGGCTGCCAGCGACGGCGCCAGCTACGTGACCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCg                                                                                                                                  <  2:881560/150‑1 (MQ=255)
                       tCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGa                                                                                                                               >  2:3392600/1‑150 (MQ=255)
                        ttcgcTGGCAGCCTGGGCGGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGag                                                                                                                              >  2:2999580/5‑149 (MQ=255)
                            ccggCCGCCCCGGCGGCCAGCGACGGCGCCAGCTCCGTGGCCTTCGGCCGCTTCTTCAATTGTGTCACCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                                          <  2:3547916/148‑1 (MQ=255)
                            cTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCACGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGt                                                                                                                          <  2:2373318/150‑1 (MQ=255)
                             tGGCCGCCCCGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTCCAATTCTGCCACCACGCCGGGCGCGCCCGCGGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTg                                                                                                                         <  1:2361253/150‑1 (MQ=255)
                              ggCCGCCCTGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCGGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTgg                                                                                                                        <  2:1970443/150‑1 (MQ=255)
                              ggCCGCCCCGGCTGCCCGCGCCGGCGCCCGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGCCACCACGCGGGGCGCGCCCGCTGCCACCGGCGGCCGGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTgg                                                                                                                        <  2:1418884/150‑1 (MQ=255)
                              ggCCGCCCCGGCTGCCCGCGACGGCGCCCGCTACGTGGCCTTTGGCCGCGCCTACAATTCTGTCCCCACGCCGGGCGCGCCCGCGGCCACCGCCGGCCCGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTgg                                                                                                                        <  1:1765985/150‑1 (MQ=255)
                               tcgaatcAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGGCACCAAGCCGGGCGCGCCCGCGGCCAACGTCGGCCAGCGCGAGCAGGCCCGCGCCCAGGTGCAACAGCCGATCGCCGGGATCGGTGGc                                                                                                                       >  2:45846/7‑150 (MQ=255)
                                 cgccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCCCGCCCACGTGAAAAAGCCGATCGCCGGGCTCGGTGGCAt                                                                                                                     >  1:5604871/1‑150 (MQ=255)
                                  gccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGGCGGCGCGAGCAGGCCCGCGCCCCCGGGACACTGCCCAGCGCCGGGAGCGGTGGcctg                                                                                                                    >  1:3181551/1‑147 (MQ=255)
                                  gccCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTGCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCCGGTGAAACTGCCGATCGCCGTGCTCGGGGGCCTc                                                                                                                    >  2:4345264/1‑150 (MQ=255)
                                    ccAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCAc                                                                                                                  <  2:1040672/150‑1 (MQ=255)
                                       ggCTGCCAGCGACGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCGGCCACCGCCGGCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCt                                                                                                               <  2:6226200/150‑1 (MQ=255)
                                        gCTGCCCGCGCTGGCGCCTGCTTGGTGGCCTTGGGCCGCTTCTTAAATTCTGCCCCCACGCCGGGCGCGCCCGCCCCCCCCGGCGGCCGGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTc                                                                                                              <  1:2264921/150‑1 (MQ=255)
                                          tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGGGATCGGTGGCATCACCCACGa                                                                                                            >  1:2977106/1‑150 (MQ=255)
                                          tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGGCGGCCCGCGCCCCCGGGGAACGGCCGCGCGCCGGGAGCGGGGGCGGCACCCCCGa                                                                                                            >  1:4564127/1‑150 (MQ=255)
                                          tGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGACACCACGCCGGGCGCGCCCGCCGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGa                                                                                                            <  2:1306677/150‑1 (MQ=255)
                                            ccAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGCTCGGTGGCCTCACCCTCGACa                                                                                                          >  2:4052967/1‑150 (MQ=255)
                                              aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAc                                                                                                        >  2:7284254/1‑150 (MQ=255)
                                              aGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCACGCCCGCGCCCAGGTGAAACAGCCGAGCGCCGTGATAGGTGGCATCACCCACCACAAc                                                                                                        >  1:851940/1‑150 (MQ=255)
                                                cgtcggCGCCAGCTCGGTGGCCTTTGGCGGATTCTTCACTTGTGCCACAACGCGGGGCGCGCCCGCCCCCACCGTCGGCCTGCTCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAAcgc                                                                                                      <  1:5201738/146‑1 (MQ=255)
                                                      cgcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGcccc                                                                                                >  2:3706405/1‑150 (MQ=255)
                                                       gcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCa                                                                                               <  2:3643308/150‑1 (MQ=255)
                                                       gcCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACCCCGCCCCa                                                                                               >  2:2515050/1‑150 (MQ=255)
                                                          aGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCCCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACTg                                                                                            >  1:6677981/1‑150 (MQ=255)
                                                                ttGGACTTTGGCCTATTCTTCTATTCTGTCACCCCGCGGGGGGCGCCCGCTGCCACCGTCGGCCTGCTGGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGcc                                                                                      <  2:99207/149‑1 (MQ=255)
                                                                gTGGCTTTTGGCCGATTCTTCAATTCTGTCACCCAGCCGGGCGCGCCCGCTGCCCACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGcc                                                                                      <  1:7325089/150‑1 (MQ=255)
                                                                 tGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCCCTGGTCGccc                                                                                     >  2:1458890/1‑150 (MQ=255)
                                                                     cTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATAGCCGTGATCGGTGGCATCACCCTCCACAACGCCGCCCCACTGGTCGCCCACgg                                                                                 >  2:2712613/1‑150 (MQ=255)
                                                                     cTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCAGGGGCAACAGCCGATCGCCGGGCTCGGGGGCCGCACCCTCGCCCACGCCCCCCCCCAGGTCGccccccg                                                                                 >  1:4542798/1‑146 (MQ=255)
                                                                           ccGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGAGCAGGCCCGCGCCCCCGTGAAACTGCCGATAGCCGTGAACGGTGGCAACACCCACGACAACCCCCCCCCACACGTCGCCCACCGCGCCGa                                                                           >  1:6430657/1‑150 (MQ=255)
                                                                             gCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGAcc                                                                         >  2:6701643/1‑150 (MQ=255)
                                                                               ttcttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCCGGCCCGCGCCCACGTGAAACTGCCGATCGCCGGGATCGGTGGCCTCACCCTCGACAACGCCCCCCCCCTCGTCGCCCACCGCGccccccgg                                                                       >  1:4095463/1‑144 (MQ=255)
                                                                               ttcttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctg                                                                       >  2:6085739/1‑150 (MQ=255)
                                                                                tcttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgc                                                                      >  2:4762420/1‑150 (MQ=255)
                                                                                tcttcAATTATGTCACCCAGGCGGGCGCGCCCGCTGCCAACGTCGGCGTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgc                                                                      <  2:7069247/150‑1 (MQ=255)
                                                                                 cttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGcc                                                                             >  2:6199560/1‑142 (MQ=255)
                                                                                 cttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGcc                                                                             <  1:6199560/142‑1 (MQ=255)
                                                                                 cttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgct                                                                     <  2:3495413/150‑1 (MQ=255)
                                                                                  ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACctgctg                                                                    >  2:5622709/1‑150 (MQ=255)
                                                                                  ttcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCGCGCGCAGGCCCCCGCCCAGGGGAAACTGCCGCTCGCCGGGATCGGTGGCATCACCCACGACAACCCCCCCCCACAGGTCGCCCACCGCGCCCACcagcag                                                                    >  1:3927252/1‑148 (MQ=255)
                                                                                   tcAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGGTCGAGCAGGCCCCCGGCCAGGTGAAACAGCCGATCGCCGTGATCGGTGGCTTCACCCTCGTCACCGCCGCCCCCCGGGGCGCCCACGGCGCCGGCCTGCggg                                                                   >  1:2642902/1‑150 (MQ=255)
                                                                                     aaTACTGTCACCACGCCGGGCGCGCCCGCTGCAAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCg                                                                 <  2:7428421/150‑1 (MQ=255)
                                                                                      aTTCTGCCACCCCGCCGGGCGCGCCCGCTGCCACCGGCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCgg                                                                <  1:1105265/150‑1 (MQ=255)
                                                                                        tCTGTCACCCCGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGCCCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTg                                                              <  2:211226/150‑1 (MQ=255)
                                                                                        tCTGTCACCACGCCGGGCGCGCCCGCTGCCACCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTg                                                              <  2:4158753/150‑1 (MQ=255)
                                                                                              cccACGCCGGGCGCGCCCGCGGCCCCCGGCGGCCGGCGCGCGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCCCTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                                        <  2:2880103/149‑1 (MQ=255)
                                                                                              aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCCCCCCACAGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCAc                                                        >  1:4951994/1‑150 (MQ=255)
                                                                                              aCCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCACGTGAAACTGCCGATCGCCCTGATCGGTGGCATCACCCCCGACAACGCCGCCCCACAGGTCGCCCACCGCGCCGCCCTGCTGGCGGTGATCCAc                                                        >  1:2036779/1‑150 (MQ=255)
                                                                                               ccACGCCGGGCGCGCCCGCTGCCACCGGCGGCCTGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACg                                                       <  1:6132797/150‑1 (MQ=255)
                                                                                               ccAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCGGCTCGAGCAGGCCCGCGCCCCGGTGAAACATCCGATCGCCGTGATCGGTGGCCTCACCCACGACAACGCCGCCCCACAGGTCGCCCACCGCCCCGACCAGCTGGCGGTGATCCAcc                                                       >  1:5772119/1‑149 (MQ=255)
                                                                                                ccacccGGCCGCCCCCCCTGCCCCCGTCGGCCGGCTCGCCCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                                      <  1:3436875/146‑1 (MQ=255)
                                                                                                cAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACgg                                                      >  1:4309185/1‑150 (MQ=255)
                                                                                                       ggCGCGCCCGCCGCCAACGCCGGCCTGCGGTAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTc                                               <  2:67109/150‑1 (MQ=255)
                                                                                                        gcgcgcCGGCGCCCCGCGGCGGCCGGCTCGCCCCGCCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCg                                              <  1:3807929/150‑1 (MQ=255)
                                                                                                            gcCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGc                                          <  2:2437686/150‑1 (MQ=255)
                                                                                                             cccGCTCCAACCGTCGGCCGGCTCGCGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGcc                                         <  2:6983961/150‑1 (MQ=255)
                                                                                                              ccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCg                                        >  2:5269697/1‑150 (MQ=255)
                                                                                                              ccGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCg                                        >  2:2897828/1‑150 (MQ=255)
                                                                                                               cGCTGCCAACGCCGGCCGGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGa                                       <  2:6972463/150‑1 (MQ=255)
                                                                                                               cGCTCCCAACGTCGCCCTGATCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGa                                       <  2:3875124/150‑1 (MQ=255)
                                                                                                               cGCTCCCAACGACCGCCTGCTCGACCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGa                                       <  2:162445/150‑1 (MQ=255)
                                                                                                                    acAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCACGACAACGCCGCCCCACAGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcg                                  >  1:3688343/2‑150 (MQ=255)
                                                                                                                     cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCGTCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc                                 >  2:3912942/1‑150 (MQ=255)
                                                                                                                     cAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgc                                 >  2:549050/1‑150 (MQ=255)
                                                                                                                      aaCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcg                                >  2:3739485/1‑150 (MQ=255)
                                                                                                                      aaCGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAgcgcg                                >  1:5874504/1‑150 (MQ=255)
                                                                                                                        cGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCAGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCa                              <  2:5366555/150‑1 (MQ=255)
                                                                                                                           cGGCCTGCGCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGa                           <  2:5935459/150‑1 (MQ=255)
                                                                                                                            ggCCTGCTCGAGCCGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGaa                          <  1:3072435/150‑1 (MQ=255)
                                                                                                                             gCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAg                         <  2:465089/150‑1 (MQ=255)
                                                                                                                               cTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGc                                                                                       <  1:1493209/86‑1 (MQ=255)
                                                                                                                                tGCTCGAGCAGGTCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCa                      >  1:6867119/1‑150 (MQ=255)
                                                                                                                                  cTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAcc                    >  1:2129306/1‑150 (MQ=255)
                                                                                                                                   tCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCAccc                   >  2:2562357/1‑150 (MQ=255)
                                                                                                                                     gAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCg                                              >  1:120154/1‑121 (MQ=255)
                                                                                                                                     gAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgc                 >  1:4773766/1‑150 (MQ=255)
                                                                                                                                     gAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgc                 >  1:3614894/1‑150 (MQ=255)
                                                                                                                                     gAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgc                 >  1:6468647/1‑150 (MQ=255)
                                                                                                                                      aGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACccgcc                <  2:5754995/150‑1 (MQ=255)
                                                                                                                                        cAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCcgccgc              >  1:5362916/1‑150 (MQ=255)
                                                                                                                                           gCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcc           >  1:6033853/1‑150 (MQ=255)
                                                                                                                                           gCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcc           >  2:3194730/1‑150 (MQ=255)
                                                                                                                                            cccGCGCCCCGGTGAAACGGCCGCTCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgccc          <  1:1636357/150‑1 (MQ=255)
                                                                                                                                             ccGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcccg         >  1:5068505/1‑150 (MQ=255)
                                                                                                                                             ccGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcccg         >  2:1156925/1‑150 (MQ=255)
                                                                                                                                             ccGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcccg         >  1:2802857/1‑150 (MQ=255)
                                                                                                                                              cgcgCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcccgc        <  2:6367121/150‑1 (MQ=255)
                                                                                                                                               gcgcCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGccgcgcccgcg       <  1:7280871/150‑1 (MQ=255)
                                                                                                                                                 gcCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGCcgcgcccgcgcc     <  1:3239495/150‑1 (MQ=255)
                                                                                                                                                    cAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGCCGCGCCCGCGCCTTc  <  1:2495674/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGCTCCACCTGCCACGCCAGCCTCGAACTGGCCGCCCAGGCTGCCAGCGAAGGCGCCAGCTACGTGGCCTTTGGCCGCTTCTTCAATTCTGTCACCAAGCCGGGCGCGCCCGCTGCCAACGTCGGCCTGCTCGAGCAGGCCCGCGCCCAGGTGAAACTGCCGATCGCCGTGATCGGTGGCATCACCCTCGACAACGCCGCCCCACTGGTCGCCCACGGCGCCGACCTGCTGGCGGTGATCCACGGCTTGTTCGGTGCCGACAGCGCGCAGGAAGTCACCCGCCGCGCCCGCGCCTTC  >  NC_002947/5451158‑5451455

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