Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,088,892 G→C 5.4% P316R (CCG→CGG)  astA‑I ← arginine N‑succinyltransferase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,088,8920GC 5.4% 246.2 / 4.6 74P316R (CCG→CGG) astA‑Iarginine N‑succinyltransferase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base C (2/2);  total (37/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.27e-01

GTTGAGCGACTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTAC  >  NC_002947/5088744‑5089037
                                                                                                                                                     |                                                                                                                                                 
gTTGAGCGACTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACg                                                                                                                                                   <  6:244973/149‑1 (MQ=255)
gTTGAGCGACTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGACGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACg                                                                                                                                                   >  8:78195/1‑149 (MQ=255)
 ttGAGCGACTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACgg                                                                                                                                                  >  5:31522/1‑149 (MQ=255)
       gACTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGac                                                                                                                                            >  3:133367/1‑149 (MQ=255)
        aCTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCTGCAGGCGTTTGGCGGTCAACGGGTCGac                                                                                                                                            <  1:24946/149‑1 (MQ=255)
            cAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGaccac                                                                                                                                         <  8:28982/147‑1 (MQ=255)
                    cGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAgg                                                                                                                                     >  2:224624/1‑143 (MQ=255)
                    cGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAgg                                                                                                                                     <  1:224624/143‑1 (MQ=255)
                       ccTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCcggcg                                                                                                                            >  1:154435/1‑149 (MQ=255)
                        cTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCcggcgg                                                                                                                           <  5:138923/149‑1 (MQ=255)
                         tGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGa                                                                                                                                             >  3:168159/1‑130 (MQ=255)
                         tGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGa                                                                                                                                             <  4:168159/130‑1 (MQ=255)
                                  tGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTg                                                                                                                  >  4:181522/1‑148 (MQ=255)
                                      aaTTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGt                                                                                                             >  7:168762/1‑149 (MQ=255)
                                        ttGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGc                                                                                                           <  6:31522/149‑1 (MQ=255)
                                        ttGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGc                                                                                                           <  2:239264/149‑1 (MQ=255)
                                           ccTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCgg                                                                                                         >  2:175222/1‑148 (MQ=255)
                                           ccTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCgg                                                                                                         >  2:181318/1‑148 (MQ=255)
                                            cTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCa                                                                                                       <  3:36907/149‑1 (MQ=255)
                                             tGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAg                                                                                                      >  2:45243/1‑149 (MQ=255)
                                                gATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGa                                                                                                   >  2:22843/1‑149 (MQ=255)
                                                 aTGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGAt                                                                                                  <  1:45243/149‑1 (MQ=255)
                                                    tAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCTACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGc                                                                                                <  1:22843/148‑1 (MQ=255)
                                                      aTGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGc                                                                                                                 <  7:96844/129‑1 (MQ=255)
                                                      aTGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGc                                                                                                                 >  8:96844/1‑129 (MQ=255)
                                                      aTGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGa                                                                                              >  6:364479/1‑148 (MQ=255)
                                                      aTGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGa                                                                                              <  1:268186/148‑1 (MQ=255)
                                                                aTTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGac                                                                                                                                            >  7:382499/1‑92 (MQ=255)
                                                                aTTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGac                                                                                                                                            <  8:382499/92‑1 (MQ=255)
                                                                aTTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCa                                                                                                       <  7:352010/129‑1 (MQ=255)
                                                                aTTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCa                                                                                                       >  8:352010/1‑129 (MQ=255)
                                                                 tttAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCgcgc                                                                                  <  2:154435/149‑1 (MQ=255)
                                                                 tttAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCgcgc                                                                                  >  6:356745/1‑149 (MQ=255)
                                                                  ttAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTTGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCa                                                                                                       >  6:361350/1‑127 (MQ=255)
                                                                  ttAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTTGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCa                                                                                                       <  5:361350/127‑1 (MQ=255)
                                                                   tAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACa                                                                                            <  6:231852/137‑1 (MQ=255)
                                                                   tAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACa                                                                                            >  5:231852/1‑137 (MQ=255)
                                                                    aaTGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACCGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCgcgc                                                                                  <  8:32574/146‑1 (MQ=255)
                                                                    aaTGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACCGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCtag                                                                               >  7:32574/1‑147 (MQ=255)
                                                                     aTGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGc                                                                              >  8:116393/1‑149 (MQ=255)
                                                                       gggCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCgg                                                                            >  6:110198/1‑149 (MQ=255)
                                                                         gCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGtt                                                                          >  2:335999/1‑149 (MQ=255)
                                                                           cTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTtgt                                                                        <  4:133367/149‑1 (MQ=255)
                                                                            tCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCCACGGGTCGCCCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGACGTCGCGCTAGCGGTtgag                                                                       >  3:96296/1‑147 (MQ=255)
                                                                             cccGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGc                                                                                                >  7:61230/1‑123 (MQ=255)
                                                                             cccGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGc                                                                                                <  8:61230/123‑1 (MQ=255)
                                                                              ccGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGga                                                                     <  8:168762/149‑1 (MQ=255)
                                                                              ccGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGga                                                                     <  1:28911/149‑1 (MQ=255)
                                                                               cGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGgat                                                                    <  7:78195/149‑1 (MQ=255)
                                                                                 gCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTTCGGTTGTGgatga                                                                  <  5:356745/149‑1 (MQ=255)
                                                                                     gCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGaa                                                              <  7:116393/149‑1 (MQ=255)
                                                                                       cGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAg                                                            <  3:181522/149‑1 (MQ=255)
                                                                                          cAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTtgtg                                                                       <  1:355428/135‑1 (MQ=255)
                                                                                          cAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTtgtg                                                                       >  2:355428/1‑135 (MQ=255)
                                                                                            gCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTc                                                       <  1:181318/149‑1 (MQ=255)
                                                                                            gCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTc                                                       <  1:175222/149‑1 (MQ=255)
                                                                                                          aCATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTgg                                         >  2:374070/1‑149 (MQ=255)
                                                                                                               gTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGta                                     <  5:364479/148‑1 (MQ=255)
                                                                                                                   ccAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGtacctac                                >  5:86790/1‑149 (MQ=255)
                                                                                                                    cAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTg                                                                               >  1:304562/1‑101 (MQ=255)
                                                                                                                    cAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTg                                                                               <  2:304562/101‑1 (MQ=255)
                                                                                                                        gCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCa                           <  1:374070/149‑1 (MQ=255)
                                                                                                                         cccATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAg                          <  6:86790/149‑1 (MQ=255)
                                                                                                                          ccATGCGCAGGCGTTTGGTGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAg                          >  4:112880/1‑148 (MQ=255)
                                                                                                                            aTGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCAc                       >  6:19967/1‑149 (MQ=255)
                                                                                                                               cgcAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCa                     <  5:110198/148‑1 (MQ=255)
                                                                                                                                     gCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCt               >  4:218642/1‑149 (MQ=255)
                                                                                                                                     gCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCt               >  2:372019/1‑149 (MQ=255)
                                                                                                                                       gTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTgg             <  1:372019/149‑1 (MQ=255)
                                                                                                                                        cctggtggTCGACCCGTTGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGc                                             <  4:319896/114‑1 (MQ=255)
                                                                                                                                        cctggtggTCGACCCGTTGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGc                                             >  3:319896/3‑116 (MQ=255)
                                                                                                                                            gCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTc                                                  <  8:244537/107‑1 (MQ=255)
                                                                                                                                            gCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTc                                                  >  7:244537/1‑107 (MQ=255)
                                                                                                                                              ggTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTcgc      >  8:381989/1‑149 (MQ=255)
                                                                                                                                               gTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTcgcg     >  7:263494/1‑149 (MQ=255)
                                                                                                                                                tCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTcgcg     >  6:250928/1‑148 (MQ=255)
                                                                                                                                                  aaCGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTAc  >  6:424472/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTTGAGCGACTGCAGGGTTTCGCCCTGGCCGGCCTGCCAATTGCCTGCGATGTAATGCGTGGTCATTTAATGGGCCTCCCGGCTGGCCGACAGCGGCACCGCGCGGACATTGTCGCCAACGCCCATGCGCAGGCGTTTGGCGGTCAACGGGTCGACCACCAGGGTGCCGGCGGCCAACCGTGCAGGTGCGGCAGTGATGCGACAGTCGTCGCGCTTGCGGTTGTGGATGATGAAAGGGGTCGCGTCATCGCCTGGCGTACCTACCGCCAGCACCAGGGTCTGGCTCTCGCGTAC  >  NC_002947/5088744‑5089037

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: