Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,651,374 T→G 13.1% intergenic (‑182/+31) PP_4959 ← / ← fba diguanylate cyclase/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,651,3740TG13.1% 259.9 / 27.9 98intergenic (‑182/+31)PP_4959/fbadiguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (45/40);  new base G (7/6);  total (52/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 = 1.00e+00

CGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGC  >  NC_002947/5651227‑5651518
                                                                                                                                                   |                                                                                                                                                
cGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAt                                                                                                                                                  <  1:330380/148‑1 (MQ=255)
   ttGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCt                                                                                                                                              >  4:197059/1‑149 (MQ=255)
   ttGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCgat                                                                                                                                              >  8:269948/1‑147 (MQ=255)
       aTCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                          >  2:249236/1‑149 (MQ=255)
          cATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGt                                                                                                                                        <  3:257809/148‑1 (MQ=255)
          cATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtt                                                                                                                                       <  4:370233/149‑1 (MQ=255)
           aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttt                                                                                                                                      <  7:269948/149‑1 (MQ=255)
              aGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttttt                                                                                                                                    <  7:143828/148‑1 (MQ=255)
               gTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATTATCAAACTGTAAGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCa                                                                                                                                  <  2:196415/149‑1 (MQ=255)
                tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTACCGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                 >  4:313690/1‑149 (MQ=255)
                tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                 >  3:154311/1‑149 (MQ=255)
                tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttttttta                                                                                                                                 <  4:286703/149‑4 (MQ=255)
                tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttttt                                                                                                                                    <  3:59757/146‑1 (MQ=255)
                tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttttt                                                                                                                                    >  4:59757/1‑146 (MQ=255)
                  ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAggg                                                                                                                               >  8:108581/1‑149 (MQ=255)
                  ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAgg                                                                                                                                <  3:299375/148‑1 (MQ=255)
                    ctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCAGCATCGCTGCGGGTTTTTCAGGGt                                                                                                                              >  7:209892/1‑148 (MQ=255)
                       gACCTGAGGCTGCAACACGCAATGCACTCAGCACGCTATCAATGATCAACCTGTACGGTAGCCGTCTAACGGTCGGTCTGCCGTAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAc                                                                                                                           <  8:345315/148‑1 (MQ=255)
                       gACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTa                                                                                                                            <  5:122770/147‑1 (MQ=255)
                       gACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAcc                                                                                                                          >  4:163432/1‑149 (MQ=255)
                          cTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGt                                                                                                                       >  6:240627/1‑149 (MQ=255)
                            gAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTc                                                                                                                      >  3:16787/1‑148 (MQ=255)
                                cTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCa                                                                                                                 <  7:2274/149‑1 (MQ=255)
                                    atcacGCAATGCAATCAGAACGCTATCAATGATCATACTGTTCGGTAGCCGACTATCTGTCGGTCTGCCGGAATTATTTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGTGTTTTTCAGGGTTa                                                                                                                            <  4:230812/132‑1 (MQ=255)
                                    aaCACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                              >  8:189956/1‑148 (MQ=255)
                                    aaCACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                              >  4:318415/1‑148 (MQ=255)
                                     acacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtttttttat                                                                                                                                <  4:79063/129‑5 (MQ=255)
                                     acacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtttttttat                                                                                                                                >  3:79063/1‑125 (MQ=255)
                                       acGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGAcc                                                                                                          >  5:157580/1‑149 (MQ=255)
                                             tgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTCCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTggcg                                                                                                    >  3:303490/1‑149 (MQ=255)
                                                aatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTggcggc                                                                                                  >  4:297693/1‑148 (MQ=255)
                                                aatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTggcggc                                                                                                  >  2:40710/1‑148 (MQ=255)
                                                    aGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGtt                                                                                             >  8:71666/1‑149 (MQ=255)
                                                     gAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGtt                                                                                             >  6:275556/1‑148 (MQ=255)
                                                     gAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTc                                                                                            <  3:318415/149‑1 (MQ=255)
                                                       aCGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCAc                                                                                          <  7:108581/149‑1 (MQ=255)
                                                       aCGCTATCAATGATCAAACTGAACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAACACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCAc                                                                                          <  8:177274/149‑1 (MQ=255)
                                                           tATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCttt                                                                                      >  8:259744/1‑149 (MQ=255)
                                                             tCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCa                                                                                    <  7:71666/149‑1 (MQ=255)
                                                               aaTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                  <  1:249236/149‑1 (MQ=255)
                                                                aTGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                  >  1:10935/1‑148 (MQ=255)
                                                                 tGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAgg                                                                                  >  7:226062/1‑147 (MQ=255)
                                                                      aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCt                                                                           >  3:221909/1‑149 (MQ=255)
                                                                       aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTg                                                                          >  4:263480/1‑149 (MQ=255)
                                                                       aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACATGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTg                                                                          <  1:87262/149‑1 (MQ=255)
                                                                            tACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACa                                                                     >  8:43521/1‑149 (MQ=255)
                                                                              cGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACAt                                                                    <  2:10935/148‑1 (MQ=255)
                                                                              cGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACAt                                                                    <  6:157580/148‑1 (MQ=255)
                                                                                 tAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCt                                                                <  4:303490/149‑1 (MQ=255)
                                                                                 tAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCt                                                                <  7:189956/149‑1 (MQ=255)
                                                                                     cGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                            <  7:43521/149‑1 (MQ=255)
                                                                                      gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTTACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGATGGTACATGCATTc                                                              >  7:18230/1‑146 (MQ=255)
                                                                                       aCTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                          <  3:225694/149‑1 (MQ=255)
                                                                                       aCTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                            >  7:238827/1‑147 (MQ=255)
                                                                                        ccgttagtcggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagag                                                         <  2:103397/143‑1 (MQ=255)
                                                                                        cTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagag                                                         <  3:297693/149‑1 (MQ=255)
                                                                                            cggtcggtCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAATAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                          >  1:262040/1‑144 (MQ=255)
                                                                                            cggtcggtCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCaga                                                          >  7:178103/1‑144 (MQ=255)
                                                                                                  gtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCtt                                               <  3:197059/149‑1 (MQ=255)
                                                                                                     tGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCa                                                                                    <  2:151723/109‑1 (MQ=255)
                                                                                                         ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACAt                                                                    <  1:246323/121‑1 (MQ=255)
                                                                                                         ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACAt                                                                    >  2:246323/1‑121 (MQ=255)
                                                                                                         ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                        <  3:263480/149‑1 (MQ=255)
                                                                                                         ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                        <  4:16787/149‑1 (MQ=255)
                                                                                                            attattTAAATGTTTTCATATCCCATAAAAAAAAAATCATCGCTTCGGTTTTTTCAGGGTTACGTGTCGATCAGTTGACCTTGGCGGTCATTTCAAATTTCAGTTAGCGCTGGTACATGACTTCCAGAGAGATAGGCTTTATCTTCGa                                      <  8:18230/148‑1 (MQ=255)
                                                                                                            aTTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAg                                     <  8:209892/149‑1 (MQ=255)
                                                                                                                 gTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTg                                              <  4:208937/135‑1 (MQ=255)
                                                                                                                 gTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTg                                              >  3:208937/1‑135 (MQ=255)
                                                                                                                  tACCGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCa                                  <  8:66980/146‑1 (MQ=255)
                                                                                                                  tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTg                               >  2:42881/1‑149 (MQ=255)
                                                                                                                  tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTg                               >  6:127806/1‑149 (MQ=255)
                                                                                                                     aGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGccc                            >  7:70121/1‑149 (MQ=255)
                                                                                                                      ggTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCcg                           <  2:100230/149‑1 (MQ=255)
                                                                                                                       gTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCcgc                          >  4:265685/1‑149 (MQ=255)
                                                                                                                       gTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCcgc                          >  2:203428/1‑149 (MQ=255)
                                                                                                                        tttCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCcgcg                         <  3:265685/149‑1 (MQ=255)
                                                                                                                           ccATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGt                       <  4:221909/148‑1 (MQ=255)
                                                                                                                            cATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                     <  7:259744/149‑1 (MQ=255)
                                                                                                                            cATCTCCCAT‑AAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                         >  3:81374/1‑96 (MQ=255)
                                                                                                                            cATCTCCCAT‑AAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                         <  4:81374/96‑1 (MQ=255)
                                                                                                                              tctcCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGccc                            >  2:127139/1‑140 (MQ=255)
                                                                                                                              tctcCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGccc                            <  1:127139/140‑1 (MQ=255)
                                                                                                                                tcCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCg                   >  7:201088/1‑147 (MQ=255)
                                                                                                                                  ccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGaa                 <  8:238827/147‑1 (MQ=255)
                                                                                                                                  ccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg               >  2:256115/1‑149 (MQ=255)
                                                                                                                                   cATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                         >  2:322294/1‑90 (MQ=255)
                                                                                                                                      tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc           >  6:372009/3‑149 (MQ=255)
                                                                                                                                      tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc           >  8:160428/3‑149 (MQ=255)
                                                                                                                                      tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc           >  4:340392/3‑149 (MQ=255)
                                                                                                                                      tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTc           >  3:127457/3‑149 (MQ=255)
                                                                                                                                       aaaaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTTATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg          >  6:304558/1‑149 (MQ=255)
                                                                                                                                       aaaaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg          <  3:313690/149‑1 (MQ=255)
                                                                                                                                       aaaaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg          >  5:23336/1‑149 (MQ=255)
                                                                                                                                         aaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTAAATGCCTTCAAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGt         >  7:175147/1‑148 (MQ=255)
                                                                                                                                          aaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGTGGCCAGTTCACCTTTCAGGTGGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTaa       >  8:289587/1‑149 (MQ=255)
                                                                                                                                          aaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTaa       >  2:325342/1‑149 (MQ=255)
                                                                                                                                           aaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAc      <  8:201088/149‑1 (MQ=255)
                                                                                                                                              ccGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGTATGCTTCGTAAcgcg   <  8:70121/149‑1 (MQ=255)
                                                                                                                                                gCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc  >  3:137504/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGGTTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGC  >  NC_002947/5651227‑5651518

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