Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,535,425 G→C 5.5% E148Q (GAG→CAG)  pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,535,4250GC 5.5% 281.1 / 4.3 87E148Q (GAG→CAG) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (39/43);  new base C (2/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

CGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTGGCAAGGG  >  NC_002947/5535278‑5535563
                                                                                                                                                      |                                                                                                                                          
cGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGa                                                                                                                                           <  7:51614/149‑1 (MQ=255)
      gTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCTCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTCGGCCAGgcgca                                                                                                                                     <  6:535/149‑2 (MQ=255)
      gTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGcgcg                                                                                                                                     >  5:157297/1‑149 (MQ=255)
       tGCACGTCGGCTTCGAACATGATCAGCTGATCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCGGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCgg                                                                                                                                    <  6:225953/149‑1 (MQ=255)
            gTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGAc                                                                                                                                 >  6:240486/1‑147 (MQ=255)
            gTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCa                                                                                                                               <  1:406191/149‑1 (MQ=255)
             tCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAg                                                                                                                              <  6:242841/149‑1 (MQ=255)
             tCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAg                                                                                                                              <  3:358990/149‑1 (MQ=255)
             tCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGACTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAg                                                                                                                              <  2:2209/149‑1 (MQ=255)
              cGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATAATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGt                                                                                                                             >  2:361740/1‑149 (MQ=255)
                gCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCCTGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACCTGCGCCTGGCCGAGGCGCGGGACCAGTTg                                                                                                                           <  3:281182/149‑1 (MQ=255)
                  aTCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCTCATGCGCCTGGCCGAGGCGCGGGACCAGTTGTa                                                                                                                         <  1:52928/149‑1 (MQ=255)
                  aTCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTa                                                                                                                         >  4:35600/1‑149 (MQ=255)
                  aTCGAAAATAATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTa                                                                                                                         <  7:147135/149‑1 (MQ=255)
                     gAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCg                                                                                                                      >  2:232117/1‑149 (MQ=255)
                      aaaaTGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCg                                                                                                                      >  1:196922/1‑148 (MQ=255)
                      aaaaTGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCg                                                                                                                      >  2:208707/1‑148 (MQ=255)
                        aaTGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCaa                                                                                                                   <  3:35600/149‑1 (MQ=255)
                              cAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGa                                                                                                             >  3:57638/1‑149 (MQ=255)
                               aGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGAc                                                                                                            <  2:72554/149‑1 (MQ=255)
                                gCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACt                                                                                                           >  5:293079/1‑149 (MQ=255)
                                 cTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTg                                                                                                          >  8:287767/1‑149 (MQ=255)
                                    ttCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTgg                                                                                                         <  5:262598/147‑1 (MQ=255)
                                    ttCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTgg                                                                                                         <  6:7101/147‑1 (MQ=255)
                                    ttCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGc                                                                                                        <  1:208707/148‑1 (MQ=255)
                                       cTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCtgc                                                                                                     <  2:156121/148‑1 (MQ=255)
                                           gcCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCtgcgtgcg                                                                                                >  8:109113/1‑149 (MQ=255)
                                           gcCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCtgc                                                                                                     <  8:323397/144‑1 (MQ=255)
                                           gcCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCtgc                                                                                                     >  7:323397/1‑144 (MQ=255)
                                            ccTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTCGGCCAGgcgcat                                                                                                                                    <  3:373369/112‑3 (MQ=255)
                                            ccTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTCGGCCAGgcgcat                                                                                                                                    >  4:373369/1‑110 (MQ=255)
                                             cTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCa                                                                                              >  6:55513/1‑149 (MQ=255)
                                             cTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGAGTGCGCa                                                                                              >  6:61845/1‑149 (MQ=255)
                                                  aGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCCCATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGc                                                                                           >  4:39584/1‑147 (MQ=255)
                                                  aGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCa                                                                                         <  5:348011/149‑1 (MQ=255)
                                                  aGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCa                                                                                         <  5:240486/149‑1 (MQ=255)
                                                     tcggTGGCTGGATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATGATCTGCGAAGTACGCAAGCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCTGGACCATTTGTACCGCATGTTCGACTGGCTGCGTGCGCATGCCAGcg                                                                                      <  6:55106/146‑1 (MQ=255)
                                                       cGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGcgac                                                                                    >  6:418806/1‑149 (MQ=255)
                                                          tggcgCATGTGCTTCCTTTTCGATGTGCCGCTGCGGGCATCATTATCGCATGTAGCCAACGCCATCTTGCCCCCCCAATTGGGCCTGGCGCAGGGGGGGTACCTGTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGcgacgac                                                                                 <  4:19125/145‑1 (MQ=255)
                                                             tGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAg                                                                              >  7:61732/1‑149 (MQ=255)
                                                               cATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAg                                                                              <  7:109113/147‑1 (MQ=255)
                                                                         cTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAAc                                                                    <  7:147520/148‑1 (MQ=255)
                                                                              cGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAg                                                               <  8:267960/148‑1 (MQ=255)
                                                                               gAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTgg                                                                                                         >  4:93569/1‑105 (MQ=255)
                                                                               gAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTgg                                                                                                         <  3:93569/105‑1 (MQ=255)
                                                                               gAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGt                                                             >  1:412755/1‑149 (MQ=255)
                                                                                 aGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTgg                                                           <  8:61732/149‑1 (MQ=255)
                                                                                    gCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGcc                                                         <  4:57638/148‑1 (MQ=255)
                                                                                     ccGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTgg                                                           <  7:343338/145‑1 (MQ=255)
                                                                                     ccGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTgg                                                           >  8:343338/1‑145 (MQ=255)
                                                                                       gctgctGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGCCt                                                     >  4:360582/1‑149 (MQ=255)
                                                                                         tgctgGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGt                                                             >  4:250893/1‑139 (MQ=255)
                                                                                         tgctgGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGt                                                             <  3:250893/139‑1 (MQ=255)
                                                                                           ctgGCCATCATCAGCGAAGTACTCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCg                                                  <  5:61845/148‑1 (MQ=255)
                                                                                           ctgGCCATCATCAGCGAAGTACGCTACCGCAATGTGCACCCGCACATGCGCCTGGCCGAGGGGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGGGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTAGCCGACTTCg                                                  <  1:361740/148‑1 (MQ=255)
                                                                                           ctgGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCg                                                  <  5:55513/148‑1 (MQ=255)
                                                                                           ctgGCAATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCgg                                                 <  2:196922/149‑1 (MQ=255)
                                                                                             ggCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCa                                               >  5:221678/1‑149 (MQ=255)
                                                                                             ggCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCa                                               >  5:104208/1‑149 (MQ=255)
                                                                                             ggCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCa                                               >  3:232968/1‑149 (MQ=255)
                                                                                                catcatCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCAc                                              >  1:17265/1‑147 (MQ=255)
                                                                                                catcatCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCAc                                              <  2:17265/147‑1 (MQ=255)
                                                                                                      caGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGgc                                      >  8:232055/1‑149 (MQ=255)
                                                                                                          gAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCtt                                  <  2:377697/149‑1 (MQ=255)
                                                                                                          gAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCtt                                  >  3:11046/1‑149 (MQ=255)
                                                                                                                 cGCAACCGAAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCa                             <  1:232117/148‑1 (MQ=255)
                                                                                                                  gCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTCCGGCACCCGCCGGCGCTTTTCCAGc                           >  5:190313/1‑149 (MQ=255)
                                                                                                                  gCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACCTCGGCACCCGCCGGCGCTTTTCCAGc                           >  5:54435/1‑149 (MQ=255)
                                                                                                                    aaCCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCcg                         >  3:394870/1‑149 (MQ=255)
                                                                                                                     aCCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCcgc                        <  1:263540/149‑1 (MQ=255)
                                                                                                                        gCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTg                     <  7:297612/149‑1 (MQ=255)
                                                                                                                          aaTCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTg                     >  4:185400/1‑147 (MQ=255)
                                                                                                                          aaTCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCa                   <  5:418806/149‑1 (MQ=255)
                                                                                                                           aTCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCag                  <  6:293079/149‑1 (MQ=255)
                                                                                                                               tGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGAATGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCaggag               >  2:304852/1‑149 (MQ=255)
                                                                                                                                 cACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACCTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCaggagg              <  6:54435/148‑1 (MQ=255)
                                                                                                                                  aCCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCaggagg              >  5:56547/1‑147 (MQ=255)
                                                                                                                                  aCCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCaggagg              >  3:12889/1‑147 (MQ=255)
                                                                                                                                    ccGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCtt                                  >  7:261262/1‑125 (MQ=255)
                                                                                                                                    ccGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCtt                                  <  8:261262/125‑1 (MQ=255)
                                                                                                                                    ccGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGt          >  7:39246/1‑149 (MQ=255)
                                                                                                                                    ccGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGt          >  7:226026/1‑149 (MQ=255)
                                                                                                                                     tggtcccgcgcCTCGGCCAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTg         >  3:419767/8‑149 (MQ=255)
                                                                                                                                       cacaTGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTGGc       <  3:39584/149‑1 (MQ=255)
                                                                                                                                       cacaTGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTCCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTGGc       <  6:190313/149‑1 (MQ=255)
                                                                                                                                           tGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTGGCAAgg   <  7:30327/149‑1 (MQ=255)
                                                                                                                                            gcgcCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTGGCAAggg  >  4:165241/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
CGTTACGTGCACGTCGGCATCGAAAATGATCAGCTGTTCCTGCGCCTGAAAGGTCCGTGGCTGCATGTGATCCTGTTCGAAGTGCCGCTGCTGGCCATCATCAGCGAAGTACGCAACCGCAATCTGCACCCGCACATGCGCCTGGCCGAGGCGCGGGACCAGTTGTACCGCAAGTTCGACTGGCTGCGTGCGCATGCCAGCGACGACGAGCTGGCCGAACTGCAGGTGGCCGACTTCGGCACCCGCCGGCGCTTTTCCAGCCGCGTGCAGGAGGAAGTGGCAAGGG  >  NC_002947/5535278‑5535563

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: