Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 926,438 C→G 6.7% S64T (AGT→ACT)  PP_0805 ← outer membrane efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947926,4380CG 6.7% 252.4 / 6.3 87S64T (AGT→ACT) PP_0805outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base C (54/26);  new base G (2/4);  total (56/31)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

ATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCAGGAGCCCGA  >  NC_002947/926291‑926585
                                                                                                                                                    |                                                                                                                                                   
aTACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTg                                                                                                                                                   <  1:3035710/150‑1 (MQ=255)
  aCCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGc                                                                                                                                                 >  1:1210061/1‑150 (MQ=255)
    cgccgGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGa                                                                                                                                               >  2:1634785/1‑150 (MQ=255)
     gccgGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCATGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGAc                                                                                                                                              >  1:2503315/1‑150 (MQ=255)
        ggAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCgg                                                                                                                                           <  1:1997092/150‑1 (MQ=255)
         gAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCggg                                                                                                                                          <  1:2702723/150‑1 (MQ=255)
           aaTCTGTGGGTACTACCACGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAACAGGTGCCGATGGCCTCAGCGATGCTAGGATGCCATTGCGCTGCCTGCAGCACCGCAGCGGTAAGGTCCAGCCGTTCGCCACTCGCGACCGGGGc                                                                                                                                        >  2:397937/1‑150 (MQ=255)
            aTCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGcc                                                                                                                                       >  1:1319325/1‑150 (MQ=255)
            aTCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGGTCGCCAGTGGCGAACggctga                                                                                                                                       <  2:340495/150‑5 (MQ=255)
             tCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCGTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGACTGCAGCACCGCGTCGGTAAGGCTCAGCCGTCCGCCACTGGCCAACGGCGCCa                                                                                                                                      >  2:2840873/1‑150 (MQ=255)
              cTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCAt                                                                                                                                     <  2:309173/150‑1 (MQ=255)
              cTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCAt                                                                                                                                     >  2:1827257/1‑150 (MQ=255)
               tgtgGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATg                                                                                                                                    >  1:2818110/1‑150 (MQ=255)
                    gTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGAt                                                                                                                               >  1:268357/1‑150 (MQ=255)
                       aTACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTTCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATggg                                                                                                                            <  1:370017/150‑1 (MQ=255)
                        tACCCCGCACGGGCGCCGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGGTGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATgggg                                                                                                                           >  2:1289526/1‑150 (MQ=255)
                        tACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATgggg                                                                                                                           <  1:2668771/150‑1 (MQ=255)
                        tACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATgggg                                                                                                                           >  1:899296/1‑150 (MQ=255)
                          ccccGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATgggggg                                                                                                                         >  2:2342120/1‑150 (MQ=255)
                           cccGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATggggggg                                                                                                                        >  1:2902448/1‑150 (MQ=255)
                           cccGCACGGGCGACGTTAATGCCTTCCCCCTCCACGTAGAGTCAGCCGATGGCCTCAGCGATTCACGTGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTCGCGACCGGGGCCATGCCGAtgggtggg                                                                                                                        >  2:1274334/1‑150 (MQ=255)
                            ccGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggc                                                                                                                       <  1:388603/150‑1 (MQ=255)
                               cACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCATTAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcggg                                                                                                                    <  1:188842/150‑1 (MQ=255)
                                      gACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGAt                                                                                                                               >  2:1942944/1‑132 (MQ=255)
                                      gACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGAt                                                                                                                               <  1:1942944/132‑1 (MQ=255)
                                      gACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcgg                                                                                                                     >  1:1987551/1‑142 (MQ=255)
                                      gACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcgg                                                                                                                     <  2:1987551/142‑1 (MQ=255)
                                         gTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGc                                                                                                          >  2:1549166/1‑150 (MQ=255)
                                            gATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTc                                                                                                       >  1:197865/1‑150 (MQ=255)
                                             aTGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCg                                                                                                      >  1:2056139/1‑150 (MQ=255)
                                               gCCTTCGCCGTGGTGGTAAAGGTGGCCGATGGCTTCAGGGAGTCTGGGGTGCCTTGGCGCTGCCGGAAGACCGGCTTCGGTATGGCTCAGCTTTCCGCCAGTGGCGCCCGGGCCCATCCGGTTGGGGGGGCGGGCGGGCACCGCTTCggg                                                                                                    <  2:588537/150‑1 (MQ=255)
                                                 cTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcgggcg                                                                                                                  >  2:1839582/1‑134 (MQ=255)
                                                 cTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcgggcg                                                                                                                  <  1:1839582/134‑1 (MQ=255)
                                                     cccctttttGTAAGGGTTGCGGAGGGCCTCAGCGGTGCGGGGGCCCCTTTGCGCTGCCTGCACCCCCGGGTGGGTAGGGCTCCGCCTTCCCCCATGGGCGCCGGGGGCCATGCCGGGGGGGGGGGGGGCGGGCACCGCTTCGGGCGCATc                                                                                              <  2:2503315/143‑1 (MQ=255)
                                                       ccTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATcgc                                                                                            >  1:641132/1‑150 (MQ=255)
                                                       ccTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTACCATTGCGCTGCCTGCACCACCGCGTCGGTAAGGCTCATCCGGTCGCCACTGGCGACCGGGGCCATGCCGCTGGGGGGGCGGGCGGGCACCGCTTAGGGCGCATcgc                                                                                            >  1:1107764/1‑150 (MQ=255)
                                                       ccTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCCCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATcgc                                                                                            >  2:701678/1‑150 (MQ=255)
                                                        tgtgttGTAAAGGGTTGCGGAGGGCCTCGGCGTGGCTCGGGTGCTGTGGCGCTCCCGGCCGCACCGCGTCGGGAGGCCTCGCCCTTCCCCCCCGGGCGCCCGGGCCCATGCCGAGGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATcgc                                                                                            <  2:2878394/146‑1 (MQ=255)
                                                        cTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGTCGGGCGGGCACCGCTTCGGGCGCATCGCGt                                                                                          <  2:1210061/150‑1 (MQ=255)
                                                         tGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATcgcga                                                                                          >  1:2878394/1‑149 (MQ=255)
                                                           cTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTag                                                                                        >  2:984661/1‑150 (MQ=255)
                                                              gTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTagcag                                                                                     >  1:2295199/1‑150 (MQ=255)
                                                               tAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATg                                                                                                                              <  1:2721639/108‑1 (MQ=255)
                                                               tAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATg                                                                                                                              >  2:2721639/1‑108 (MQ=255)
                                                               tAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTagcagc                                                                                    >  1:2470140/1‑150 (MQ=255)
                                                               tAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTagcagc                                                                                    >  1:2274468/1‑150 (MQ=255)
                                                                aaagggTTGCGGAGGGCCTGGGCGAGGCTGGGGTGCCGTTGCGGGGCCGGCAGCAGCGCGTCGGTAAGGCTCAGCCGTTGGCCACTGGCGACCGGGGCCAGGCGGAGGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTagcagc                                                                                    <  2:1348562/146‑1 (MQ=255)
                                                                aaaGGTTGCCGATGGCCTCAGCGATGCTCGGGTCCCATTGCGGTGCCTGCAGCCCGGCGTCGGTAAGGCTCAGCCGTTCCCCACTGGCGCCGGGGGCCATGCCGAGGGGGGGGGGGGGGGGCACCGCTTCGGGCGCATCGCGTAGCAGCg                                                                                   <  1:984661/150‑1 (MQ=255)
                                                                 aaGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcggg                                                                                                                    <  1:500359/116‑1 (MQ=255)
                                                                 aaGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGgggcggg                                                                                                                    >  2:500359/1‑116 (MQ=255)
                                                                  aGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGGACCGCTTCGGGCGCATCGCGTAGCAGCGag                                                                                 >  2:1521835/1‑149 (MQ=255)
                                                                  aGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGGTCGCCAGTGGCGAACggctga                                                                                                                                       >  1:850579/1‑92 (MQ=38)
                                                                  aGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGGTCGCCAGTGGCGAACggctga                                                                                                                                       <  2:850579/96‑5 (MQ=38)
                                                                     gTTGCCGATGGCCTCAGCGATGCTCTGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCAc                                                                                                                                                     <  1:373676/80‑1 (MQ=255)
                                                                     gTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCAc                                                                                                                                                     >  2:373676/1‑80 (MQ=255)
                                                                     gTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACgg                                                                               >  1:2199608/1‑150 (MQ=255)
                                                                        gCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCGGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGccc                                                                            <  1:2310204/150‑1 (MQ=255)
                                                                                ccTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGgc                                                                    <  1:2342120/150‑1 (MQ=255)
                                                                                 cTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGAc                                                                                 >  2:1552210/1‑136 (MQ=255)
                                                                                 cTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGAc                                                                                 <  1:1552210/136‑1 (MQ=255)
                                                                                 cTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGgcg                                                                   >  2:2535189/1‑150 (MQ=255)
                                                                                 cTCAGCGATGCTCGGGTGCCATTGCGCTGAGCGCAGCACGGCGTCGGTGGTGCGAAGTCGAGCGCTACTTGGAGCCGGGGCCGTGTCGATTGGTGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGgcg                                                                   >  2:349789/1‑150 (MQ=255)
                                                                                  tCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGAGCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCGCAGCGTCGCAGCGACGGCCCGATCTGgcgc                                                                  >  2:1282155/1‑150 (MQ=255)
                                                                                  tCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGgcgc                                                                  >  1:2880687/1‑150 (MQ=255)
                                                                                        aTGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTgg                                                            >  2:1985612/1‑150 (MQ=255)
                                                                                              ggggGCCTTTGCGCTCCCTGCAGCCCCGCGTCGGAAAGGCCCGGCCTTTGCCCCCGGGGGACCGGGGCCAGGCCGGGGGGGGGGGGGGGGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCt                                                      <  1:716688/150‑1 (MQ=255)
                                                                                                gTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCtat                                                    >  2:1726688/1‑150 (MQ=255)
                                                                                                    cATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATatca                                                >  1:1005771/1‑150 (MQ=255)
                                                                                                    cATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAACCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATatca                                                >  1:2488045/1‑150 (MQ=255)
                                                                                                            tgcctgcAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCt                                                      <  1:1010889/136‑1 (MQ=255)
                                                                                                            tgcctgcAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCt                                                      >  2:1010889/1‑136 (MQ=255)
                                                                                                                tgcAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTgg                                    >  1:2818571/1‑150 (MQ=255)
                                                                                                                    gcaCCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCaa                                >  2:2286282/1‑150 (MQ=255)
                                                                                                                         gcgTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGGAGGTGGGGGTCTATATCAGCGGCTGCCTGGGCAACCggg                           >  2:1430805/1‑150 (MQ=255)
                                                                                                                              ggTAAGGCTCAGCGGTTCGCCACTGGCGACGGGGGCAATGCCGAGGGGGGGGGGGGGGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGaa                      <  2:197865/150‑1 (MQ=255)
                                                                                                                               gaaaGGCTCCGCCGTCCGCCCCGGGCGCCGGGGGCTAGGCGGGGGGGGGGGGGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAg                     <  2:2056139/148‑1 (MQ=255)
                                                                                                                               gTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCCCCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAg                     >  2:1142848/1‑150 (MQ=255)
                                                                                                                               gTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGGCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAg                     >  1:2956430/1‑150 (MQ=255)
                                                                                                                                tAAGGCGTAGCCGTTCGCCAGGGGCTCCGGGGGCAATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGt                    <  1:927687/150‑1 (MQ=255)
                                                                                                                                  gggCTCAGCCGTCCGCCACGGGCGCCGGGGGCCAGGCGGATGGGGGGGGGGGGGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTc                  <  2:2488045/149‑1 (MQ=255)
                                                                                                                                  aGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTc                  >  2:1092237/1‑150 (MQ=255)
                                                                                                                                      tCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCTCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAggg              >  1:588546/1‑150 (MQ=255)
                                                                                                                                        aGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGcc            >  1:1504030/1‑150 (MQ=255)
                                                                                                                                        aGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGcc            >  2:2768159/1‑150 (MQ=255)
                                                                                                                                         cccggtCGCCAGTGGCGAACGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGGCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCa           >  1:2266127/6‑150 (MQ=255)
                                                                                                                                           cGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCAgg         >  1:1181322/1‑150 (MQ=255)
                                                                                                                                            gTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCAGGa        >  1:1410613/1‑150 (MQ=255)
                                                                                                                                              tCGCCACTGGGCACGGGGGCCCGGCGGAGGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCAGGAGc      <  2:641132/150‑1 (MQ=255)
                                                                                                                                                  cACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCAGGAGCCCGa  >  2:1034794/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATACCGCCGGAAATCTGTGGGTAATACCCCGCACGGGCGACGTTGATGCCTTCGCCCTGCTTGTAAAGGTTGCCGATGGCCTCAGCGATGCTCGGGTGCCATTGCGCTGCCTGCAGCACCGCGTCGGTAAGGCTCAGCCGTTCGCCACTGGCGACCGGGGCCATGCCGATGGGGGGGCGGGCGGGCACCGCTTCGGGCGCATCGCGTAGCAGCGACGGCCCGATCTGGCGCAGGTGGGGGTCTATATCATCGGCTGCCTGGGCAACCGGGTGGAAGTTCAGGGCCAGGAGCCCGA  >  NC_002947/926291‑926585

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