Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,065,511 A→C 6.1% S151R (AGC→CGC)  PP_2677 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,065,5110AC 6.1% 203.8 / 4.6 66S151R (AGC→CGC) PP_2677hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base C (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AGCGCGTGTTGCCCTGGTTGTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCT  >  NC_002947/3065371‑3065659
                                                                                                                                            |                                                                                                                                                    
aGCGCGTGTTGCCCTGGTTGTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCgtggtg                                                                                                                                              <  7:132223/149‑1 (MQ=255)
       gtTGCCCTGGTTGTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCAcc                                                                                                                                        >  3:314698/1‑148 (MQ=255)
          gctcTGGTTGTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTCTACTGACGCGCGATCGCCAGAGGCACGTCGGATCGACCCTGATCCATGCGGCCGCCGGCGGTAGTACCGCACCCAGCGTGGTGCGCACCCAGc                                                                                                                                    >  2:167456/4‑149 (MQ=255)
          gCCCTGGTTGTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAg                                                                                                                                     <  8:245809/148‑1 (MQ=255)
                   gTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTgg                                                                                                                            <  1:167456/148‑1 (MQ=255)
                   gTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTgg                                                                                                                            <  4:28637/148‑1 (MQ=255)
                    tCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGgga                                                                                                                          <  5:187046/149‑1 (MQ=255)
                      gATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCTCCCAGccc                                                                                                                                  <  3:156222/139‑1 (MQ=255)
                      gATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCTCCCAGccc                                                                                                                                  >  4:156222/1‑139 (MQ=255)
                      gATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGggag                                                                                                                         >  4:54943/1‑148 (MQ=255)
                       atccgcatccgcATCGAACAGGCCACACGGCTGATCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGtgagga                                                                                                                       <  6:167915/149‑1 (MQ=255)
                       atccgcatccgcATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGggag                                                                                                                         >  1:275281/1‑147 (MQ=255)
                        tccgcatccgcaTCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGgaggag                                                                                                                      >  8:164743/1‑149 (MQ=255)
                          cgcatccgcatcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGc                                                                                                                     <  4:314698/148‑1 (MQ=255)
                            catccgcatcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCGGGGAGGAGCAcc                                                                                                                  <  1:298018/149‑1 (MQ=255)
                               ccgcatcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCg                                                                                                               >  8:186835/1‑149 (MQ=255)
                                 gcatcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTg                                                                                                                             >  2:8650/1‑133 (MQ=255)
                                 gcatcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTg                                                                                                                             <  1:8650/133‑1 (MQ=255)
                                 gcatcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCCGGGAGGAGCACCTCggg                                                                                                             >  7:8387/1‑149 (MQ=255)
                                  catcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCa                                                                                                                                          <  2:266397/119‑1 (MQ=255)
                                  catcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCa                                                                                                                                          >  1:266397/1‑119 (MQ=255)
                                    tcGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCCGGGAGGAGCACCTCGGGGAg                                                                                                          <  8:8387/149‑1 (MQ=255)
                                         cAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCt                                                                                                     <  7:186835/149‑1 (MQ=255)
                                          aGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTg                                                                                                    <  2:275281/149‑1 (MQ=255)
                                          aGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTg                                                                                                    <  7:164743/149‑1 (MQ=255)
                                            gCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCg                                                                                                                                 >  2:31979/1‑118 (MQ=255)
                                            gCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCg                                                                                                                                 <  1:31979/118‑1 (MQ=255)
                                            gCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTggg                                                                                                  >  3:149616/1‑149 (MQ=255)
                                              cacaCCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCACGCTGGGTgcggt                                                                                                                                          >  1:307097/1‑105 (MQ=255)
                                              cacaCCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCACGCTGGGTgcggt                                                                                                                                          <  2:307097/107‑3 (MQ=255)
                                              cacaCCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCACGCTGGGTgcggt                                                                                                                                          >  1:302079/1‑105 (MQ=255)
                                              cacaCCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCACGCTGGGTgcggt                                                                                                                                          <  2:302079/107‑3 (MQ=255)
                                               acacCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTg                                                                                                       >  3:259286/1‑141 (MQ=255)
                                               acacCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTg                                                                                                       <  4:259286/141‑1 (MQ=255)
                                                            ccgcTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGt                                                                                  <  3:54943/149‑1 (MQ=255)
                                                            ccgcTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGt                                                                                  <  5:140349/149‑1 (MQ=255)
                                                                     tGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGAcag                                                                          >  1:311380/1‑148 (MQ=255)
                                                                      gACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGAcag                                                                          >  2:127248/1‑147 (MQ=255)
                                                                         gcgcgATGGCCTGTGGCACGTCGGTTCGACCCTGCTCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGt                                                                                  >  2:92890/1‑136 (MQ=255)
                                                                         gcgcgATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGt                                                                                  <  1:92890/136‑1 (MQ=255)
                                                                             gATGGCCTGTGGCACGGCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCcgcct                                                                 <  7:324970/149‑1 (MQ=255)
                                                                                  ccTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCcgcctgcgcc                                                            >  5:48422/1‑149 (MQ=255)
                                                                                   cTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCcgcctgcgcct                                                           >  6:134068/1‑149 (MQ=255)
                                                                                   cTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCcgcctgcgcct                                                           >  3:17881/1‑149 (MQ=255)
                                                                                        gCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCGGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAg                                                       <  5:134068/148‑1 (MQ=255)
                                                                                         cACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGt                                                     <  1:127248/149‑1 (MQ=255)
                                                                                          aCGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGAcag                                                                          <  6:161467/127‑1 (MQ=255)
                                                                                          aCGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGAcag                                                                          >  5:161467/1‑127 (MQ=255)
                                                                                          aCGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTg                                                    >  5:253141/1‑149 (MQ=255)
                                                                                           cGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTgg                                                   <  2:311380/149‑1 (MQ=255)
                                                                                               ggTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCCGGTGGCCCa                                               >  2:344408/1‑149 (MQ=255)
                                                                                                 ttCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCAcc                                             <  4:17881/149‑1 (MQ=255)
                                                                                                       ccTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATg                                        >  7:201632/1‑148 (MQ=255)
                                                                                                       ccTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATg                                        >  2:307345/1‑148 (MQ=255)
                                                                                                        cTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGa                                      >  7:248184/1‑149 (MQ=255)
                                                                                                         tGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAg                                                       <  7:329064/131‑1 (MQ=255)
                                                                                                         tGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAg                                                       >  8:329064/1‑131 (MQ=255)
                                                                                                          gATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGAGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGAt                                     >  5:250916/1‑148 (MQ=255)
                                                                                                              gATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACgg                                >  6:297147/1‑149 (MQ=255)
                                                                                                                              gTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTgg                                                                                >  3:142292/1‑85 (MQ=255)
                                                                                                                               ttGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTc               <  1:307345/149‑1 (MQ=255)
                                                                                                                                  tACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCt             <  8:248184/148‑1 (MQ=255)
                                                                                                                                    ccGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGa           <  6:48422/148‑1 (MQ=255)
                                                                                                                                    ccGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGAGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGa           <  6:250916/148‑1 (MQ=255)
                                                                                                                                       cACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGttc       <  8:201632/149‑1 (MQ=255)
                                                                                                                                          ccAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGttcttt    <  1:136946/149‑1 (MQ=255)
                                                                                                                                            aGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGttctttct  >  4:112781/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
AGCGCGTGTTGCCCTGGTTGTCGATCCGCATCCGCATCGAACAGGCCACACCGCTGCGCGCCGCTGTACTGACGCGCGATGGCCTGTGGCACGTCGGTTCGACCCTGATCGATGCGGCCGGCGGCGGTTGTACCGCACCCAGCGTGGTGCGCACCCAGCCCGGCTGGGAGGAGCACCTCGGGGAGGTGCTGGGCGGGCGCTACCCGCGTGGCGACAGCAGCCGCCTGCGCCTGCAGGTGGCCCACCCCATGGATAACGGCCTGGTCAGCGGCATTCCTGAGTTCTTTCT  >  NC_002947/3065371‑3065659

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