Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,261,548 Δ1 bp 7.4% coding (1087/1116 nt) PP_1104 ← AotO‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,261,5480A. 7.4% 211.6 / 10.8 68coding (1087/1116 nt)PP_1104AotO‑related protein
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base . (3/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

CTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTT  >  NC_002947/1261409‑1261691
                                                                                                                                            |                                                                                                                                               
cTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTc                                                                                                                                        <  2:92561/149‑1 (MQ=255)
cTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTc                                                                                                                                        >  5:185562/1‑149 (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                                                                                                              >  2:215544/1‑135 (MQ=255)
                                        aTGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGc                                                                                                              <  1:215544/135‑1 (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                                                                       <  2:39331/145‑1 (MQ=255)
                                                                     gCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATAcagca                                                                       >  1:39331/1‑145 (MQ=255)
                                                                      cTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGcccc                                                                  <  3:110493/149‑1 (MQ=255)
                                                                      cTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGccc                                                                   <  6:185562/148‑1 (MQ=255)
                                                                         cTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATa                                                                            >  8:135941/1‑136 (MQ=255)
                                                                         cTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATa                                                                            <  7:135941/136‑1 (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                                         >  1:253961/1‑149 (MQ=255)
                                                                                               aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                         >  8:107878/1‑149 (MQ=255)
                                                                                               aaTGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGa                                         >  3:325448/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                           >  4:297673/1‑141 (MQ=255)
                                                                                                                     gCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGt                           <  3:297673/141‑1 (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)
                                                                                                                                        cGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGtt  <  2:320015/148‑1 (MQ=255)
                                                                                                                                        cGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGtt  <  2:314739/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
CTTGGCAAACCAGGCATTGGCGAGGTCGGTAGCAGCAGGCATGAGGGCAGGCGGTCCTTCGAAGGCTGTGCTTCTACAGACCACTTGCCGGCTGCAATGGTTTTACCGGATCAATTGGCCAACAGGTA‑GCCGCTGCGGCATACCTGTTCGCCGGCCACATGGGCGATGACCATGCCAGCGGTGGCCATGCGCCGGACGCTACGCGCATACAGCAGCCCCGGCTGGCTGTTGCGCAGGGGGGTGACTCGGTCACTGAGTGGGTCGATGACTTCGGCGATCAGTT  >  NC_002947/1261409‑1261691

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