Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 327,335 G→C 5.5% P133R (CCG→CGG)  PP_0269 ← glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947327,3350GC 5.5% 237.9 / 3.7 73P133R (CCG→CGG) PP_0269glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base C (2/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01

CGCCGCGGTTGAGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACC  >  NC_002947/327189‑327480
                                                                                                                                                    |                                                                                                                                                 
cgccgcGGTTGAGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGt                                                                                                                                                 <  1:253838/149‑1 (MQ=255)
 gccgcGGTTGAGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCgg                                                                                                                                                  >  4:281983/1‑147 (MQ=255)
 gccgcGGTTGAGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTa                                                                                                                                                <  5:47380/149‑1 (MQ=255)
           aGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAAt                                                                                                                                       <  6:218389/148‑1 (MQ=255)
           aGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCCGTAGCATGGAATg                                                                                                                                      >  4:133152/1‑149 (MQ=255)
           aGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCCGTAGCATGGAATg                                                                                                                                      >  4:127558/1‑149 (MQ=255)
            gcgcGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATgg                                                                                                                                     <  1:18932/149‑1 (MQ=255)
            gcgcGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATgg                                                                                                                                     >  2:311572/1‑149 (MQ=255)
            gcgcGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATgg                                                                                                                                     >  7:134102/1‑149 (MQ=255)
             cgcgGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGc                                                                                                                                    <  7:188750/149‑1 (MQ=255)
              gcggcgATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCt                                                                                                                                   >  7:164212/1‑149 (MQ=255)
                ggcgATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCtga                                                                                                                                 <  3:203890/149‑1 (MQ=255)
                       gCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTc                                                                                                                          <  8:335931/149‑1 (MQ=255)
                          ttGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGTGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCTTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGaa                                                                                                                       <  2:66237/149‑1 (MQ=255)
                              cgcTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGTCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCa                                                                                                                   <  5:127041/149‑1 (MQ=255)
                               gcTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAg                                                                                                                  >  3:50064/1‑149 (MQ=255)
                                    gggCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTg                                                                                                             >  2:94454/1‑149 (MQ=255)
                                        gcCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACGACCAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTa                                                                                                           <  1:111354/149‑1 (MQ=255)
                                        gcCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAgg                                                                                                         >  3:136546/1‑149 (MQ=255)
                                           gAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGcgt                                                                                                      >  4:130939/1‑149 (MQ=255)
                                               cTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGAATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGtt                                                                                                  <  4:293136/149‑1 (MQ=255)
                                                  aTCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAg                                                                                               >  8:327213/1‑149 (MQ=255)
                                                   tCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGc                                                                                              <  8:164212/149‑1 (MQ=255)
                                                   tCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGc                                                                                              >  7:366201/1‑149 (MQ=255)
                                                       cGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGc                                                                                          >  6:147255/1‑149 (MQ=255)
                                                       cGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGc                                                                                          <  3:330368/149‑1 (MQ=255)
                                                             gTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACg                                                                                    >  5:176368/1‑149 (MQ=255)
                                                                gAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCtt                                                                                 >  3:138906/1‑149 (MQ=255)
                                                                gAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGAGGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGTACGCtt                                                                                 >  7:140325/1‑149 (MQ=255)
                                                                 aaGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCttt                                                                                <  6:332636/149‑1 (MQ=255)
                                                                 aaGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGGCGTTGAGCGTGCCGAACGCtt                                                                                 <  8:134102/148‑1 (MQ=255)
                                                                   gATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTtct                                                                              <  8:366201/149‑1 (MQ=255)
                                                                   gATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTtct                                                                              >  1:222279/1‑149 (MQ=255)
                                                                     tCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGc                                                                                              <  6:175588/131‑1 (MQ=255)
                                                                     tCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTtctc                                                                             <  7:327213/148‑1 (MQ=255)
                                                                     tCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCCGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCg                                                                            <  3:133152/149‑1 (MQ=255)
                                                                     tCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCCGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCg                                                                            <  3:127558/149‑1 (MQ=255)
                                                                        aGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGgcg                                                                         >  6:51999/1‑149 (MQ=255)
                                                                             gAGTAGGGCATTCGGCTTTGCGGCCCTCCGATGGCTATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTc                                                                    <  4:136546/149‑1 (MQ=255)
                                                                               gTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCg                                                                   <  1:94454/148‑1 (MQ=255)
                                                                                   ggCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCtt                                                              >  7:225560/1‑149 (MQ=255)
                                                                                      aTTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGATTGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGgc                                                           >  2:67853/1‑149 (MQ=255)
                                                                                         cGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGcgcg                                                        >  3:64889/1‑149 (MQ=255)
                                                                                          ggCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAg                                                                                               <  6:342630/109‑1 (MQ=255)
                                                                                          ggCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAg                                                                                               >  5:342630/1‑109 (MQ=255)
                                                                                          ggCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCgg                                                       <  2:254961/149‑1 (MQ=255)
                                                                                           gCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCggg                                                      >  5:105337/1‑149 (MQ=255)
                                                                                                  ggCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGAc                                               <  5:147255/149‑1 (MQ=255)
                                                                                                   gCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGAcc                                              <  8:225560/149‑1 (MQ=255)
                                                                                                      cgccgATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCagc                                           >  1:188670/1‑149 (MQ=255)
                                                                                                         cgATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGc                                                     <  6:208220/136‑1 (MQ=255)
                                                                                                         cgATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGc                                                     >  5:208220/1‑136 (MQ=255)
                                                                                                          gATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTtct                                                                              >  6:63534/1‑110 (MQ=255)
                                                                                                          gATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTtct                                                                              <  5:63534/110‑1 (MQ=255)
                                                                                                          gATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCagcgagc                                       <  3:211808/149‑1 (MQ=255)
                                                                                                                     ggAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCgg                            <  4:50064/149‑1 (MQ=255)
                                                                                                                     ggAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCgg                            <  4:138906/149‑1 (MQ=255)
                                                                                                                      gAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCgg                            >  7:323512/1‑148 (MQ=255)
                                                                                                                      gAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCgg                            >  3:264771/1‑148 (MQ=255)
                                                                                                                           acg‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAAc                      >  5:385722/1‑149 (MQ=255)
                                                                                                                           acg‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAAc                      <  1:61448/149‑1 (MQ=255)
                                                                                                                             g‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACgg                    <  3:178934/149‑1 (MQ=255)
                                                                                                                                    gTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTg                                   <  3:20114/129‑1 (MQ=255)
                                                                                                                                    gTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTg                                   >  4:20114/1‑129 (MQ=255)
                                                                                                                                     tCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCtt              <  3:281983/149‑1 (MQ=255)
                                                                                                                                     tCTGGCGTGGCCACCGGTAGCATGGAATGGCTGAGGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCtt              <  8:189737/149‑1 (MQ=255)
                                                                                                                                         gCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCg                             <  5:225395/130‑1 (MQ=255)
                                                                                                                                         gCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCg                             >  6:225395/1‑130 (MQ=255)
                                                                                                                                         gCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTgtcg          <  2:188670/149‑1 (MQ=255)
                                                                                                                                           gTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTgtcgtc        >  7:70009/1‑149 (MQ=255)
                                                                                                                                              gCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGtt                                                                                                  <  4:192923/56‑1 (MQ=255)
                                                                                                                                              gCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGtt                                                                                                  >  3:192923/1‑56 (MQ=255)
                                                                                                                                                 aCCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCAcc  >  5:168744/1‑149 (MQ=255)
                                                                                                                                                 aCCGGTAGCATGGAATGGCTGATGAATCCGAAACCAGGGTTGTAGGCGTCGTTGAGCGTGCCGAACGCATTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGCCCAGCGAGCGTTGCGAGCGGGAGCACGGCAGCTTGTCGTAGTcccc  >  1:215504/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCCGCGGTTGAGCGCGGCGATGGCGTTGGCGCTGAGGGCGCCGAAACTCATCGCCGAGATGTTGAAGATCGAGGCCGAGTAGGGCATTCGGCATTGCGGCCCGCCGATGGCAATGCGGAACGACG‑‑CAGGGTCTGGCGTGGCCACCGGTAGCATGGAATGGCTGATGAATTCGAAACCAGGTTTGTAGGCGTCGTTGAGCGTGCCGAACGCTTTCTCGGCGCTTTCGTTCTTGGCGCGGGCGTAGACCAGCGAGCGTTGCGAGCGGGAGAACGGCAGCTTGTCGTCGTCACC  >  NC_002947/327189‑327480

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