Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 326,065 C→A 7.2% *556Y (TAG→TAT PP_0269 ← glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947326,0650CA 7.2% 228.5 / 2.0 75*556Y (TAG→TATPP_0269glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base A (3/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAACTCACC  >  NC_002947/325917‑326211
                                                                                                                                                    |                                                                                                                                                  
gTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGcc                                                                                                                                                    >  1:21312/1‑149 (MQ=255)
   gATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTa                                                                                                                                                  <  8:48512/148‑1 (MQ=255)
    aTCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATg                                                                                                                                                >  6:55806/1‑149 (MQ=255)
        aGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCgg                                                                                                                                            <  8:149448/149‑1 (MQ=255)
         gCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGCATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCggg                                                                                                                                           >  5:114446/1‑149 (MQ=255)
           tGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGc                                                                                                                                         >  1:145086/1‑149 (MQ=255)
           tGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGc                                                                                                                                         >  1:139652/1‑149 (MQ=255)
           tGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGc                                                                                                                                         >  8:255088/1‑149 (MQ=255)
            gCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCt                                                                                                                                        <  7:255088/149‑1 (MQ=255)
             cGTGGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCt                                                                                                                                        <  4:79173/148‑1 (MQ=255)
              gTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGt                                                                                                                                      >  7:6898/1‑149 (MQ=255)
                 gCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGttt                                                                                                                                    <  5:55806/148‑1 (MQ=255)
                  cTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCAGCGTATTAGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCt                                                                                                                                  >  3:290421/1‑149 (MQ=255)
                     gTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCATAGGCCGGGGCTGttcgata                                                                                                                                <  8:355221/148‑6 (MQ=255)
                      ttCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTc                                                                                                                              <  8:16295/149‑1 (MQ=255)
                        cTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCgggg                                                                                                                                          <  2:131022/135‑1 (MQ=255)
                        cTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCg                                                                                                                             <  8:227025/148‑1 (MQ=255)
                        cTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCg                                                                                                                             <  3:90045/148‑1 (MQ=255)
                             gTACCAACAACTCCGTGCATCAGAACGTTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCTGGCTTATGTCGGGGCTGTTTCTTTTCGGCGAAc                                                                                                                       <  8:56407/149‑1 (MQ=255)
                               aCCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCgggg                                                                                                                                          >  3:305510/1‑128 (MQ=255)
                               aCCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCgggg                                                                                                                                          <  4:305510/128‑1 (MQ=255)
                               aCCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACgg                                                                                                                     >  8:171598/1‑149 (MQ=255)
                                ccAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGCATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGt                                                                                                                    <  6:114446/149‑1 (MQ=255)
                                   acaacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTAc                                                                                                                 >  6:311818/1‑149 (MQ=255)
                                    caacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACt                                                                                                                <  6:149609/149‑1 (MQ=255)
                                     aacGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTg                                                                                                               >  3:162178/1‑149 (MQ=255)
                                      acGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGt                                                                                                              <  8:263795/149‑1 (MQ=255)
                                            tGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGc                                                                                                        <  2:139652/149‑1 (MQ=255)
                                            tGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGc                                                                                                        <  2:145086/149‑1 (MQ=255)
                                                    aaCGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTg                                                                                                >  2:178636/1‑149 (MQ=255)
                                                     aCGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTgg                                                                                               <  2:21312/149‑1 (MQ=255)
                                                           aCAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTa                                                                                          <  8:6898/148‑1 (MQ=255)
                                                              acgacgacgGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGt                                                                                      >  4:326175/1‑149 (MQ=255)
                                                                    acgGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCgg                                                                                <  1:178636/149‑1 (MQ=255)
                                                                     cgGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGt                                                                               <  8:247632/149‑1 (MQ=255)
                                                                          actaattCCTCGTATTCTTCGAGTACCAGATCAGCATCTTCTTATTCTGATCTGATGTATCGAACAGCCCCGGCATATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCaaa                                                                          <  4:290421/143‑1 (MQ=255)
                                                                                  cgcgTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACg                                                                                        >  4:209082/1‑127 (MQ=255)
                                                                                  cgcgTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACg                                                                                        <  3:209082/127‑1 (MQ=255)
                                                                                           gtcgAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCa                                                         >  8:106157/1‑149 (MQ=255)
                                                                                                 tACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCtt                                                                                           >  4:49975/1‑109 (MQ=255)
                                                                                                 tACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGcc                                                   <  7:171598/149‑1 (MQ=255)
                                                                                                 tACCCGATCAGCATCTTCTGATTCTGATCTGAAGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGc                                                    >  6:319181/1‑148 (MQ=255)
                                                                                                   cccGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTGTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCca                                                  >  5:18480/1‑148 (MQ=255)
                                                                                                   cccGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCt                                                                       >  6:265607/1‑127 (MQ=255)
                                                                                                   cccGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCt                                                                       <  5:265607/127‑1 (MQ=255)
                                                                                                          aGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCcgcg                                          >  5:280922/1‑149 (MQ=255)
                                                                                                            cATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCa                                        >  8:30739/1‑149 (MQ=255)
                                                                                                             aTCTTCTGATTCTTATATGATGTATCGAACAGCCCCGGCCTATTCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGa                                                                                                          <  6:327651/82‑1 (MQ=255)
                                                                                                             aTCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTAACTGTTTGa                                                                                                          >  5:327651/1‑82 (MQ=255)
                                                                                                                ttctgattctgaTCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATaa                                     >  5:177265/1‑148 (MQ=255)
                                                                                                                  ctgattctgatCTGATGTATCGAACAGCCCCGGCTTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATaaaa                                   <  4:162178/148‑1 (MQ=255)
                                                                                                                   tgattctgatCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACt                                 >  5:50501/1‑149 (MQ=255)
                                                                                                                      ttctgatCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGc                              >  6:322181/1‑149 (MQ=255)
                                                                                                                      ttctgatCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCAGAAAACTCGc                              >  1:26774/1‑149 (MQ=255)
                                                                                                                       tctgatctgaTGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCt                             >  7:153914/1‑149 (MQ=255)
                                                                                                                           atctgatGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCttt                           <  6:280922/147‑1 (MQ=255)
                                                                                                                           atctgatGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCttt                           <  6:177265/147‑1 (MQ=255)
                                                                                                                               gatGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGTTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAg                                                                         <  5:172376/97‑1 (MQ=255)
                                                                                                                               gatGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAg                                                                         >  6:172376/1‑97 (MQ=255)
                                                                                                                                   tagcGAACAGCCCCGGCCAATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCttt                           <  8:239101/136‑1 (MQ=255)
                                                                                                                                   tATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTGTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGGGCATAAAACTCGCttt                           >  7:239101/1‑139 (MQ=255)
                                                                                                                                   tATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGc                  <  7:106157/148‑1 (MQ=255)
                                                                                                                                       gAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACt              >  1:347597/1‑148 (MQ=255)
                                                                                                                                       gAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACt              >  3:8510/1‑148 (MQ=255)
                                                                                                                                        aaCAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACt              >  7:185360/1‑147 (MQ=255)
                                                                                                                                         aCAGCCCCGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAg           >  5:29180/1‑149 (MQ=255)
                                                                                                                                         aCAGCCCCGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAg           >  5:366213/1‑149 (MQ=255)
                                                                                                                                          cAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGt          <  5:311818/149‑1 (MQ=255)
                                                                                                                                           aGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTa         <  3:326175/149‑1 (MQ=255)
                                                                                                                                            gCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTaa        <  1:85319/149‑1 (MQ=255)
                                                                                                                                            gCCCCGGCATAGGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTaa        >  6:310057/1‑149 (MQ=255)
                                                                                                                                             ccccGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTaa        >  8:20123/1‑148 (MQ=255)
                                                                                                                                             ccccGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAAc       >  8:25660/1‑149 (MQ=255)
                                                                                                                                               ccGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCa                                                                            >  7:128543/1‑78 (MQ=255)
                                                                                                                                               ccGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCa                                                                            <  8:128543/78‑1 (MQ=255)
                                                                                                                                                  gCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAACTCAcc  <  8:153914/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTCGATCAAGCTGCGTAGCTCGTTCTTGCGTACCAACAACGCCGTGCAGCAGAACGGTTACAACGACGACGGCAACGAAGTCCGCGTATTCGTCGAGTACCCGATCAGCATCTTCTGATTCTGATCTGATGTATCGAACAGCCCCGGCCTATGCCGGGGCTGTTTCTTTTCGGCGAACGGTTACTGTTTGAGCCGGCGCTGGGCTTACGGTCTCCGGTGCAAAGCTGTCGCTACGCGCCATCCGCCACATCCGCGCATAAAACTCGCTTTCGATCGAGCCACTGAGTAACTCACC  >  NC_002947/325917‑326211

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: