Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,378,320 G→A 5.1% G405G (GGC→GGT dnaK ← chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,378,3200GA 5.1% 496.1 / 10.3 156G405G (GGC→GGTdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (64/84);  new base A (4/4);  total (68/88)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 161 total aligned reads displayed.

CCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAA  >  NC_002947/5378172‑5378466
                                                                                                                                                    |                                                                                                                                                  
ccAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCAcgccgc                                                                                                                                                   >  1:1318060/1‑150 (MQ=255)
     gACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATg                                                                                                                                              <  2:5011943/150‑1 (MQ=255)
      aCTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATgg                                                                                                                                             >  1:1491682/1‑150 (MQ=255)
         tGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGttt                                                                                                                                          >  1:266141/1‑150 (MQ=255)
         tGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGttt                                                                                                                                          <  1:575045/150‑1 (MQ=255)
          gTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTc                                                                                                                                         <  1:2967868/150‑1 (MQ=255)
          gTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTc                                                                                                                                         <  1:3274744/150‑1 (MQ=255)
          gTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTc                                                                                                                                         <  1:3274801/150‑1 (MQ=255)
          gTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTc                                                                                                                                         <  1:2087792/150‑1 (MQ=255)
           ttCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCg                                                                                                                                        >  2:4131879/1‑150 (MQ=255)
              tGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATa                                                                                                                                     >  1:3456544/1‑150 (MQ=255)
               gcgcAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATAc                                                                                                                                    >  1:2156371/1‑150 (MQ=255)
                 gcAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATAccc                                                                                                                                  <  2:591486/150‑1 (MQ=255)
                    gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                                                                                               <  2:1059392/150‑1 (MQ=255)
                    gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                                                                                               >  1:4979264/1‑150 (MQ=255)
                    gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                                                                                               >  1:4597939/1‑150 (MQ=255)
                    gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                                                                                               <  2:3175843/150‑1 (MQ=255)
                    gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                                                                                               <  2:324049/150‑1 (MQ=255)
                    gCCTGCTTACGCTCACCTTGCAGCACGTGAACGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                                                                                                               <  1:4430337/150‑1 (MQ=255)
                         cTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAg                                                                                                                          >  1:2558484/1‑150 (MQ=255)
                          ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                                                                                                                         >  2:2434770/1‑150 (MQ=255)
                          ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                                                                                                                         >  2:3615569/1‑150 (MQ=255)
                          ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                                                                                                                         <  1:729601/150‑1 (MQ=255)
                            aCGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCgg                                                                                                                       <  2:3983868/150‑1 (MQ=255)
                            aCGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCgg                                                                                                                       >  2:4409924/1‑150 (MQ=255)
                             cGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCggg                                                                                                                      <  2:2705490/150‑1 (MQ=255)
                             cGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCggg                                                                                                                      >  1:2138916/1‑150 (MQ=255)
                              gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                     >  2:1864030/1‑150 (MQ=255)
                              gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                     >  2:4348437/1‑150 (MQ=255)
                              gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                     >  2:4348531/1‑150 (MQ=255)
                              gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                     >  1:1525341/1‑150 (MQ=255)
                              gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                                                                                                     >  2:59882/1‑150 (MQ=255)
                                 cACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGa                                                                                                                  <  1:1117191/150‑1 (MQ=255)
                                  aCCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGAc                                                                                                                 <  2:799997/150‑1 (MQ=255)
                                    catGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGACAACACCTGCGACTTCTTGGTCGGGATGCAGCTCTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCTATAACCAGTGTCAGCGGGCTGACGt                                                                                                               >  2:3026450/3‑150 (MQ=255)
                                     ttGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                                                                                              <  2:2556631/150‑1 (MQ=255)
                                     gtGCAGCACGTGACTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTc                                                                                                              <  2:2844680/149‑1 (MQ=255)
                                       gcagcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCca                                                                                                            <  2:2403828/150‑1 (MQ=255)
                                              gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                                                                                                     <  1:4131879/150‑1 (MQ=255)
                                              gTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                                                                                                     >  1:4121886/1‑150 (MQ=255)
                                                 aaTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGAc                                                                                                  >  1:4514407/1‑150 (MQ=255)
                                                 aaTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGAc                                                                                                  <  2:4298642/150‑1 (MQ=255)
                                                   tGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGTGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGt                                                                                                <  1:3821940/150‑1 (MQ=255)
                                                       cACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACTCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCttt                                                                                            <  2:3294997/150‑1 (MQ=255)
                                                       cACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCttt                                                                                            <  2:2311095/150‑1 (MQ=255)
                                                        aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc                                                                                           >  1:138938/1‑150 (MQ=255)
                                                        aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc                                                                                           >  1:597905/1‑150 (MQ=255)
                                                           gcgcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACg                                                                                        >  2:2027438/1‑150 (MQ=255)
                                                           gcgcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACg                                                                                        <  2:3980200/150‑1 (MQ=255)
                                                             gcTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTc                                                                                      >  2:3186892/1‑150 (MQ=255)
                                                             gcTCGGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCGTTCACGTc                                                                                      >  1:106887/1‑150 (MQ=255)
                                                               tctGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCAc                                                                                    >  1:1634241/1‑150 (MQ=255)
                                                                    ttGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGTCCTTGACGTCGccctgc                                                                               >  2:2677996/1‑146 (MQ=255)
                                                                    ttGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGcc                                                                               <  1:2791644/150‑1 (MQ=255)
                                                                     tGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCa                                                                              <  2:3498937/150‑1 (MQ=255)
                                                                     tGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCa                                                                              <  2:347614/150‑1 (MQ=255)
                                                                       tcgtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTTTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCATGTCACCGGCCAGt                                                                            >  2:2381026/1‑150 (MQ=255)
                                                                         gtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATAc                                                                                                                                    >  1:3000257/1‑92 (MQ=255)
                                                                         gtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATAc                                                                                                                                    <  2:3000257/92‑1 (MQ=255)
                                                                          tcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAcc                                                                         <  1:2199144/150‑1 (MQ=255)
                                                                          tcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAcc                                                                         <  1:2869276/150‑1 (MQ=255)
                                                                               aGTCGAGAACACCTGCGACTTCTTAGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCAcc                                                                    <  1:4666834/150‑1 (MQ=255)
                                                                                   gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTgg                                                                >  1:3770570/1‑150 (MQ=255)
                                                                                   gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTgg                                                                >  2:1118944/1‑150 (MQ=255)
                                                                                     gaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGAt                                                              <  2:4571466/150‑1 (MQ=255)
                                                                                       acacCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATgg                                                            <  2:469152/150‑1 (MQ=255)
                                                                                         acCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                          <  1:1221719/150‑1 (MQ=255)
                                                                                          ccTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcag                                                         <  2:1566880/150‑1 (MQ=255)
                                                                                           cTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcagc                                                        <  1:2929225/150‑1 (MQ=255)
                                                                                           cTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcagc                                                        >  1:1844019/1‑150 (MQ=255)
                                                                                            tGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCTCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcagca                                                       >  2:1154128/1‑150 (MQ=255)
                                                                                            tGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgcagca                                                       >  2:3417090/1‑150 (MQ=255)
                                                                                                 cttcttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCAt                                                  <  1:5019528/150‑1 (MQ=255)
                                                                                                 cttcttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCAt                                                  <  1:4032677/150‑1 (MQ=255)
                                                                                                  ttcttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATg                                                 <  2:4424048/150‑1 (MQ=255)
                                                                                                       ggTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAAc                                            >  2:2720840/1‑150 (MQ=255)
                                                                                                         tCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACgg                                          >  2:1409580/1‑150 (MQ=255)
                                                                                                            ggATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCtt                                       <  1:232066/150‑1 (MQ=255)
                                                                                                            ggATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCtt                                       >  1:214415/1‑150 (MQ=255)
                                                                                                                ggtggtGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGGTTCGAGACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTc                                   <  2:1491682/150‑1 (MQ=255)
                                                                                                                      gtTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCGACGGCTTCGTCCGGGtt                             >  1:626489/1‑150 (MQ=255)
                                                                                                                      gtTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGtt                             >  1:4967446/1‑150 (MQ=255)
                                                                                                                         cttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGAc                          >  1:530248/1‑150 (MQ=255)
                                                                                                                          ttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                          <  2:1945466/117‑1 (MQ=255)
                                                                                                                          ttctCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACg                         <  2:4113857/150‑1 (MQ=255)
                                                                                                                            ctcGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTc                       <  2:1247760/150‑1 (MQ=255)
                                                                                                                               gATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCtt                    <  1:2582633/150‑1 (MQ=255)
                                                                                                                               gATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCtt                    >  2:5098158/1‑150 (MQ=255)
                                                                                                                                 tCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgc                  >  2:1044428/1‑150 (MQ=255)
                                                                                                                                  cAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgcg                 <  1:607702/150‑1 (MQ=255)
                                                                                                                                         gTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGTCTTCGTCCGGGTTGACGTCCTTGCGCGcttct          >  1:3931811/1‑150 (MQ=255)
                                                                                                                                         gTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGcttct          >  1:3931880/1‑150 (MQ=255)
                                                                                                                                            atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                          <  2:1602590/99‑1 (MQ=255)
                                                                                                                                            atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTa       >  1:1615552/1‑150 (MQ=255)
                                                                                                                                            atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTa       <  1:4543832/150‑1 (MQ=255)
                                                                                                                                            atcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTa       >  1:4413941/1‑150 (MQ=255)
                                                                                                                                             tcaCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAc      >  1:919305/1‑150 (MQ=255)
                                                                                                                                              caCGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAcc     >  1:2033246/1‑150 (MQ=255)
                                                                                                                                                 gccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa  <  2:4597939/150‑1 (MQ=255)
                                                                                                                                                 gccgccCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaa  <  2:4979264/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAA  >  NC_002947/5378172‑5378466

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: