Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,153 T→C 7.4% D461G (GAC→GGC)  dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,1530TC 7.4% 667.6 / 27.5 217D461G (GAC→GGC) dnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (111/90);  new base C (7/9);  total (118/99)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 219 total aligned reads displayed.

CTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGA  >  NC_002947/5378004‑5378300
                                                                                                                                                     |                                                                                                                                                   
ctTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGt                                                                                                                                                     >  1:2693374/1‑150 (MQ=255)
   cGTCCGACAGTCCTTAGTTCGCCATGATCACCATCGACTGAGTCTTGCCCGTCGCCTTGTCTTACGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGCCACGCCGCGCGGAGCTGGCGCAATGTCAg                                                                                                                                                  <  2:442441/150‑1 (MQ=255)
         aCAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGt                                                                                                                                            >  2:4766858/1‑150 (MQ=255)
          cAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTc                                                                                                                                           <  1:1511875/150‑1 (MQ=255)
           aGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCg                                                                                                                                          <  2:1255278/150‑1 (MQ=255)
             gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                        >  2:5030180/1‑150 (MQ=255)
              ccGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAAc                                                                                                                                       >  1:1458800/1‑150 (MQ=255)
               cGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACt                                                                                                                                      >  2:3031742/1‑150 (MQ=255)
                ggAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGAGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACtt                                                                                                                                     >  1:4088311/1‑150 (MQ=255)
                    ttGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGccc                                                                                                                                 >  1:636455/1‑150 (MQ=255)
                       gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                              >  2:960691/1‑150 (MQ=255)
                       gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                              <  2:2773813/150‑1 (MQ=255)
                       gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                              >  2:351911/1‑150 (MQ=255)
                        ccTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCg                                                                                                                             >  2:1022832/1‑150 (MQ=255)
                        ccTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCg                                                                                                                             >  1:2240812/1‑150 (MQ=255)
                           tGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACt                                                                                                                          >  1:2134212/1‑150 (MQ=255)
                            gATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACtt                                                                                                                         >  2:4081090/1‑150 (MQ=255)
                             aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                                        <  2:1519438/150‑1 (MQ=255)
                             aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                                        <  2:623984/150‑1 (MQ=255)
                               cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                                      >  1:955358/1‑150 (MQ=255)
                                  gATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTg                                                                                                                   >  2:3924795/1‑150 (MQ=255)
                                  gATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTg                                                                                                                   >  1:5448859/1‑150 (MQ=255)
                                   ctCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                                  <  1:3148498/149‑1 (MQ=255)
                                    tCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgcg                                                                                                                 >  1:195408/1‑150 (MQ=255)
                                      gACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCa                                                                                                               <  1:2377662/150‑1 (MQ=255)
                                      gACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCa                                                                                                               <  1:4077108/150‑1 (MQ=255)
                                        gtgttccTTGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGc                                                                                                             >  1:4976651/7‑150 (MQ=255)
                                         tgttccTTGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGcc                                                                                                            >  1:904020/6‑150 (MQ=255)
                                         tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                                            >  2:3675436/1‑150 (MQ=255)
                                          gAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCt                                                                                                           <  2:4679394/150‑1 (MQ=255)
                                             tCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCt                                                                                                        >  2:4078634/1‑150 (MQ=255)
                                               ttGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGccatttca                                                                                                      >  2:2424572/1‑144 (MQ=255)
                                               ttGCCGGTGGCCTTGCCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTa                                                                                                      >  2:4462894/1‑150 (MQ=255)
                                                        gCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCt                                                                                             >  2:2441368/1‑150 (MQ=255)
                                                          cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTg                                                                                                                                    >  2:3825701/1‑109 (MQ=255)
                                                          cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                           >  1:100934/1‑150 (MQ=255)
                                                          cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                           >  1:4275496/1‑150 (MQ=255)
                                                          cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                           >  1:4867909/1‑150 (MQ=255)
                                                           ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                          >  1:98178/1‑150 (MQ=255)
                                                            tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                         <  2:4981605/150‑1 (MQ=255)
                                                            tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                         <  1:2956362/150‑1 (MQ=255)
                                                              tCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGTTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgcagc                                                                                       >  1:1062923/1‑150 (MQ=255)
                                                              tCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgcagc                                                                                       >  1:442630/1‑150 (MQ=255)
                                                                tttCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCAc                                                                                     >  1:1666580/1‑150 (MQ=255)
                                                                     cACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGaa                                                                                <  2:1823170/150‑1 (MQ=255)
                                                                         tACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGt                                                                            <  2:2719184/150‑1 (MQ=255)
                                                                           cGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCa                                                                          <  1:1836426/150‑1 (MQ=255)
                                                                            gTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCAc                                                                         >  1:3796675/1‑150 (MQ=255)
                                                                             tGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                          <  2:198669/132‑1 (MQ=255)
                                                                              gCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACgg                                                                       >  1:676058/1‑150 (MQ=255)
                                                                              gCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACgg                                                                       <  2:4140806/150‑1 (MQ=255)
                                                                                acgATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGgcc                                                                     >  1:4335881/3‑149 (MQ=255)
                                                                                  gATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGct                                                                   >  2:2199384/1‑150 (MQ=255)
                                                                                  gATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGct                                                                   >  1:3928419/1‑150 (MQ=255)
                                                                                     gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                                  >  2:5303098/1‑100 (MQ=255)
                                                                                     gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                                  <  1:5303098/100‑1 (MQ=255)
                                                                                       cGTCGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGt                                                              <  1:3422753/150‑1 (MQ=255)
                                                                                        gTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGtt                                                             <  2:5406617/150‑1 (MQ=255)
                                                                                           ggCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtc                                                          <  2:3611848/150‑1 (MQ=255)
                                                                                            gCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtcg                                                         >  1:4151997/1‑150 (MQ=255)
                                                                                               tCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGtcgtcg                                                      <  2:5240399/150‑1 (MQ=255)
                                                                                               tCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTgtcttca                                                      <  2:1653256/150‑2 (MQ=255)
                                                                                                     tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                <  2:2693374/150‑1 (MQ=255)
                                                                                                     tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                >  2:4015637/1‑150 (MQ=255)
                                                                                                     tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                <  1:3534671/150‑1 (MQ=255)
                                                                                                     tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                <  2:609203/150‑1 (MQ=255)
                                                                                                      cGAAGGTCACTTCGATGTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCg                                               <  2:4352931/150‑1 (MQ=255)
                                                                                                       gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTc                                                                                                                     >  2:1052796/1‑79 (MQ=255)
                                                                                                       gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                         >  2:4874392/1‑107 (MQ=255)
                                                                                                       gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                         <  1:4874392/107‑1 (MQ=255)
                                                                                                       gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                              >  2:1192355/1‑150 (MQ=255)
                                                                                                       gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                              >  1:3326904/1‑150 (MQ=255)
                                                                                                            tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                         >  2:1312848/1‑150 (MQ=255)
                                                                                                            tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                         >  1:4489858/1‑150 (MQ=255)
                                                                                                            tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                         >  1:5370355/1‑150 (MQ=255)
                                                                                                             cACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAacac                                        >  1:4679099/1‑150 (MQ=255)
                                                                                                                ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                     >  1:3686955/1‑150 (MQ=255)
                                                                                                                ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                     >  1:4770822/1‑150 (MQ=255)
                                                                                                                ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGACCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACATTGCAGCACGTGAATGGTCACGGGGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                     >  2:3672717/1‑150 (MQ=255)
                                                                                                                 tCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                                    <  1:1031645/150‑1 (MQ=255)
                                                                                                                   gATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGa                                  <  2:3948635/150‑1 (MQ=255)
                                                                                                                    aTTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGAc                                 <  2:2134212/150‑1 (MQ=255)
                                                                                                                         cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtctt                            <  1:2019787/150‑1 (MQ=255)
                                                                                                                          ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg                           >  1:398183/1‑150 (MQ=255)
                                                                                                                          ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg                           <  2:2245278/150‑1 (MQ=255)
                                                                                                                           gCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                >  1:11832/1‑128 (MQ=255)
                                                                                                                               gccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCgg                      >  2:518669/1‑150 (MQ=255)
                                                                                                                                 cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGGCTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa                    <  2:2959103/150‑1 (MQ=255)
                                                                                                                                   cgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtg                  <  1:2412622/150‑1 (MQ=255)
                                                                                                                                      ggTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGAtcgtc               <  2:3293581/150‑2 (MQ=255)
                                                                                                                                       gTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtgg              <  2:2890598/150‑1 (MQ=255)
                                                                                                                                         gCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGgtggtg            <  1:130684/150‑1 (MQ=255)
                                                                                                                                           tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                         <  1:2905534/71‑1 (MQ=255)
                                                                                                                                           tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                         >  2:2905534/1‑71 (MQ=255)
                                                                                                                                            ggCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttc         <  1:5030180/150‑1 (MQ=255)
                                                                                                                                             gCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttct        >  1:4700344/1‑150 (MQ=255)
                                                                                                                                             gCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgc        <  1:1080889/150‑3 (MQ=255)
                                                                                                                                             gCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgc        >  1:2789669/1‑148 (MQ=255)
                                                                                                                                                 aaTGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgcgttc    <  1:2767954/150‑7 (MQ=255)
                                                                                                                                                  aTGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCg   <  2:1385888/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTCGTCCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGA  >  NC_002947/5378004‑5378300

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