Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,542 T→A 7.6% S365C (AGC→TGC)  PP_1104 ← AotO‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,5420TA 7.6% 195.0 / 5.5 66S365C (AGC→TGC) PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base T (34/27);  new base A (3/2);  total (37/29)
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.85e-01

CTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCA  >  NC_002947/1261409‑1261688
                                                                                                                                      |                                                                                                                                                  
cTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTc                                                                                                                                     >  5:185562/1‑149 (MQ=255)
cTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTc                                                                                                                                     <  2:92561/149‑1 (MQ=255)
     cAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGCCCAACAGGTA‑GCCACTGCGGCATACCTGTTCGCCgg                                                                                                                                >  3:177093/1‑149 (MQ=255)
          cAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCCGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCcaca                                                                                                                           >  3:100044/1‑149 (MQ=255)
             gCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATgg                                                                                                                        >  6:14546/1‑149 (MQ=255)
               aTTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGc                                                                                                                      >  2:184956/1‑149 (MQ=255)
                                  gcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTg                                                                                                   >  4:110493/1‑149 (MQ=255)
                                  gcagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGAAGACCATGCCAGCGGTg                                                                                                   <  7:105549/149‑1 (MQ=255)
                                   tagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTgg                                                                                                  <  5:215802/148‑1 (MQ=255)
                                   cagGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTgg                                                                                                  <  1:184956/149‑1 (MQ=255)
                                     ggCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTgg                                                                                                  >  6:215802/1‑147 (MQ=255)
                                       cATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCAt                                                                                              <  6:35813/149‑1 (MQ=255)
                                        aTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                           <  1:215544/135‑1 (MQ=255)
                                        aTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                           >  2:215544/1‑135 (MQ=255)
                                        aTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCAATTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATg                                                                                             >  4:305465/1‑149 (MQ=255)
                                         tGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATgg                                                                                                                        >  8:28557/1‑121 (MQ=255)
                                         tGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATgg                                                                                                                        <  7:28557/121‑1 (MQ=255)
                                                   cGGTCCTTCGAAGGCCGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCt                                                                                  <  4:100044/149‑1 (MQ=255)
                                                     gTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCTTGCGCCGGACGCTAc                                                                                >  3:38861/1‑149 (MQ=255)
                                                     gTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAc                                                                                >  1:358084/1‑149 (MQ=255)
                                                     gTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCAGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAc                                                                                >  4:59077/1‑149 (MQ=255)
                                                        cTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTAcgc                                                                              <  6:249732/148‑1 (MQ=255)
                                                          tCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GACGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGACATGCGCCGGACGCTTCGCGCa                                                                           >  5:41896/1‑149 (MQ=255)
                                                              aGGCTGTGCTTCTACAGACCAATTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAca                                                                       <  3:305465/149‑1 (MQ=255)
                                                                 cTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagca                                                                    <  5:14546/149‑1 (MQ=255)
                                                                     gCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagca                                                                    >  1:39331/1‑145 (MQ=255)
                                                                     gCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagca                                                                    <  2:39331/145‑1 (MQ=255)
                                                                      cTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGcccc                                                               <  3:110493/149‑1 (MQ=255)
                                                                      cTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGccc                                                                <  6:185562/148‑1 (MQ=255)
                                                                         cTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATa                                                                         <  7:135941/136‑1 (MQ=255)
                                                                         cTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATa                                                                         >  8:135941/1‑136 (MQ=255)
                                                                            cAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggct                                                           <  2:358084/147‑1 (MQ=255)
                                                                             aGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctggc                                                        >  4:144572/1‑149 (MQ=255)
                                                                             aGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCggctggc                                                        >  2:73680/1‑149 (MQ=255)
                                                                              gACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATggg                                                                                                                       >  3:39393/1‑85 (MQ=255)
                                                                              gACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATggg                                                                                                                       <  4:39393/85‑1 (MQ=255)
                                                                               aCCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGgctggctg                                                      >  3:349342/1‑149 (MQ=255)
                                                                                    ttGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTgcg                                                 >  7:122086/1‑149 (MQ=255)
                                                                                          gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGTATACAGCAGCCCCGGCTGGCTGTTGCGCAgg                                             <  1:32590/147‑1 (MQ=255)
                                                                                          gCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAgg                                             >  2:32590/1‑147 (MQ=255)
                                                                                              cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGt                                        <  4:349342/148‑1 (MQ=255)
                                                                                              cAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGAAGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAgggggg                                         <  1:73680/147‑1 (MQ=255)
                                                                                               aaaggTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTg                                       >  2:92440/4‑148 (MQ=255)
                                                                                               aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                      >  8:107878/1‑149 (MQ=255)
                                                                                               aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                      >  3:325448/1‑149 (MQ=255)
                                                                                               aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                      >  1:253961/1‑149 (MQ=255)
                                                                                                aTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGTGTGAc                                     <  8:209506/149‑1 (MQ=255)
                                                                                                           ggATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACt                            >  5:174348/1‑147 (MQ=255)
                                                                                                             aTCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                        <  2:253961/149‑1 (MQ=255)
                                                                                                                aaTTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTgg                      <  3:144572/148‑1 (MQ=255)
                                                                                                                  ccggCGAACAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                          <  6:98855/124‑1 (MQ=37)
                                                                                                                  ccggCGAACAGGTATGCCGCAGCGGC‑TACCTGTTGGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAggggg                                          >  5:98855/3‑126 (MQ=37)
                                                                                                                     gCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                        <  3:297673/141‑1 (MQ=255)
                                                                                                                     gCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                        >  4:297673/1‑141 (MQ=255)
                                                                                                                     gCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGAt                <  3:59077/149‑1 (MQ=255)
                                                                                                                      ccAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATg               <  8:122086/149‑1 (MQ=255)
                                                                                                                       cAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                    <  7:101289/111‑1 (MQ=255)
                                                                                                                       cAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGtt                                                    >  8:101289/1‑111 (MQ=255)
                                                                                                                       cAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGAt                >  7:183927/1‑147 (MQ=255)
                                                                                                                       cAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGa              >  5:279709/1‑149 (MQ=255)
                                                                                                                          cAGGTATGCCGCAGCGGC‑TACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACtt           >  6:40448/1‑149 (MQ=255)
                                                                                                                           aggtatgCCGCAGCGGC‑TACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTc          >  6:60832/7‑149 (MQ=255)
                                                                                                                           aGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTc          >  1:314739/1‑149 (MQ=255)
                                                                                                                           aGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTc          >  1:320015/1‑149 (MQ=255)
                                                                                                                                  ccGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGa     >  6:87161/1‑148 (MQ=255)
                                                                                                                                   cGCAGCGGC‑TACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCa  <  5:40448/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
CTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCA  >  NC_002947/1261409‑1261688

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: ‑