Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 97,592 G→T 10.5% intergenic (+21/+15) PP_0090 → / ← PP_0091 hypothetical protein/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294797,5920GT10.5% 170.6 / 5.6 68intergenic (+21/+15)PP_0090/PP_0091hypothetical protein/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base T (3/5);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCAAGGGGTTTCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAGTTGTCAA  >  NC_002947/97446‑97740
                                                                                                                                                  |                                                                                                                                                     
gCCCAAGGGGTTTCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGaa                                                                                                                                                     >  4:263142/1‑149 (MQ=255)
    aaGGGGTTTCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAagcag                                                                                                                                                 >  8:102918/1‑149 (MQ=255)
      ggggTTTCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGATCATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGcc                                                                                                                                               <  8:323540/149‑1 (MQ=255)
         gTTTCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCca                                                                                                                                            >  8:166430/1‑149 (MQ=255)
          tttCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCTTCGCAGccccgg                                                                                                                                           <  7:288089/149‑3 (MQ=255)
            tCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCTTCGCAGccccggg                                                                                                                                          <  5:50308/148‑4 (MQ=255)
            tCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCcaca                                                                                                                                          <  7:229790/148‑1 (MQ=255)
               tCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTcc                                                                                                                                      >  5:164896/1‑148 (MQ=255)
               tCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTc                                                                                                                                       >  5:151364/1‑148 (MQ=255)
               tCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTc                                                                                                                                       <  6:151364/148‑1 (MQ=255)
                  tACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACgg                                                                                                                                   >  5:339439/1‑149 (MQ=255)
                      cGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATg                                                                                                                               <  3:263142/149‑1 (MQ=255)
                        cTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCa                                                                                                                             <  7:166430/149‑1 (MQ=255)
                           gagaTCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGaa                                                                                                                          <  7:102918/149‑1 (MQ=255)
                           gagaTCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCCGCGAAGCAGCCCCACACTCACGGGATGCAGaa                                                                                                                          <  6:319302/149‑1 (MQ=255)
                              aTCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCCCACTCACGGGATGCAGAacac                                                                                                                       >  3:208852/1‑149 (MQ=255)
                                       aaCAGCGAGGTCAAGGACCAGAACATGGTCGTAGGGACCTGCGATGGCGGGACGAAGAAGATCATATATCAGAAGAACCATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAgg                                                                                                             >  1:36298/1‑149 (MQ=255)
                                            cGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGt                                                                                                        <  6:164896/149‑1 (MQ=255)
                                              aGGTCAAGCACCAGAACATGGTCGGCGGGACCTGCGATGGCGGGACGAATAATATCATCTATCACAAGAACGATCTATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTc                                                                                                      <  6:137099/149‑1 (MQ=255)
                                                 tCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTcac                                                                                                    >  2:43328/1‑148 (MQ=255)
                                                  cAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTcacca                                                                                                  >  2:202563/1‑149 (MQ=255)
                                                   aaGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCGGCGAAGCCGGCCCCCCCTCACCGGCTGCAGAACCAGTCCCTAGGCTCGGTcacccc                                                                                                 >  6:26717/1‑147 (MQ=255)
                                                    aGGACCAGAACATGGTCGTCGGGACCTGCCATGGCGGGACGAAGAAGATCATCTATCAGAAGAACCATCGATAAGCATCAGGGCCCGGGGCTGGGAAGCAGCACCACACTCACGGGATGCATAACAATTCGCTACGCACGTTcaacac                                                                                                 >  7:79341/1‑148 (MQ=255)
                                                     ggACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGTGCCCGGGGCTGCGAAGCAGCCCCACACTCATGGGATGCAGAACACTTCGCTAGGCTCGTTcaccacca                                                                                               <  1:62257/149‑1 (MQ=255)
                                                          aGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGGGAAGCCGCCCCCCACTGCCCGGGTGCAGAACACGTGCCGATGCTCGTTGACCACCCCCGTg                                                                                          >  8:174768/1‑149 (MQ=255)
                                                           gAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCAAACTCACGGGATGCAGAACACATCGCTAGGCTCGTTCACCACCACCTTGc                                                                                         >  5:323843/1‑149 (MQ=255)
                                                               ttggtggTCGGGACCTGCGATGGCGGTACTAATAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGtt                                                                                     <  6:332896/148‑1 (MQ=255)
                                                                  ggcggcgGGACGTGCGATGGCGTCACTAAGAAGAACATCTATCAGAAGAACGATCGATAAGCATCCGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGc                                                                                  <  1:257150/144‑1 (MQ=255)
                                                                   tcgtcgTGTCCTCCGATGGCGCGACGAAGTAGATCATATATTAGATCAAGGATGGATTAGCTGCAGGGGCCGGGGGTGCGAAGCTGGCCCACACTCACGGGATGCAGAACACTTCGCTAGTGTCGTTCACCACCACCTTGCGGTTGGCg                                                                                 <  4:221849/149‑1 (MQ=255)
                                                                   tcgtcgGGACCTGCGTTGGCGGGACGAAGAAGATCATCTATCAGAAGAACTATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCAGCTAGGCTCGTTCACCACCACCTTGCGGTTGGc                                                                                  <  1:220329/149‑1 (MQ=255)
                                                                     gtcgGGACCTGCGCTGGGGGTACTATGACCACCATATAACATAACACCTATCTAAACGCATCAGGGCCCGGGGCTCCGACCCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTc                                                                               <  1:136157/149‑1 (MQ=255)
                                                                      tcgGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTg                                                                                    <  7:100811/143‑1 (MQ=255)
                                                                      tcgGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTg                                                                                    >  8:100811/1‑143 (MQ=255)
                                                                       cgGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCGGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCa                                                                              >  7:203969/1‑148 (MQ=255)
                                                                        gggACCTGCGATGGCGGTACGAACAAGATCATCTATCAGAACAACCATCGATATGCATCAGGGCCCGGGGCTGCTAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCata                                                                            <  7:174768/149‑1 (MQ=255)
                                                                          gACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCATGCGGTTGGCGGCata                                                                            >  7:322001/1‑147 (MQ=255)
                                                                            ccTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATagc                                                                        <  1:43328/149‑1 (MQ=255)
                                                                                             aagaagATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCat                                                                             >  7:21646/1‑127 (MQ=255)
                                                                                             aagaagATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACAACCACCATGCGGTTGGCGTCATATATCAGCACCTGCGGCTCCa                                                        >  3:37457/1‑148 (MQ=255)
                                                                                               gaagaTCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGGTGGGAAGCAGCCCCCCACACACGGGCTGGACAACACATGGCGAGGGTCGTTAGACACCCCCATGCGGGTGGGGGCATTTAGAAAAACCTGCGGCTcacac                                                      >  2:142838/1‑144 (MQ=255)
                                                                                                  catcatcTATCACACTAACGAGCGATAAGCCTCCGGCCCGGGGGGTGCGACCCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCg                                                  <  4:307521/148‑1 (MQ=255)
                                                                                                    tcattTATCATAAGAACTATCGATAAGCATCTGGGCCCGGGGGTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGgc                                                <  8:322001/149‑1 (MQ=255)
                                                                                                     catcTATCAGAAGAACTATCCATAAGCATCAGGGCCCGGGGCTGCTAAGCAGCCCCACACTCACGGGAt                                                                                                                                <  2:6980/69‑1 (MQ=255)
                                                                                                     catcTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGAt                                                                                                                                >  1:6980/1‑69 (MQ=255)
                                                                                                      atcTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGgcgc                                              >  6:309818/1‑149 (MQ=255)
                                                                                                         tATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCcg                                           >  8:53921/1‑149 (MQ=255)
                                                                                                         tATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCcg                                           >  4:121994/1‑149 (MQ=255)
                                                                                                          aTCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAACCCCCCACTCAAGGGGTTGAAAACACATACCGAGGCGGGTTAACCACCCCCATTCGGTTGGCGTCCTATAACAGCCCCAGCGGCTCCAAAACCCGCGccccc                                          >  4:318567/1‑147 (MQ=255)
                                                                                                            cAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCgcgc                                        >  3:289939/1‑149 (MQ=255)
                                                                                                                  aaCGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAg                                  <  6:339439/149‑1 (MQ=255)
                                                                                                                   aCGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCCGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAgg                                 <  8:203969/149‑1 (MQ=255)
                                                                                                                        cgatAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCAc                                                                   <  6:87307/110‑1 (MQ=255)
                                                                                                                        cgatAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCAc                                                                   >  5:87307/1‑110 (MQ=255)
                                                                                                                        cgatAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCAc                                                                   >  5:82553/1‑110 (MQ=255)
                                                                                                                        cgatAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCAc                                                                   <  6:82553/110‑1 (MQ=255)
                                                                                                                          atAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACACACGGGATGCAGAACACATCGCTAGGCGCGGTAACCACCACCATGCGGGTGGCGGGATATAGCAGCACCAGCGGCGACCTCACCGGCGCCCGCGCGTCCCGGCGATAgc                          >  4:49168/1‑149 (MQ=255)
                                                                                                                          atAAGCATCAGGGCCCGGGGCGGGGAAGCAGCCCCACACACACCGGGTGCAGAACACATAGCGAGGCTCCTTCACCCCCCCCGTGCGGGTGGGGGGCAATAACAGCACCCGGGGCGCCAAAACCCGCGCCCGCGCGTCCAGGCGCTAgc                          >  6:217319/1‑149 (MQ=255)
                                                                                                                                 tCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGcc                   >  6:148458/1‑149 (MQ=255)
                                                                                                                                   aGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGACg                  <  7:53921/148‑1 (MQ=255)
                                                                                                                                       tgtggggCTGCTTCGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCtt             >  4:37296/4‑149 (MQ=255)
                                                                                                                                       tggggggCTGCTTCGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAgcg                         <  7:119693/134‑1 (MQ=255)
                                                                                                                                       cccGGGGCTGCGAAGCAGCCCCACACTCCCGGGATGCAGAACACTTCGCTCGGCGCGTTCACCACCGCCTTGCGGGTGGCGGCAGTTAGCAGCGCCTGCGGCTCCGTCCCCGGCGCCCGCGCGTCCCGGCGAGAGCGCGCGCCGGcggg             >  2:248193/1‑146 (MQ=255)
                                                                                                                                        gtggggCTGCTTCGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTg            >  8:118952/3‑149 (MQ=255)
                                                                                                                                           gggCTGCTTCGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAg         >  2:120317/1‑149 (MQ=255)
                                                                                                                                           gggCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAg         <  2:36298/149‑1 (MQ=255)
                                                                                                                                            ggCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAGt        <  5:309818/149‑1 (MQ=255)
                                                                                                                                             gCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAGtt       >  6:350111/1‑149 (MQ=255)
                                                                                                                                                 cGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAGTTGTc    >  4:308551/1‑148 (MQ=255)
                                                                                                                                                  gAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAGTTGTCaa  <  4:289939/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCCAAGGGGTTTCGTCCTACACGCTGGAGATCGTGCCCAACAGCGAGGTCAAGGACCAGAACATGGTCGTCGGGACCTGCGATGGCGGGACGAAGAAGATCATCTATCAGAAGAACGATCGATAAGCATCAGGGCCCGGGGCTGCGAAGCAGCCCCACACTCACGGGATGCAGAACACTTCGCTAGGCTCGTTCACCACCACCTTGCGGTTGGCGTCATATAGCAGCACCTGCGGCTCCATCACCGGCGCCCGCGCTTCCAGGCGATAGCGCTCGCCGGCCTTGAAGTTGTCAA  >  NC_002947/97446‑97740

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