Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,942,042 G→T 5.9% T175N (ACT→AAT)  gcvA ← glycine cleavage system transcriptional activator; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,942,0420GT 5.9% 235.7 / 3.3 85T175N (ACT→AAT) gcvAglycine cleavage system transcriptional activator; autorepressor
Reads supporting (aligned to +/- strand):  ref base G (33/47);  new base T (2/3);  total (35/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GGCGCTATGGCTAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGATCAACCGCCTGGAT  >  NC_000913/2941894‑2942187
                                                                                                                                                    |                                                                                                                                                 
ggCGCTATGGCTAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCCCGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAg                                                                                                                                                   <  2:339706/149‑1 (MQ=255)
  cgtagtGGCTAAATATTGGCCCTTGCTGAACGTTGATATGATTGAAGTCGTATTGTTTGGTATATGTAAGCCAGTCACGCCGCGAGTCATCATGTAGTAACGTAGGCTTTGCGTTGGCGACTGGTGTCTTCAAGGGTGTTTCGCGAGTc                                                                                                                                                 <  2:353779/144‑1 (MQ=255)
   gcTATGGCTAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTca                                                                                                                                                >  1:811467/1‑149 (MQ=255)
    cTATGGCTAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGCCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCAGTGTCTTCAAGGGTTTTTCGCCAGTcag                                                                                                                                               <  1:761357/149‑1 (MQ=255)
       tGGCTAAAAATTGGCCCTTGCTGAACGTTGATATGGTTTAACCCCAACTGTCGGGTATATGTCTGCCAGTAACGGCGCGCAGCATCATGTAATAAGGTCTGTTTAGCCAGATCTTCCGGTGCCTTCAAGGGTTTTTCGACTGTcagcag                                                                                                                                            >  2:666568/1‑149 (MQ=255)
         gCTAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTa                                                                                                                                          >  1:768927/1‑149 (MQ=255)
         gCTAAAAATTGGCCCCTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTCTGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTa                                                                                                                                          <  2:1018761/149‑1 (MQ=255)
           tagaaaTTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTCgc                                                                                                                                        >  2:1246459/4‑149 (MQ=255)
           tAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCTGGGATATGTCTGCCAGTCACGGCGCGAAGCATCCAGCAATTTCCTACGGGTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTAAGCAGTAgg                                                                                                                                        >  1:314929/1‑148 (MQ=255)
              aaaTTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCATGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAgcggc                                                                                                                                     <  1:437131/149‑1 (MQ=255)
                atcggCCCTTGCTGAACGTTGATATGCTTTAACCCCAACTGTCGGCTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGa                                                                                                                                   <  2:950755/146‑1 (MQ=255)
                aTTGGCCCGTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGa                                                                                                                                   <  2:45780/149‑1 (MQ=255)
                 ttGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGaa                                                                                                                                  <  2:86336/149‑1 (MQ=255)
                 tgggCCCTTGATGAACGTTGACATGATTTAACACCAACTGGCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGaa                                                                                                                                  <  1:1070332/147‑1 (MQ=255)
                 gtGGCCCTTGCTGAACGTTGCTATGATTTAACCCCACCTGTCGGGTATATGTCTGCCCGTCACGGCGCGAAGCATCATGTAATAACGTATGTTCAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGac                                                                                                                                  <  1:1057486/148‑2 (MQ=255)
                     cccTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTGGCGGCGAacaca                                                                                                                              >  1:1127374/1‑149 (MQ=255)
                     cccTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAacaca                                                                                                                              >  2:1329146/1‑149 (MQ=255)
                     cccTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAACCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGCCGAacaca                                                                                                                              >  2:332864/1‑149 (MQ=255)
                      ccTTGCTGGACGTTGATATGATTTAACCCCACCTGCCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAacacac                                                                                                                             <  2:1171803/149‑1 (MQ=255)
                      ccTTGCTGAACGTTGATATGATTTAACCCCCACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCATAGCGCAGTAGCGGCGAacacac                                                                                                                             <  2:387510/149‑1 (MQ=255)
                      ccTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCCGTCAGCAGTAGCGGCGAacacac                                                                                                                             >  2:1127247/1‑149 (MQ=255)
                      ccTTGCTGAACGTTCATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCATAGCGCAGTAGCGGCGAacacac                                                                                                                             <  2:387684/149‑1 (MQ=255)
                      ccTTGCTGAACGGTGATATGATTTAACCCCAACTGCCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTGCCGGTGTCTTCAAGGGTTTTTCGCCATAGCGCAGTAGCGGCGAacacac                                                                                                                             <  2:517755/149‑1 (MQ=255)
                             gAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTCGCCCGCGTTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAAAACACAGGCaat                                                                                                                      >  1:649602/1‑149 (MQ=255)
                                     atatGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTGGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTc                                                                                                              >  1:295743/1‑149 (MQ=255)
                                        tGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGc                                                                                                           >  1:736359/1‑149 (MQ=255)
                                         gATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCg                                                                                                          >  2:357115/1‑149 (MQ=255)
                                             tAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACa                                                                                                      >  2:1269526/1‑149 (MQ=255)
                                             tAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACa                                                                                                      <  1:420166/149‑1 (MQ=255)
                                                 cccAGCGGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGAGAATTCAAGGTGTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGttt                                                                                                  <  1:260141/149‑1 (MQ=255)
                                                   cgtcTGTCTGGTATCTGTCTGCCAGTCACGTCTCGAAGCATCATGGAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTGAATGTTTTGTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGttttt                                                                                                <  1:945787/146‑1 (MQ=255)
                                                    accTGCCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTc                                                                                               <  1:407074/147‑1 (MQ=255)
                                                      cTGTCGGGTATTTGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTTAAGGGTTTTTCGTCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCa                                                                                             <  2:3759/149‑1 (MQ=255)
                                                       cgccGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCAc                                                                                            <  1:1106708/146‑1 (MQ=255)
                                                          ctggTATCTGTCCGCCAGCCTCGGCCCGAAGCATCATGTAATAACGTATCTGTCGCCAGAGCTTCCGGTGTCTTCAAGGTTTTTTCGCCAGTCAGCAGTATCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCg                                                                                         <  2:48342/147‑1 (MQ=255)
                                                          cGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCg                                                                                         <  2:1063682/149‑1 (MQ=255)
                                                           gggTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGGCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCCACAAATATTCGGCGTACAGTTTTTCCACACGt                                                                                        >  2:90144/1‑149 (MQ=255)
                                                             gTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAg                                                                                      >  1:1149948/1‑149 (MQ=255)
                                                                  tGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCCGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCg                                                                                 <  2:224948/149‑1 (MQ=255)
                                                                   gTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCgg                                                                                >  1:276827/1‑149 (MQ=255)
                                                                   gTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCgg                                                                                >  2:143185/1‑149 (MQ=255)
                                                                   gTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAACCCCgg                                                                                >  1:398840/1‑149 (MQ=255)
                                                                        ccAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTGCCACCCGTAGCCCCGGCCAGt                                                                           <  2:768927/149‑1 (MQ=255)
                                                                          aGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCGCCCGTATCCCCGGCCAGTTg                                                                         >  1:146807/1‑149 (MQ=255)
                                                                          aGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTg                                                                         <  1:781764/149‑1 (MQ=255)
                                                                           gTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGc                                                                        >  2:509221/1‑149 (MQ=255)
                                                                            tCACGGCGCGAAGCAGCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGcc                                                                       <  1:779169/149‑1 (MQ=255)
                                                                              aCGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGAAGTCTACAAGGGCCTCTTGTCAGTCAGCAGCAACGGCGCACACCGCGTCAATAAATATTCCGCGTACATTGTTTCCACCCGTAGCCCCGGCCAGTTGcccc                                                                     >  2:1332696/1‑149 (MQ=255)
                                                                                 gcgcGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCTCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGAc                                                                  >  1:716041/1‑149 (MQ=255)
                                                                                 gcgcGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCACTAAATATTCGGCGTACAGCTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGAc                                                                  <  2:525603/149‑1 (MQ=255)
                                                                                  cgcgAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTCCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGCCGTACAGTTTTTCCACCCGACGCCCCGGCCCGTTGCACCGagc                                                                 >  2:964118/1‑147 (MQ=255)
                                                                                     gAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATa                                                              <  1:816039/149‑1 (MQ=255)
                                                                                     gAAACATCATGTAATCACGGATCTTTAGCCAGATCTGCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATa                                                              <  2:314929/149‑1 (MQ=255)
                                                                                         catcatGTAAGAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCGTGACCATAAAat                                                          <  2:1169528/149‑1 (MQ=255)
                                                                                            catTTAATAACGTATGTTTACCCAGATCTTCCGGTGTCTTCAAGCGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCACTAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATc                                                       <  1:621596/149‑1 (MQ=255)
                                                                                            catGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATc                                                       <  1:1127247/149‑1 (MQ=255)
                                                                                            catGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATc                                                       <  1:451272/149‑1 (MQ=255)
                                                                                            catGTAATAACGTATGTTGAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATc                                                       <  2:202738/149‑1 (MQ=255)
                                                                                              tGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGc                                                     <  2:828459/149‑1 (MQ=255)
                                                                                              tGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGc                                                     <  1:145327/149‑1 (MQ=255)
                                                                                              ggcaagaaCGTCTGTTTCGCCAGATCTTCCGGTGTCTTAAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGc                                                     <  2:102852/143‑1 (MQ=255)
                                                                                               gTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGcc                                                    >  1:1195279/1‑149 (MQ=255)
                                                                                                 aataaCGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCcac                                                  <  1:176330/149‑1 (MQ=255)
                                                                                                       gTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTACTTGCCATGCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCacatcaac                                            >  1:1336930/1‑149 (MQ=255)
                                                                                                        tATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTTCCCCGACCCTAAAATATCGCCgcatcaacc                                           >  2:83472/1‑148 (MQ=255)
                                                                                                        tATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTACCCCGCGGAGCCCCGGCCAGTTGCCCCGACCATACAATATCGCCacagcagct                                           >  1:518128/1‑146 (MQ=255)
                                                                                                              tAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATcc                                     <  1:663304/149‑1 (MQ=255)
                                                                                                                cccAGATCTTCCGGCGTCCTCACGGGTTTTTTGCCAGTCAGCAGTAGCGGCGAACCCACCGGCAATAAATATTCCGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGc                                   <  2:449561/148‑1 (MQ=255)
                                                                                                                 ccAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATCTCGCCACATCAACATCATCCGcc                                  >  2:1055607/1‑149 (MQ=255)
                                                                                                                  cAGGTCTTCCGGTGCCCTCAAGGTCGTTTCGCCAGTCAGCAGTAGCGGCGTACATACCTCCAATAAATATACGGCGTACAGTTTTTCCACCCGTATCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCa                                 <  2:39400/149‑1 (MQ=255)
                                                                                                                   aGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACGGTTTATCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAg                                <  1:593402/149‑1 (MQ=255)
                                                                                                                   aGATCTTCCGGTGTCTTAAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGTCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAg                                <  2:649602/149‑1 (MQ=255)
                                                                                                                   aGATCTCCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCCACATCATCCGCCAg                                >  1:748154/1‑149 (MQ=255)
                                                                                                                    gATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTCCAGTTTTTCCAACCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCACCATCATCCGCCAga                               >  2:472898/1‑148 (MQ=255)
                                                                                                                        ttCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTAt                           <  1:632157/149‑1 (MQ=255)
                                                                                                                         tCCGGTGTCGTCAAGGGGTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATGCGGCGTACAGTTTTGCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATc                          <  1:709607/149‑1 (MQ=255)
                                                                                                                           cGGTGTCTTCAAGGGTTTTTTTTTTTTTTTCAGTAGCGGCGAACACACCGGCAATAAATATTGGGCGTACAGTTTTTCCACCCGTAGCCCCTGCCACTTGCCCCGACCATAAATTTTCGCCCCATTAACATCATCCGCCCGCTTCTCtt                        >  2:4422/1‑149 (MQ=255)
                                                                                                                             gtgtCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTcc                      <  1:808825/149‑1 (MQ=255)
                                                                                                                             gtgtCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTcc                      <  1:897908/149‑1 (MQ=255)
                                                                                                                                 cTTCAAGGGTTTTTGGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGAc                  <  1:863940/149‑1 (MQ=255)
                                                                                                                                    cAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAACTATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGAt               <  1:1066936/149‑1 (MQ=255)
                                                                                                                                      aGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGATCa             >  1:18951/1‑149 (MQ=255)
                                                                                                                                       gggTTTTACGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGATCaa            <  2:805049/149‑1 (MQ=255)
                                                                                                                                          tttttCGCCAGTCAGCAGTAGCGGCGAACATACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGATCAACCg         >  1:708516/1‑149 (MQ=255)
                                                                                                                                             ttCGCCAGCCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGATCAACCGCCt      <  2:1033853/149‑1 (MQ=255)
                                                                                                                                                 ccAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCCGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCGCCTTCTCTTCCTGACGATCAACCGCATGGAt  >  2:1289132/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCGCTATGGCTAAAAATTGGCCCTTGCTGAACGTTGATATGATTTAACCCCAACTGTCGGGTATATGTCTGCCAGTCACGGCGCGAAGCATCATGTAATAACGTATGTTTAGCCAGATCTTCCGGTGTCTTCAAGGGTTTTTCGCCAGTCAGCAGTAGCGGCGAACACACCGGCAATAAATATTCGGCGTACAGTTTTTCCACCCGTAGCCCCGGCCAGTTGCCCCGACCATAAAATATCGCCACATCAACATCATCCGCCAGCTTATCTTCCTGACGATCAACCGCCTGGAT  >  NC_000913/2941894‑2942187

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: