Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,373,823 G→A 6.4% S44F (TCT→TTT) ‡ PP_1199 ← toxin‑antitoxin system HigB‑like toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,373,8230GA 6.4% 236.8 / 3.6 78S44F (TCT→TTT) ‡PP_1199toxin‑antitoxin system HigB‑like toxin
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base A (3/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

ACGATGCCCTTGCGATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGCTGGATAATA  >  NC_002947/1373679‑1373971
                                                                                                                                                |                                                                                                                                                    
aCGATGCCCTTGCGATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTCTCTCCTGGAGGAAAcc                                                                                                                                                  >  2:44164/1‑149 (MQ=255)
    tGCCCTTGCGATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTCTCTCCTGGAGGAAACCGGcg                                                                                                                                              <  1:44164/149‑2 (MQ=255)
     gCCCTTGCGATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCtt                                                                                                                                             >  6:61373/1‑149 (MQ=255)
           gCGATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGa                                                                                                                                        <  6:277253/148‑1 (MQ=255)
             gATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTCTCTCCTGGAGGAAACCGGc                                                                                                                                               >  8:224587/1‑139 (MQ=255)
             gATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTCTCTCCTGGAGGAAACCGGc                                                                                                                                               <  7:224587/139‑1 (MQ=255)
             gATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCCCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATg                                                                                                                                      >  2:206633/1‑148 (MQ=255)
                  aTCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCGTCGGGATGCTCGGc                                                                                                                                >  4:225225/1‑149 (MQ=255)
                       aaTCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGa                                                                                                                            >  2:97549/1‑148 (MQ=255)
                         tCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGTCAGATAg                                                                                                                         <  5:239116/149‑1 (MQ=255)
                           aTACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAgg                                                                                                                        <  1:165230/148‑1 (MQ=255)
                            tACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAAcc                                                                                                                                                  >  7:318850/1‑121 (MQ=255)
                            tACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAAcc                                                                                                                                                  <  8:318850/121‑1 (MQ=255)
                             aCGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCt                                                                                                                     >  2:230479/1‑149 (MQ=255)
                                       tggGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCg                                                                                                           <  4:89735/149‑1 (MQ=255)
                                       tggGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCg                                                                                                           >  6:353430/1‑149 (MQ=255)
                                       tggGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCg                                                                                                           <  5:61373/149‑1 (MQ=255)
                                            tGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTAtt                                                                                                      >  3:274056/1‑149 (MQ=255)
                                              cAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTAttct                                                                                                    <  4:274056/149‑1 (MQ=255)
                                                        gtagtaGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCa                                                                                          <  3:205304/149‑1 (MQ=255)
                                                           gtagATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGt                                                                                       <  7:122109/149‑1 (MQ=255)
                                                                 cACTCGTACACCCCCACTCTTGCCCATGCCCTAGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGCTTTCCTCAAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGTTAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAggagatcat                                                                                 >  1:331398/1‑142 (MQ=255)
                                                                   cTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCCGTATTCTTCATCATCCAGGTGCTCCTTg                                                                               <  1:230479/149‑1 (MQ=255)
                                                                     cGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCtt                                                                                >  3:25577/1‑146 (MQ=255)
                                                                          acCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTcc                                                                                  >  2:280892/1‑139 (MQ=255)
                                                                          acCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTcc                                                                                  <  1:280892/139‑1 (MQ=255)
                                                                          acCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCt                                                                                 <  7:179049/140‑1 (MQ=255)
                                                                          acCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCt                                                                                 >  8:179049/1‑140 (MQ=255)
                                                                           cccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCTGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGa                                                                       <  5:362700/149‑1 (MQ=255)
                                                                           cccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTtcatca                                                                                              <  3:50093/126‑1 (MQ=255)
                                                                           cccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTtcatca                                                                                              >  4:50093/1‑126 (MQ=255)
                                                                            ccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCtt                                                                                >  2:50448/1‑139 (MQ=255)
                                                                            ccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCtt                                                                                >  4:200501/1‑139 (MQ=255)
                                                                            ccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCggtg                                                                      >  5:364874/1‑147 (MQ=255)
                                                                            ccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCg                                                                         >  3:37865/1‑146 (MQ=255)
                                                                            ccccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCg                                                                         <  4:37865/146‑1 (MQ=255)
                                                                              ccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGt                                                                    >  4:21534/1‑149 (MQ=255)
                                                                              ccACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGt                                                                    >  3:133620/1‑149 (MQ=255)
                                                                               cACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTg                                                                   <  4:25577/149‑1 (MQ=255)
                                                                                aCTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTg                                                                   >  5:213590/1‑148 (MQ=255)
                                                                                  tctTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGt                                                                    <  4:173744/145‑1 (MQ=255)
                                                                                  tctTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGt                                                                    >  3:173744/1‑145 (MQ=255)
                                                                                  tctTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGaa                                                                 >  8:39035/1‑148 (MQ=255)
                                                                                    ttgcccttgcccTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGTGGCGAATAc                                                              >  8:372342/1‑149 (MQ=255)
                                                                                      gcccttgcccttGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTAttcttc                                                                                                  >  1:384189/1‑111 (MQ=255)
                                                                                      gcccttgcccttGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTAttcttc                                                                                                  <  2:384189/111‑1 (MQ=255)
                                                                                       cccttgcccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTg                                                            <  2:247609/148‑1 (MQ=255)
                                                                                        ccttgccctttGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTggg                                                          <  1:30211/149‑1 (MQ=255)
                                                                                        ccttgcccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTggg                                                          <  1:97549/149‑1 (MQ=255)
                                                                                        ccttgcccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTggg                                                          >  2:148088/1‑149 (MQ=255)
                                                                                        ccttgcccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTggg                                                          <  3:200501/149‑1 (MQ=255)
                                                                                         cttgcccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTgggg                                                         >  7:314387/1‑149 (MQ=255)
                                                                                         cttgcccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTgggg                                                         >  7:213103/1‑149 (MQ=255)
                                                                                              ccttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGtctc                                                     >  1:411418/1‑148 (MQ=255)
                                                                                               cttgGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGa                                                   <  8:213103/149‑1 (MQ=255)
                                                                                                  ggCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGAt                                                  <  5:40819/147‑1 (MQ=255)
                                                                                                   gCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGaa                                               <  3:21534/149‑1 (MQ=255)
                                                                                                     tGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGa                                             <  1:206633/149‑1 (MQ=255)
                                                                                                       cccAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAg                                                                      >  5:82288/1‑122 (MQ=255)
                                                                                                       cccAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAg                                                                      <  6:82288/122‑1 (MQ=255)
                                                                                                        ccAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTga                                                                              >  3:239785/1‑113 (MQ=255)
                                                                                                        ccAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTga                                                                              <  4:239785/113‑1 (MQ=255)
                                                                                                        ccAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCa                                          <  6:213590/149‑1 (MQ=255)
                                                                                                          aGCTGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGTTGCTCGGCCAGTTAGGCCTGCAAAGCCCGGTATTCTTCATCTTCGAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTGTCGTTGAAGATCATg                                        <  2:331398/149‑1 (MQ=255)
                                                                                                              gATCTTTCTGAGGCCACCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGGCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGa                                                <  6:415938/137‑1 (MQ=255)
                                                                                                              gATCTTTCTGAGGCCACCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGa                                                >  5:415938/1‑137 (MQ=255)
                                                                                                                 cTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTCCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGCAGGTGAATACTGGGGTCTCGAAGAAGATCATGCAGTca                                  >  4:44091/1‑147 (MQ=255)
                                                                                                                   ttCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTg                               >  5:242767/1‑149 (MQ=255)
                                                                                                                     cTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGc                              <  7:39035/148‑1 (MQ=255)
                                                                                                                          gccgccGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCg                        >  7:239423/1‑149 (MQ=255)
                                                                                                                           ccgccgGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGt                       <  1:148088/149‑1 (MQ=255)
                                                                                                                                  tttCCTCCTGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCTTCCATGTGCTCCTTGAGATCGGAGGTGAATACTGGGTTCTCGATGAAGTTCATGCAGTCCCTGCTATGCGTCAATGGc                <  3:44091/149‑1 (MQ=255)
                                                                                                                                    ctctccTGGAGGAAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCg               >  1:370404/3‑148 (MQ=255)
                                                                                                                                     cctccAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGt              <  8:239423/148‑1 (MQ=255)
                                                                                                                                      ctccAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAg            <  6:364874/149‑1 (MQ=255)
                                                                                                                                       tccAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGc           >  2:185692/1‑149 (MQ=255)
                                                                                                                                       tccAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGc           >  6:170238/1‑149 (MQ=255)
                                                                                                                                         cAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGCTg         >  2:206009/1‑149 (MQ=255)
                                                                                                                                                gaACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGCTGGataata  <  1:50448/149‑1 (MQ=255)
                                                                                                                                                gaACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGCTGGataat   <  8:314387/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACGATGCCCTTGCGATAGATCAGAATCATACGGATCTGGTGGGCTGCAGTGACGTGGTAGTAGATCACTCGTACACCCCCACTCTTGCCCTTGCCCTTGGCTGCCCAGCGGATCTTTCTGAGGCCGCCGGTTTCCTCCAGGAGAGAACCGGCTTCGGGATGCTCGGCCAGATAGGCCTGCAAAGCCCGGTATTCTTCATCATCCAGGTGCTCCTTGAGATCGGAGGTGAATACTGGGGTCTCGATGAAGATCATGCAGTCCCTGCTATGCGTCAATGGCGTAGCTGGATAATA  >  NC_002947/1373679‑1373971

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