Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,541,781 C→A 7.3% Y166* (TAC→TAA) ‡ nhaB → Na(+)/H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,541,7810CA 7.3% 217.7 / 5.1 68Y166* (TAC→TAA) ‡nhaBNa(+)/H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base A (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 = 9.43e-01

TGTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCG  >  NC_002947/4541633‑4541923
                                                                                                                                                    |                                                                                                                                              
tgttcTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTac                                                                                                                                                >  3:368316/1‑149 (MQ=255)
       aCGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcg                                                                                                                                         >  6:23996/1‑149 (MQ=255)
         gCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtc                                                                                                                                       >  1:336782/1‑149 (MQ=255)
           aTTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgc                                                                                                                                     >  1:144143/1‑149 (MQ=255)
             tttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACcgcgtcgcgt                                                                                                                                   >  5:171078/1‑149 (MQ=255)
              ttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgc                                                                                                                                          >  1:138749/1‑141 (MQ=255)
              ttGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGGTACACcgc                                                                                                                                          <  2:138749/141‑1 (MQ=255)
                   cGGTGTGCGTTCCAAGGCAATTCTGGCGATGCTTTTTTGCTTACTTTCGGCGTTTCTATCCGCGTTTCTAGATGCCCTGACCGGCACGGCGGTTATCATCAGCGCTGGGGTGGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                              <  2:187567/148‑1 (MQ=255)
                    ggggTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                              <  1:85074/147‑1 (MQ=255)
                    ggggTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCggg                                                                                                                              <  6:221583/147‑1 (MQ=255)
                            ttCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcgc                                                                                                                    <  3:361015/149‑1 (MQ=255)
                            ttCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCcgc                                                                                                                    >  7:142526/1‑149 (MQ=255)
                                aacgCCAGTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTAACCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCAGGGGCAATCCGTGCGa                                                                                                                <  4:317575/146‑1 (MQ=255)
                                aaGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGa                                                                                                                <  6:297789/149‑1 (MQ=255)
                                     cATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGa                                                                                                             >  3:44872/1‑147 (MQ=255)
                                       ttCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAgc                                                                                                         <  4:368316/149‑1 (MQ=255)
                                              gctgctGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTg                                                                                                   >  7:199258/1‑148 (MQ=255)
                                                     ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGc                                                                                                                            >  7:14718/1‑116 (MQ=255)
                                                     ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGc                                                                                                                            <  8:14718/116‑1 (MQ=255)
                                                     ttttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGa                                                                                           <  5:23996/149‑1 (MQ=255)
                                                       ttGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGAc                                                                                          <  6:96055/148‑1 (MQ=255)
                                                          cGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAg                                                                                       >  8:34970/1‑148 (MQ=255)
                                                                 tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGc                                                                               >  3:308708/1‑149 (MQ=255)
                                                                 tCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGc                                                                               >  4:323003/1‑149 (MQ=255)
                                                                   ggCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCt                                                                              <  3:323003/148‑1 (MQ=255)
                                                                   ggCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTc                                                                             <  2:144143/149‑1 (MQ=255)
                                                                    gCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCa                                                                            >  7:261554/1‑149 (MQ=255)
                                                                     cGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCAc                                                                           >  8:340981/1‑149 (MQ=255)
                                                                       ggtctctcCGCGGGGCTCGCGGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTCCGCGGTGTACCACCGCGTCGCGTCCGGGGCCACTCCGCGCGAGGACAGCGCACTGGTGAGCGACTTGCTGGTCGCTctcnn                                                                         <  2:336782/147‑5 (MQ=255)
                                                                        ttCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGAGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGATAGCGACCAGCAAGTCGCTCACCTg                                                                        >  3:245030/1‑149 (MQ=255)
                                                                           tATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACGGCGACGCGTCCGGGGACAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAc                                                                     >  1:119578/1‑149 (MQ=255)
                                                                           tATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCAc                                                                     >  6:145952/1‑149 (MQ=255)
                                                                              ccGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                             >  3:108503/1‑90 (MQ=255)
                                                                              ccGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCgggg                                                                                                                             <  4:108503/90‑1 (MQ=255)
                                                                                gcgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCggg                                                                 >  1:13349/1‑148 (MQ=255)
                                                                                 cgTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGaa                                                               <  8:199258/149‑1 (MQ=255)
                                                                                         gATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGa                                                       <  5:145952/149‑1 (MQ=255)
                                                                                          aTGCCCTGACCGTCACGGCTGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGTTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGAc                                                      <  4:210692/149‑1 (MQ=255)
                                                                                             cccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCa                                                    <  7:34970/148‑1 (MQ=255)
                                                                                             cccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCa                                                    <  7:340981/148‑1 (MQ=255)
                                                                                              ccTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGAGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGt                                                  >  1:136052/1‑149 (MQ=255)
                                                                                                 gACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTcc                                               >  5:211192/1‑149 (MQ=255)
                                                                                                  aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCg                                                        >  6:353599/1‑139 (MQ=255)
                                                                                                  aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCg                                                        <  5:353599/139‑1 (MQ=255)
                                                                                                  aCCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTcc                                               <  8:261554/148‑1 (MQ=255)
                                                                                                     gTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                             <  2:362271/147‑1 (MQ=255)
                                                                                                     gTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCcgc                                             >  1:362271/1‑147 (MQ=255)
                                                                                                      tCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGAc                                                                                          >  2:192609/1‑101 (MQ=255)
                                                                                                              gTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCg                                  >  4:43570/1‑149 (MQ=255)
                                                                                                                 atcatcAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGACGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTa                                <  2:119578/148‑1 (MQ=255)
                                                                                                                   catcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGcc                             >  7:162301/1‑149 (MQ=255)
                                                                                                                     tcaGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCct                            <  4:44872/148‑1 (MQ=255)
                                                                                                                       aGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCctgc                          <  6:58401/148‑1 (MQ=255)
                                                                                                                           cTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGc                              >  4:115679/1‑140 (MQ=255)
                                                                                                                           cTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGc                              <  3:115679/140‑1 (MQ=255)
                                                                                                                            tGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGc                    <  6:211192/149‑1 (MQ=255)
                                                                                                                              cGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGc                    <  6:171078/147‑1 (MQ=255)
                                                                                                                                 ccggacgcTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgc               >  5:24460/8‑149 (MQ=255)
                                                                                                                                   ggCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACCGCGCACTGGAGAGCGACCAGCAACTCGCTCAGCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgccc             >  5:117404/1‑147 (MQ=255)
                                                                                                                                   ggCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGgcgc             >  3:264771/1‑149 (MQ=255)
                                                                                                                                       gcgaCGCGGTGGTACACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCtgtg         >  5:149277/4‑149 (MQ=255)
                                                                                                                                       gcgaCGCGGTGGTACACCGCGTAGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGa                       >  6:185223/4‑135 (MQ=255)
                                                                                                                                       gcgaCGCGGTGGTACACCGCGTAGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGa                       <  5:185223/132‑1 (MQ=255)
                                                                                                                                        cTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTACTGATGCATGGCGCTGTgg        >  6:192126/1‑149 (MQ=255)
                                                                                                                                         tACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGt                                 >  4:35974/1‑123 (MQ=255)
                                                                                                                                           cgcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTgg        <  3:43570/146‑1 (MQ=255)
                                                                                                                                           cgcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTa     >  1:286952/1‑149 (MQ=255)
                                                                                                                                            gcgGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTAc    >  4:178090/1‑149 (MQ=255)
                                                                                                                                              ggTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCg  >  8:225502/1‑149 (MQ=255)
                                                                                                                                                   accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTa     <  4:120589/141‑1 (MQ=255)
                                                                                                                                                   accaccGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTa     >  3:120589/1‑141 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGTTCTCACGCATTTTGCTCGGGGTGCGTTCCAAGGCCATTCTGGCGCTGCTGTTTTGCGTACTGTCGGCGTTTCTATCCGCGTTTCTCGATGCCCTGACCGTCACGGCGGTGATCATCAGCGCTGCGGTCGGCTTCTACGCGGTGTACCACCGCGTCGCGTCCGGGGCCAATCCGCGCGAGGACAGCGCACTGGAGAGCGACCAGCAAGTCGCTCACCTGCACCGGGAAGACCTCGACCAGTTCCGCGCCTTCCTGCGTAGCCTGCTGATGCATGGCGCTGTGGGTACCG  >  NC_002947/4541633‑4541923

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