Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,210,435 Δ1 bp 6.0% coding (109/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,4350C. 6.0% 276.5 / 10.4 83coding (109/1014 nt)tsaDN(6)‑L‑threonylcarbamoyladenine synthase, TsaD subunit
Reads supporting (aligned to +/- strand):  ref base C (34/43);  new base . (2/3);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

CCGACTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACC  >  NC_000913/3210287‑3210584
                                                                                                                                                    |                                                                                                                                                      
ccGACTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCg                                                                                                                                                       <  1:595092/150‑1 (MQ=255)
 cGACTAATCCAGGGCCGGCGGCATAGGCCACAGCATCAATGTCTTTTGCCGTAAAACCCGACTCCTTTAGCGCCGCCGGGCTCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGGTCAGGCACGACGCCGCCGTAGTCAGCGt                                                                                                                                                      <  1:1093900/150‑1 (MQ=255)
  gACTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTg                                                                                                                                                     >  1:1261538/1‑150 (MQ=255)
      aaTCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAAt                                                                                                                                                 >  2:1115754/1‑150 (MQ=255)
          cAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCa                                                                                                                                             <  1:1296202/150‑1 (MQ=255)
          cAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCa                                                                                                                                             >  1:119569/1‑150 (MQ=255)
           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:920015/150‑1 (MQ=255)
                                             tttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCtttt                                                                                                         <  1:1239485/150‑1 (MQ=255)
                                              ttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCttttt                                                                                                        <  1:90859/150‑1 (MQ=255)
                                              ttGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCttttt                                                                                                        <  1:981240/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:797125/1‑150 (MQ=255)
                                                      aaaCCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCAtcgtcg                                                                                                <  2:236213/150‑1 (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                      <  2:1387370/150‑1 (MQ=255)
                                                                                                                                cACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTtta                      >  1:1412474/1‑150 (MQ=255)
                                                                                                                                  cgacgCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatt                    <  1:1410618/150‑1 (MQ=255)
                                                                                                                                  cgacgCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTttatt                    <  1:641485/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)
                                                                                                                                                    |                                                                                                                                                      
CCGACTAATCCAGGGCCTGCGGTATAGGCCACAGCATCAATGTCTTTTGCCGTTAAACCAGACTCCTTTAGCGCCGCCTGGATCAACGGTACGGTTTTACGCACATGATCGCGGGAGGCCAGTTCAGGCACGACGCCGCCGTAGTCAGCGTGCAATTTCACCTGACTATACAATTGGTTGGCTAACAAACCTTTTTCATCGTCGTAAATGGCGATGCCGGTTTCATCGCAGGAAGTTTCAATACCCAGTACACGCATGACTTGTTTTACCTCGCTTTATTACCGCGCAGTGTAGGACC  >  NC_000913/3210287‑3210584

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: