Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,651,374 T→G 7.9% 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,3740TG 7.9% 251.7 / 6.9 84intergenic (‑182/+31)PP_4959/fbadiguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base G (4/3);  total (43/41)
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

TTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACG  >  NC_002947/5651230‑5651515
                                                                                                                                                |                                                                                                                                              
ttGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCt                                                                                                                                            >  6:118742/1‑149 (MQ=255)
     tCCATTAGTACTATACACCTCACGCTGCAACACGCAATGCAATCAGAACGATATCAATGATAAAACAGTACGCCAGCCGACTAAAGCTCGGTCTGCAGGAATTATGAACAGCTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCggg                                                                                                                                       >  4:101145/1‑149 (MQ=255)
       cATTAGTACTCTACACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAATAAAACCCGCAGCGATGCGGGtt                                                                                                                                     <  7:55242/149‑1 (MQ=255)
        aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGttt                                                                                                                                    <  1:321339/149‑1 (MQ=255)
        aTTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGttt                                                                                                                                    >  1:95206/1‑149 (MQ=255)
            gTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCa                                                                                                                                <  7:277600/149‑1 (MQ=255)
             tACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGttttttta                                                                                                                               <  6:324364/149‑4 (MQ=255)
              aCTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                               <  2:240840/148‑1 (MQ=255)
               ctctAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAggg                                                                                                                             >  3:313283/1‑149 (MQ=255)
                  tAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATGAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAggg                                                                                                                             <  2:323582/146‑1 (MQ=255)
                  tAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTa                                                                                                                          >  6:232875/1‑149 (MQ=255)
                  tAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTa                                                                                                                          >  5:300967/1‑149 (MQ=255)
                   aGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAc                                                                                                                         <  4:247058/149‑1 (MQ=255)
                    gACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAc                                                                                                                         <  4:313283/148‑1 (MQ=255)
                          aGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCg                                                                                                                   <  5:232875/148‑1 (MQ=255)
                          aGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGa                                                                                                                  <  2:95206/149‑1 (MQ=255)
                            gCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGAt                                                                                                                 >  8:333068/1‑148 (MQ=255)
                              tGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGACGACTAACGGTAGGTCTGCCGTAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAACTCTGCGGGTTTTTCAGGGTAACATGTCGATAAg                                                                                                              >  1:287497/1‑149 (MQ=255)
                                cAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                            >  6:280906/1‑149 (MQ=255)
                                cAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                            <  3:183680/149‑1 (MQ=255)
                                 aaCACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                            >  2:43744/1‑148 (MQ=255)
                                  acacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑a                                                                                                         >  3:133024/1‑149 (MQ=255)
                                   cacGCAATGCAATCAGAACGGTATCAATGATCAAACTCTTCGGTAGCCTACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑Ac                                                                                                        <  7:270476/149‑1 (MQ=255)
                                   cacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑Ac                                                                                                        >  4:31219/1‑149 (MQ=255)
                                   cacGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACCGTAGCCGACTAACGGTCGGTCTGCCGGAATTATCTACAGGATTCCATCACCCATAAAAAAACCCGCATCGCTGCGGGATTTTCAGGGTTACCTGACGATCAGTTG‑a                                                                                                         >  4:224914/1‑148 (MQ=255)
                                        aatgcaatCAGAACGCTATCAATGATCAATCTGTACGGCAGCAGACTAACGGTCGGTCTGCCTGAATTATGTACATGTTTCCATGTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTgg                                                                                                   <  3:101145/149‑1 (MQ=255)
                                            caatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTggcggc                                                                                               <  5:118742/149‑1 (MQ=255)
                                             aatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTggcggc                                                                                               >  4:145657/1‑148 (MQ=255)
                                                cAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGt                                                                                           >  2:213509/1‑149 (MQ=255)
                                                   aaggCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCTTCTGGCATAAAAAAAGCCGCTTGGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTc                                                                                         <  3:224914/145‑1 (MQ=255)
                                                     cGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCAcc                                                                                      >  2:198896/1‑149 (MQ=255)
                                                        tATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCttt                                                                                   >  5:236406/1‑149 (MQ=255)
                                                        tATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCttt                                                                                   >  4:208369/1‑149 (MQ=255)
                                                        tATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCtt                                                                                    >  5:22086/1‑148 (MQ=255)
                                                          tCAATGATCAAACTGTACGGTAGCCGACTAACGGTGGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCTTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTc                                                                                  <  5:283733/148‑1 (MQ=255)
                                                          tCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTc                                                                                  <  5:280906/148‑1 (MQ=255)
                                                          tCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCa                                                                                 <  3:31219/149‑1 (MQ=255)
                                                              tGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTa                                                                             <  6:236406/149‑1 (MQ=255)
                                                               gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCAc                                                                                       <  1:43379/138‑1 (MQ=255)
                                                               gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCAc                                                                                       >  2:43379/1‑138 (MQ=255)
                                                               gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCAc                                                                                       <  8:307182/138‑1 (MQ=255)
                                                               gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCAc                                                                                       >  7:307182/1‑138 (MQ=255)
                                                                 tCAAACTGTACGGTAGCCGACTAACGGTGGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAgcg                                                                          >  1:206708/1‑149 (MQ=255)
                                                                 tCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAgcg                                                                          >  3:136666/1‑149 (MQ=255)
                                                                   aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAgcgc                                                                         <  2:206708/148‑1 (MQ=255)
                                                                    aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGATTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCtc                                                                       >  1:271088/1‑148 (MQ=255)
                                                                    aaCTGTAAGGGAGCAGACTAACGTTCTGTCTTCCGGAATTATGTACAGTTTTCCATCTCCCATAAAAAAACGCGCATCGCTGCGGGTTTTTCAGTGTTACCTGTCGTTCAGTCG‑ACATTGGCGGCCATTTCACCTTTCAGGTAGCGCTg                                                                       <  2:287497/149‑1 (MQ=255)
                                                                      cTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                      >  5:336495/1‑148 (MQ=255)
                                                                            ggTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGc                                                               <  6:300967/149‑1 (MQ=255)
                                                                                gCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTCGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                           >  5:316669/1‑149 (MQ=255)
                                                                                  cGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                         <  6:22086/149‑1 (MQ=255)
                                                                                        aCGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagagag                                                    <  7:333068/148‑1 (MQ=255)
                                                                                        aCGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagagag                                                    <  1:213509/148‑1 (MQ=255)
                                                                                        aCGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCagaga                                                     <  3:145657/147‑1 (MQ=255)
                                                                                         cggtcggtCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGAt                                                  <  3:227424/149‑1 (MQ=255)
                                                                                          ggtcggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGCTGCGGGGTTTTCAGGGTTACCTGTCGATCAGTTG‑CCCTTGGCGGCCAGTTCCCCGTTCCGGTAGCGCTGGTACATGCGTGCCagcgagc                                                   >  6:375462/1‑146 (MQ=255)
                                                                                                     cGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCtt                                      >  3:203180/1‑149 (MQ=255)
                                                                                                         aTTATGTACTGGTTTCCATCTCCGATAAAAAAACCCGCATGGCTGCGGGTTTTTCAGGGTTTCCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGTGAGATCGGCTTGATCTTCGa                                   <  2:271088/148‑1 (MQ=255)
                                                                                                          ttATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCgg                                               >  3:19191/1‑135 (MQ=255)
                                                                                                          ttATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCgg                                               <  4:19191/135‑1 (MQ=255)
                                                                                                            aTGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGc                                <  4:133024/148‑1 (MQ=255)
                                                                                                              gTACAGGTTTCCATCTCACATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑CCCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGCTCGGCTTGATCTTCGAGGCAtg                             >  6:47430/1‑148 (MQ=255)
                                                                                                               tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCAtt                             >  5:315030/1‑148 (MQ=255)
                                                                                                               tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCAtt                             <  6:315030/148‑1 (MQ=255)
                                                                                                                aCAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCTGCTTTTTCAGGGTTACCTATCGATCAGTTG‑ACCATGGCGGCCAGTTAACCTTTGAGGTAGCGCTGGTACATGCCATACAGATAGATCGGCATCATCTTAGAGGCATTg                            >  3:277883/1‑148 (MQ=255)
                                                                                                                 cAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                           <  1:198896/148‑1 (MQ=255)
                                                                                                                        ccATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGCACCTTGGCGGCCAGTTCACCTTTCAgg                                                                               <  1:16096/90‑1 (MQ=255)
                                                                                                                        ccATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAgg                                                                               >  2:16096/1‑89 (MQ=255)
                                                                                                                        ccATCTCCCATAAAAAAACCCGCATCGCGGCGGGGTTTTCAGGGTTACCTGTGGCTCAGTTG‑ACCTTGGCGGCCAGTTCACCGTTCAGGTAGCGCTGGGACATGCCTTCCGGAGAGAGCGGCTGGATCTTCGAGGCATTGCCCGCGGTg                   >  6:140042/1‑149 (MQ=255)
                                                                                                                         cATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATc                                                 <  8:34244/118‑1 (MQ=255)
                                                                                                                         cATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATc                                                 >  7:34244/1‑118 (MQ=255)
                                                                                                                         cATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                  <  1:43744/149‑1 (MQ=255)
                                                                                                                         cATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGc                  >  4:361990/1‑149 (MQ=255)
                                                                                                                          aTCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGcc                 >  5:49066/1‑149 (MQ=255)
                                                                                                                               ccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg            >  5:370611/1‑149 (MQ=255)
                                                                                                                               ccATAAAAAAACCCGCAGCCCTGCGGGTTTTTCAGGGTTACCTGTCGATAAGTTG‑ACCTTGGCGGCCAGTTCACCTTTGAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATATTCGAGGCATTGCCCGCGGTGCCGAAgg            >  2:103441/1‑149 (MQ=255)
                                                                                                                                 aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTcc                                                          <  6:360224/101‑1 (MQ=255)
                                                                                                                                 aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTAGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTcc                                                          >  5:360224/1‑101 (MQ=255)
                                                                                                                                  tAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCtt         >  8:286792/1‑149 (MQ=255)
                                                                                                                                   tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGaa              >  7:144743/3‑143 (MQ=255)
                                                                                                                                   tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGaa              <  8:144743/141‑1 (MQ=255)
                                                                                                                                    aaaaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg       <  6:370611/149‑1 (MQ=255)
                                                                                                                                        aaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCAGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAc   <  4:129621/149‑1 (MQ=255)
                                                                                                                                         aCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcg  >  3:359083/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TTGGATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTG‑ACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACG  >  NC_002947/5651230‑5651515

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