Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 555,771 C→A 5.3% H29N (CAC→AAC)  rplN → 50S ribosomal protein L14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947555,7710CA 5.3% 306.5 / 4.7 95H29N (CAC→AAC) rplN50S ribosomal protein L14
Reads supporting (aligned to +/- strand):  ref base C (45/45);  new base A (2/3);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

TGGGCACTGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCCATCAT  >  NC_002947/555622‑555916
                                                                                                                                                     |                                                                                                                                                 
tGGGCACTGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGc                                                                                                                                                   >  2:1009926/1‑150 (MQ=255)
  ggCACTGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCAc                                                                                                                                                 <  1:2083768/150‑1 (MQ=255)
        ggTTGAAGTCCTCGAACGCGCTGTTTAATTCTAAGGGCTAAGGTTTGGAGAAATTTTATGATTTATACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGt                                                                                                                                           <  2:454147/150‑1 (MQ=255)
        ggTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGt                                                                                                                                           <  1:52302/150‑1 (MQ=255)
        ggTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTGCGAAccgccga                                                                                                                                           <  2:1950128/150‑2 (MQ=255)
         gTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGtt                                                                                                                                          >  1:887359/1‑150 (MQ=255)
         gTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGtt                                                                                                                                          <  2:2264097/150‑1 (MQ=255)
             aaGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGc                                                                                                                                      <  2:1759757/150‑1 (MQ=255)
               gTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCg                                                                                                                                    <  1:653012/150‑1 (MQ=255)
                tCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTGCGAAccgccga                                                                                                                                           <  2:1785075/142‑2 (MQ=255)
                tCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTGCGAAccgccga                                                                                                                                           >  1:1785075/1‑141 (MQ=255)
                  cTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTa                                                                                                                                 >  1:250458/1‑150 (MQ=255)
                   tCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTAt                                                                                                                                >  2:952101/1‑150 (MQ=255)
                      aaCGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCCCACCGCCGTTACGCCGGTATCgg                                                                                                                             >  1:418334/1‑150 (MQ=255)
                       aCGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGt                                                                                                                            <  1:952101/150‑1 (MQ=255)
                         gcgcTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGa                                                                                                                          <  2:1678753/150‑1 (MQ=255)
                             tGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGAcatc                                                                                                                      >  1:243700/1‑150 (MQ=255)
                               ttGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGAcatcat                                                                                                                    <  1:1479990/150‑1 (MQ=255)
                               ttGAAGTCTAACGGCTAGGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTCGCCGATCACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTCCGCCGGTACCGGTGAcctcat                                                                                                                    <  2:1717/150‑1 (MQ=255)
                                 gAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACAtcatca                                                                                                                  >  2:1957492/1‑150 (MQ=255)
                                   aGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCaaa                                                                                                                <  2:250458/150‑1 (MQ=255)
                                    gTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAg                                                                                                               <  1:1957492/150‑1 (MQ=255)
                                     tCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTGCGAAccgccga                                                                                                                                           <  2:1545637/121‑2 (MQ=255)
                                     tCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTGCGAAccgccga                                                                                                                                           >  1:1545637/1‑120 (MQ=255)
                                       taagggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACGGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCg                                                                                                                                    <  2:1778854/126‑1 (MQ=255)
                                       taagggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACGGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCg                                                                                                                                    >  1:1778854/1‑126 (MQ=255)
                                        aagggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAAc                                                                                                           <  1:2010308/150‑1 (MQ=255)
                                         agggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAAcc                                                                                                          >  1:1307926/1‑150 (MQ=255)
                                         agggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAAcc                                                                                                          <  1:276427/150‑1 (MQ=255)
                                         agggataagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAAcc                                                                                                          <  1:1679628/150‑1 (MQ=255)
                                           ggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGt                                                                                                        >  1:523922/1‑150 (MQ=255)
                                           ggctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGt                                                                                                        >  1:1839700/1‑150 (MQ=255)
                                            gctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAg                                                                                                               >  2:237947/1‑142 (MQ=255)
                                            gctaagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAg                                                                                                               <  1:237947/142‑1 (MQ=255)
                                              taagggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCaa                                                                                                     <  2:1415757/150‑1 (MQ=255)
                                                 gggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGAcatca                                                                                                                     >  1:301171/1‑131 (MQ=255)
                                                 gggTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGAcatca                                                                                                                     <  2:301171/131‑1 (MQ=255)
                                                   gTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAg                                                                                                <  2:1307926/150‑1 (MQ=255)
                                                       gagaAATTTTATGATTCAGACTCAATCCATGCTAGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGa                                                                                                                          >  1:1172771/1‑120 (MQ=255)
                                                       gagaAATTTTATGATTCAGACTCAATCCATGCTAGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGa                                                                                                                          <  2:1172771/120‑1 (MQ=255)
                                                        agaAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAAtt                                                                                           >  1:334203/1‑150 (MQ=255)
                                                            aTTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCcgc                                                                                       >  2:141795/1‑150 (MQ=255)
                                                              tttATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCg                                                                                                                                    >  1:1474832/1‑103 (MQ=255)
                                                              tttATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCg                                                                                                                                    <  2:1474832/103‑1 (MQ=255)
                                                               ttATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCcgcgc                                                                                     >  2:1146862/1‑149 (MQ=255)
                                                               ttATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCcgcgc                                                                                     <  1:1146862/149‑1 (MQ=255)
                                                                   gATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTaa                                                                                <  1:670486/150‑1 (MQ=255)
                                                                    aTTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAg                                                                               <  1:141795/150‑1 (MQ=255)
                                                                      tCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGt                                                                             <  2:887359/150‑1 (MQ=255)
                                                                      tCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGt                                                                             <  2:523922/150‑1 (MQ=255)
                                                                        aGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGt                                                                                                                            >  1:1887186/1‑101 (MQ=255)
                                                                        aGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGt                                                                                                                            <  2:1887186/101‑1 (MQ=255)
                                                                        aGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTg                                                                                                                           >  2:335615/1‑102 (MQ=255)
                                                                        aGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTg                                                                                                                           <  1:335615/102‑1 (MQ=255)
                                                                        aGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCa                                                                           >  1:79832/1‑150 (MQ=255)
                                                                        aGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCa                                                                           >  2:1970950/1‑150 (MQ=255)
                                                                         gACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCaa                                                                          >  1:838273/1‑150 (MQ=255)
                                                                            tCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGACAAGGAAGCAATTCCGCGCGGTAAGGTCaaaaa                                                                       >  1:644178/1‑150 (MQ=255)
                                                                             cAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCaaaaaa                                                                      <  1:426502/150‑1 (MQ=255)
                                                                              aaTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAg                                                                     >  1:683430/1‑150 (MQ=255)
                                                                              aaTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAg                                                                     >  2:687968/1‑150 (MQ=255)
                                                                               aTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAg                                                                                                    >  2:2286387/1‑118 (MQ=255)
                                                                               aTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAg                                                                                                    <  1:2286387/118‑1 (MQ=255)
                                                                               aTCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAgg                                                                    >  2:1397386/1‑150 (MQ=255)
                                                                                   aTGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCaa                                                                <  2:79832/150‑1 (MQ=255)
                                                                                    tGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAg                                                               >  2:1193022/1‑150 (MQ=255)
                                                                                      cTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGtg                                                             >  1:439999/1‑150 (MQ=255)
                                                                                       tCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGtga                                                            >  1:208033/1‑150 (MQ=255)
                                                                                        cGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGtgat                                                           <  2:2203803/150‑1 (MQ=255)
                                                                                         gATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGt                                                                             <  2:880066/131‑1 (MQ=255)
                                                                                         gATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGt                                                                             >  1:880066/1‑131 (MQ=255)
                                                                                             tgGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCACGGTCCTCGGCGGTTCGCACCCCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCg                                                      <  1:1349841/150‑1 (MQ=255)
                                                                                               gCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCt                                                    >  2:548248/1‑150 (MQ=255)
                                                                                                 cGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCCAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGt                                                              <  1:2294465/138‑1 (MQ=255)
                                                                                                 cGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGt                                                              >  2:2294465/1‑138 (MQ=255)
                                                                                                 cGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTgt                                                  >  2:1384037/1‑150 (MQ=255)
                                                                                                         gcggcgCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTcgtac                                          >  1:1030747/1‑150 (MQ=255)
                                                                                                         gcggcgCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGAAAAAAAGAAAAAAAAGGCCAAGTGATGACCGCTGTCGTcgtac                                          <  2:243700/150‑1 (MQ=255)
                                                                                                          cggcgCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAAGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTcgtacg                                         >  1:1934847/1‑150 (MQ=255)
                                                                                                             cgcTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCgtacgtac                                      <  1:687968/150‑1 (MQ=255)
                                                                                                              gcTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTAcc                                     >  1:1549071/1‑150 (MQ=255)
                                                                                                                tcgtcgCGTCATGTGAATCAAGGTACGGGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCg                                   >  2:824780/1‑150 (MQ=255)
                                                                                                                        tCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGACAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGgtg                           <  2:644178/150‑1 (MQ=255)
                                                                                                                            gtgCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTtcg                       >  2:169349/1‑150 (MQ=255)
                                                                                                                                aTCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCaa                                                                                             <  1:509840/76‑1 (MQ=255)
                                                                                                                                aTCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCaa                                                                                             >  2:509840/1‑76 (MQ=255)
                                                                                                                                aTCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTcgc                   >  1:1436412/1‑150 (MQ=255)
                                                                                                                                aTCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTcgc                   >  2:775754/1‑150 (MQ=255)
                                                                                                                                      gTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCcgc                                                                                       >  2:472358/1‑76 (MQ=255)
                                                                                                                                      gTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCcgc                                                                                       <  1:472358/76‑1 (MQ=255)
                                                                                                                                          ccggcggTTCCCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGtt         <  1:1614335/149‑1 (MQ=255)
                                                                                                                                           cggcggTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTc        >  2:820892/1‑150 (MQ=255)
                                                                                                                                            ggcggTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTcc       <  1:820892/150‑1 (MQ=255)
                                                                                                                                             gcggTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCACTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCca      <  1:1193022/150‑1 (MQ=255)
                                                                                                                                                gTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCcatca   <  2:208033/150‑1 (MQ=255)
                                                                                                                                                 ttCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCcatcat  >  1:1899109/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGGGCACTGGTTGAAGTCCTCGAACGCGCTGTTGAAGTCTAAGGGCTAAGGGTCGGAGAAATTTTATGATTCAGACTCAATCCATGCTCGATGTGGCCGATAACAGCGGCGCTCGTCGCGTCATGTGCATCAAGGTACTCGGCGGTTCGCACCGCCGTTACGCCGGTATCGGTGACATCATCAAAGTAACCGTCAAGGAAGCAATTCCGCGCGGTAAGGTCAAAAAAGGCCAAGTGATGACCGCTGTCGTCGTACGTACCCGTCACGGTGTTCGTCGCGCTGACGGTTCCATCAT  >  NC_002947/555622‑555916

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: