Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,568,736 T→A 5.5% Q285H (CAA→CAT PP_2252 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,568,7360TA 5.5% 283.9 / 5.5 91Q285H (CAA→CATPP_2252hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (50/35);  new base A (2/3);  total (53/38)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

ACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTC  >  NC_002947/2568596‑2568885
                                                                                                                                            |                                                                                                                                                     
aCCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGAcc                                                                                                                                              >  2:1819559/1‑150 (MQ=255)
     ccGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTg                                                                                                                                         >  2:138062/1‑150 (MQ=255)
        cggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCgg                                                                                                                                      <  1:1819559/150‑1 (MQ=255)
        cggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACACCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCCGATGGTCGAACAACTGCgg                                                                                                                                      <  2:1108575/150‑1 (MQ=255)
            gcACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGACCACCCGCAGATGGTCGAACAACTGCGggtga                                                                                                                                  <  2:1612149/150‑2 (MQ=255)
                       cAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATg                                                                                                                       <  2:655960/150‑1 (MQ=255)
                        aaGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                      <  2:1954235/150‑1 (MQ=255)
                         aGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATggg                                                                                                                     <  1:514316/150‑1 (MQ=255)
                           cTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtc                                                                                                                   >  1:2325607/1‑150 (MQ=255)
                           cTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtc                                                                                                                   <  1:358707/150‑1 (MQ=255)
                              tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatc                                                                                                                >  1:326091/1‑150 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                      >  1:534129/1‑142 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                      <  2:534129/142‑1 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                              >  2:1558076/1‑150 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                              >  2:810427/1‑150 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                              >  1:1190804/1‑150 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                              >  1:1391230/1‑150 (MQ=255)
                                 tCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcc                                                                                                             >  1:253784/1‑150 (MQ=255)
                                    aCCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCg                                                                                                          >  2:1690351/1‑150 (MQ=255)
                                      cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGCTCACCCCCAGTTGTTCGACCATCTGCGGGTGTGCACGGTCATGGGTCATCACCCCGTa                                                                                                        >  1:208736/1‑150 (MQ=255)
                                      cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCCCCCGCAGTTGTTCGa                                                                                                                                                <  2:1473440/110‑1 (MQ=255)
                                      cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGa                                                                                                                                                >  1:1473440/1‑110 (MQ=255)
                                      cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                      <  2:2176858/136‑1 (MQ=255)
                                      cTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATgg                                                                                                                      >  1:2176858/1‑136 (MQ=255)
                                       tGCCCCGACAGGTAATAGACCCCCGCCATCGCCACGCTCTGGCGGTTCTCGATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGCGTGGTCACGGTCATGGGTCATCACCCCGTaa                                                                                                       >  2:147591/1‑150 (MQ=255)
                                        gCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAAc                                                                                                      <  2:620286/150‑1 (MQ=255)
                                           ccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTc                                                                                                                                                  >  1:2050245/1‑103 (MQ=255)
                                           ccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTc                                                                                                                                                  <  2:2050245/103‑1 (MQ=255)
                                           ccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCgg                                                                                                   >  2:1969720/1‑150 (MQ=255)
                                           ccGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCgg                                                                                                   >  1:14992/1‑150 (MQ=255)
                                               cAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCATGTGAGATAGGAAGCGCAGTTGTTCGACCATCTGCGGGTGGTAACGGTGGTGGGGCATCAGCGCGTAACCGGCTGc                                                                                               >  1:1022249/1‑150 (MQ=255)
                                               cAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                             >  1:400324/1‑120 (MQ=255)
                                               cAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACg                                                                                                                             <  2:400324/120‑1 (MQ=255)
                                                 ggTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                             <  2:326091/150‑1 (MQ=255)
                                                    aaTAGCCCTCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCATCCTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCaa                                                                                          <  1:637235/150‑1 (MQ=255)
                                                    aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCaa                                                                                          <  1:1558076/150‑1 (MQ=255)
                                                     aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                           >  1:4463/1‑132 (MQ=255)
                                                     aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                           <  2:4463/132‑1 (MQ=255)
                                                      tAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATc                                                                                        <  1:2226482/150‑1 (MQ=255)
                                                         aCCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaa                                                                                     >  1:1097378/1‑150 (MQ=255)
                                                          cccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaaa                                                                                    <  1:1988991/150‑1 (MQ=255)
                                                           ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                           <  1:1585686/126‑1 (MQ=255)
                                                           ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAcccc                                                                                                           >  2:1585686/1‑126 (MQ=255)
                                                           ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGATGGTCGAACAACTGCgg                                                                                                                                      <  2:339368/99‑1 (MQ=38)
                                                           ccccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGATGGTCGAACAACTGCgg                                                                                                                                      >  1:339368/1‑99 (MQ=38)
                                                           ccccGCCATGGCCCCGCTCTGGCGGCACTGCATAACCTTCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAat                                                                                   >  1:1171408/1‑150 (MQ=255)
                                                            cccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCATGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAatc                                                                                  >  2:593430/1‑150 (MQ=255)
                                                             ccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGggtgg                                                                                                                                  >  1:733323/1‑101 (MQ=255)
                                                             ccGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGggtgg                                                                                                                                  <  2:733323/101‑1 (MQ=255)
                                                                  aTGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGcc                                                                            <  1:593430/150‑1 (MQ=255)
                                                                    ggCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGcccc                                                                          >  1:309890/1‑150 (MQ=255)
                                                                      ccccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCCTCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTg                                                                        >  1:633774/1‑150 (MQ=255)
                                                                          gCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCAc                                                                    >  1:1275135/1‑150 (MQ=255)
                                                                                gCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCAGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGcc                                                              <  2:1391230/150‑1 (MQ=255)
                                                                                    ttCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTGCCCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACt                                                          >  1:24780/1‑150 (MQ=255)
                                                                                      ctcCATCATCTCCTGCCATTCCTGGTTCAACGCCAAGTTCAGGATCACCCGCAGTTGTTCGACCAGCTGCGTGTGGTCACGGTAATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCATCTTGCCCACTGc                                                        <  2:208736/150‑1 (MQ=255)
                                                                                          aTCAACTCCTGCCATGCCTGGTTCACCGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAATTCATAGCACCTGCCACCCTGCCCACTCCCGaa                                                    >  2:271411/1‑149 (MQ=255)
                                                                                           tCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                                   >  1:522372/1‑150 (MQ=255)
                                                                                           ccAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATg                                                   >  2:977272/2‑150 (MQ=255)
                                                                                                   tGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGcc                                           >  1:398960/1‑149 (MQ=255)
                                                                                                    gCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCAc                                          >  2:2093511/1‑150 (MQ=255)
                                                                                                     ccATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACt                                         >  1:1545122/1‑150 (MQ=255)
                                                                                                        ttCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAatcatc                                                                               >  2:822222/1‑109 (MQ=255)
                                                                                                        ttCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAatcatc                                                                               <  1:822222/109‑1 (MQ=255)
                                                                                                          ccTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGa                                    <  2:2325607/150‑1 (MQ=255)
                                                                                                            tGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATc                                  >  2:1639892/1‑150 (MQ=255)
                                                                                                            tGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATc                                  >  2:2166931/1‑150 (MQ=255)
                                                                                                               ttCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCgg                                                                                                   >  1:684456/1‑82 (MQ=255)
                                                                                                               ttCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCgg                                                                                                   <  2:684456/82‑1 (MQ=255)
                                                                                                                tCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCcgc                              <  2:708268/150‑1 (MQ=255)
                                                                                                                  aaCGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                            <  1:1571659/150‑1 (MQ=255)
                                                                                                                  aaCGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                            <  2:1693652/150‑1 (MQ=255)
                                                                                                                  aaCGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGt                            >  1:1710465/1‑150 (MQ=255)
                                                                                                                     gCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTggg                         <  1:1948026/150‑1 (MQ=255)
                                                                                                                          gTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                             <  1:1831814/77‑1 (MQ=255)
                                                                                                                          gTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCgg                                                                                             >  2:1831814/1‑77 (MQ=255)
                                                                                                                           ttCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAgg                   >  2:1317162/1‑150 (MQ=255)
                                                                                                                              atgATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGCCTGCGGCAATCGAACTCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTg                <  2:1022249/148‑1 (MQ=255)
                                                                                                                              aGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGGCAGGGTg                >  2:1404293/1‑150 (MQ=255)
                                                                                                                                gATCACCCGCAGGTGTTGGACGATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgc              >  2:1843892/1‑150 (MQ=255)
                                                                                                                                 aTCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg             <  2:253784/150‑1 (MQ=255)
                                                                                                                                     cccGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAccc                                                                                                            <  2:2145603/51‑1 (MQ=255)
                                                                                                                                     cccGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAccc                                                                                                            >  1:2145603/1‑51 (MQ=255)
                                                                                                                                     cccGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCaa         <  2:309890/150‑1 (MQ=255)
                                                                                                                                      ccGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAg        >  2:983257/1‑150 (MQ=255)
                                                                                                                                      ccGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAg        >  1:702567/1‑150 (MQ=255)
                                                                                                                                      ccGCAGATGGTCGAACAACTGCGGTTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAg        >  1:1721001/1‑150 (MQ=255)
                                                                                                                                       cGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGc       >  1:439963/1‑150 (MQ=255)
                                                                                                                                        gCAGTTGTTCGACCATCTGCAGGAGGACACGGTCATAGGTCATCACCCCGTAACAGGCTGCGGCAAACGAAATCATCGCCCCTACAACCTTGCCCACTGCCGATGCCACGGCACTGACGATCCCGCGTGGCGCCAGGGTGgcgcacagat      >  2:976637/1‑144 (MQ=255)
                                                                                                                                          aGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGc                                                                             >  2:2228630/1‑77 (MQ=255)
                                                                                                                                          aGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGc                                                                             <  1:2228630/77‑1 (MQ=255)
                                                                                                                                           gttgttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGTCAGGGTGGCGCCAAGCtttt   <  1:1396488/150‑1 (MQ=255)
                                                                                                                                            ttgttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc  >  2:134439/1‑150 (MQ=255)
                                                                                                                                            ttgttAGAGCATTTGCGGGTGGTCACGGACAAGGCGGCGCACCCCGTAACCGGCTGCTGTAATCGAAATCATCGCCCCTGCCACATTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTc  <  2:1171408/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
ACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTC  >  NC_002947/2568596‑2568885

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