Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,447:1 +G 6.5% coding (97/1014 nt) tsaD ← N(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,210,4471.G 6.5% 307.9 / 17.4 77coding (97/1014 nt)tsaDN(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit
Reads supporting (aligned to +/- strand):  ref base . (31/41);  new base G (2/3);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

AGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGTTG  >  NC_000913/3210298‑3210595
                                                                                                                                                      |                                                                                                                                                    
aGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                       <  1:1254754/150‑1 (MQ=255)
     cTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGAc                                                                                                                                                 <  1:1039554/150‑1 (MQ=255)
        cGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACtac                                                                                                                                              <  2:470219/150‑2 (MQ=255)
          gTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACtacgg                                                                                                                                            <  1:1223672/150‑4 (MQ=255)
           tataGGCCACTGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACa                                                                                                                                           <  1:848462/150‑1 (MQ=255)
           tataGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGGTGAAATTGCACGCTGACtacggc                                                                                                                                           <  2:432476/150‑5 (MQ=255)
            ataGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACaa                                                                                                                                          >  2:1046426/1‑150 (MQ=255)
               ggCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAAttg                                                                                                                                       <  1:556914/150‑1 (MQ=255)
                           aaTTTCTTTTGCCATTAAAACAAACTCCTTTAACGCAACCTGGATCAACGCTACTTTTTTACGCACATAATCTCGTGAAGACAGTTCAGGAACGACGCCGCCGTAATCAACTTGCAATTTCACCTGACTATACAATTGGTTGGCTAACaa                                                                                                                           >  2:609537/1‑150 (MQ=255)
                           aaTGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACaa                                                                                                                           >  1:1387370/1‑150 (MQ=255)
                             tGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                       <  2:142381/121‑1 (MQ=255)
                             tGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCAc                                                                                                                                                       >  1:142381/1‑121 (MQ=255)
                                cTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCtt                                                                                                                      >  2:1337441/1‑150 (MQ=255)
                                cTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCtt                                                                                                                      >  2:373242/1‑150 (MQ=255)
                                  tttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCtttt                                                                                                                    <  1:1239485/150‑1 (MQ=255)
                                  tttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCtttt                                                                                                                    <  1:920015/150‑1 (MQ=255)
                                   ttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCttttt                                                                                                                   <  1:981240/150‑1 (MQ=255)
                                   ttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCttttt                                                                                                                   <  1:90859/150‑1 (MQ=255)
                                     gCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCa                                                                                                                 >  1:458839/1‑150 (MQ=255)
                                      ccGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCTGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAt                                                                                                                <  1:1297254/150‑1 (MQ=255)
                                      ccGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAt                                                                                                                >  1:506854/1‑150 (MQ=255)
                                       cGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtc                                                                                                               >  2:143078/1‑150 (MQ=255)
                                       cGTTAAACAAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtc                                                                                                               >  2:745726/1‑150 (MQ=255)
                                           aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                           <  2:236213/150‑1 (MQ=255)
                                           aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                           >  2:797125/1‑150 (MQ=255)
                                           aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACACCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                           <  1:923454/150‑1 (MQ=255)
                                             aCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTa                                                                                                         <  1:1316993/150‑1 (MQ=255)
                                              ccAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTaa                                                                                                        >  2:913864/1‑150 (MQ=255)
                                                aGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAAt                                                                                                      <  1:33368/150‑1 (MQ=255)
                                                   cTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGc                                                                                                   >  1:290650/1‑150 (MQ=255)
                                                      cTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGAt                                                                                                >  2:1418802/1‑150 (MQ=255)
                                                          aGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCg                                                                                            >  2:928076/1‑150 (MQ=255)
                                                           gcgcCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCgg                                                                                           >  2:1410618/1‑150 (MQ=255)
                                                            cgccgcCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGt                                                                                          >  2:1115555/1‑150 (MQ=255)
                                                             gccgccTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGtt                                                                                         <  1:1419813/150‑1 (MQ=255)
                                                               cgccTGGATCACCGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTGCAGGCACGACGCCGCCGTAGGCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTc                                                                                       <  2:576852/150‑1 (MQ=255)
                                                               cgccTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTc                                                                                       <  1:1115555/150‑1 (MQ=255)
                                                                gccTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCa                                                                                      >  2:908929/1‑150 (MQ=255)
                                                                  cTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATc                                                                                    >  1:588532/1‑150 (MQ=255)
                                                                    ggATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGc                                                                                  <  2:1070111/150‑1 (MQ=255)
                                                                     gATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCa                                                                                 >  2:876105/1‑150 (MQ=255)
                                                                     gATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCa                                                                                 >  1:808935/1‑150 (MQ=255)
                                                                               acggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCa                                                                       <  1:1321703/150‑1 (MQ=255)
                                                                               acggtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCa                                                                       <  1:1165065/150‑1 (MQ=255)
                                                                                  gtTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATa                                                                    >  1:1201738/1‑150 (MQ=255)
                                                                                    tttACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATAcc                                                                  <  1:797125/150‑1 (MQ=255)
                                                                                    tttACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATAcc                                                                  <  1:1379598/150‑1 (MQ=255)
                                                                                    tttACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATAcc                                                                  <  1:1368121/150‑1 (MQ=255)
                                                                                      tACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCa                                                                <  1:1046426/150‑1 (MQ=255)
                                                                                        cGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGt                                                              >  1:4917/1‑150 (MQ=255)
                                                                                        cGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACATTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGt                                                              >  2:489879/1‑150 (MQ=255)
                                                                                          cacaTGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTac                                                            <  2:119569/150‑1 (MQ=255)
                                                                                           acaTGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTaca                                                           <  1:537575/150‑1 (MQ=255)
                                                                                                aTCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCAt                                                      <  1:1115754/150‑1 (MQ=255)
                                                                                                  cgctGGATGACATTTCAGGAACGACTCAGCCTTAGTCAGCGTGAAATTTCACCTGACTAAACAATTTGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGa                                                    >  2:350913/1‑150 (MQ=255)
                                                                                                    cgGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACt                                                  <  1:373242/150‑1 (MQ=255)
                                                                                                      ggAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTg                                                <  1:724815/150‑1 (MQ=255)
                                                                                                         ggCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGc                                                        <  1:710938/139‑1 (MQ=255)
                                                                                                         ggCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGc                                                        >  2:710938/1‑139 (MQ=255)
                                                                                                         ggCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGttt                                             <  1:143078/150‑1 (MQ=255)
                                                                                                          gCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGtttt                                            <  2:1192585/150‑1 (MQ=255)
                                                                                                           ccAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTa                                           <  1:913864/150‑1 (MQ=255)
                                                                                                             aGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTAcc                                         <  1:1157441/150‑1 (MQ=255)
                                                                                                              gTTCAGGCACGACGCCGCCGTAGTCAGAGTGCAATTTCACCTGACTATACAATTGGTTGGATAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCt                                        >  2:505587/1‑150 (MQ=255)
                                                                                                                    gCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCACCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCttt                                  >  1:222642/1‑150 (MQ=255)
                                                                                                                     cACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTtta                                 >  1:1412474/1‑150 (MQ=255)
                                                                                                                     cACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTtta                                 <  2:1387370/150‑1 (MQ=255)
                                                                                                                       cgacgCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatt                               <  1:641485/150‑1 (MQ=255)
                                                                                                                       cgacgCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatt                               <  1:1410618/150‑1 (MQ=255)
                                                                                                                          cgccgcCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTATCAAACCTTTTTCATCGTCGTAAATGGTGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTAcc                            >  2:1235927/1‑150 (MQ=255)
                                                                                                                             cgccgTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGTCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTAc                                          <  1:613644/133‑1 (MQ=255)
                                                                                                                             cgccgTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTCGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGTCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTAc                                          >  2:613644/1‑133 (MQ=255)
                                                                                                                               gtatAGTCAGGTGAAATTGCACGCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCa                       >  2:806681/4‑150 (MQ=255)
                                                                                                                                 atAGTCAGGTGAAATTGCACGCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAgt                     >  2:5453/2‑150 (MQ=255)
                                                                                                                                  tAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAgtg                    <  1:629654/150‑1 (MQ=255)
                                                                                                                                    gTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTa                  <  1:1394544/150‑1 (MQ=255)
                                                                                                                                      cAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAgg                >  2:1091922/1‑150 (MQ=255)
                                                                                                                                         cGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGAcc             <  1:609537/150‑1 (MQ=255)
                                                                                                                                             aaaTTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATg         <  2:1261538/149‑1 (MQ=255)
                                                                                                                                                    aCCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGTTg  <  1:745726/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACCAATGCGGGTTG  >  NC_000913/3210298‑3210595

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: