Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,641 C→T 15.1% Q118Q (CAG→CAA tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,6410CT15.1% 1239.0 / 137.5 450Q118Q (CAG→CAAtktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (184/195);  new base T (31/39);  total (215/235)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 463 total aligned reads displayed.

AGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCAT  >  CP009273/3074493‑3074790
                                                                                                                                                    |                                                                                                                                                     
agCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCt                                                                                                                                                      <  1:408188/150‑1 (MQ=255)
   aagcTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGGGATGCCTTGAc                                                                                                                                                   >  2:24198/4‑150 (MQ=255)
   aaaCTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGAc                                                                                                                                                   >  1:293065/1‑150 (MQ=255)
         cGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCg                                                                                                                                             >  2:599500/1‑150 (MQ=255)
          gTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCgg                                                                                                                                            <  2:823963/150‑1 (MQ=255)
           tGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGa                                                                                                                                           >  2:944616/1‑150 (MQ=255)
              gagaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCg                                                                                                                                        <  1:1042025/150‑1 (MQ=255)
               agaTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCgg                                                                                                                                       >  2:816459/1‑150 (MQ=255)
                gaTGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGGGATGCCTTGACCTAACGGTCCggt                                                                                                                                      >  1:933516/1‑150 (MQ=255)
                 aTGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtg                                                                                                                                     >  1:336725/1‑150 (MQ=255)
                    ccTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCggtcgtc                                                                                                                                  <  1:777655/150‑2 (MQ=255)
                    ccTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCggtcgtc                                                                                                                                  >  1:264863/1‑149 (MQ=255)
                    ccTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCggtcgt                                                                                                                                   <  2:475970/149‑1 (MQ=255)
                     cTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtgg                                                                                                                                 <  2:193174/150‑1 (MQ=255)
                      ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtggt                                                                                                                                >  2:140935/1‑150 (MQ=255)
                      ttCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCggtggtggt                                                                                                                                >  2:487349/1‑150 (MQ=255)
                         cattaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTc                                                                                                                             <  2:450063/144‑1 (MQ=255)
                         cattaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTc                                                                                                                             >  2:154995/7‑150 (MQ=255)
                         catcatGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGttt                                                                                                                              <  1:569617/149‑1 (MQ=255)
                          attaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTcg                                                                                                                            <  1:355078/145‑2 (MQ=255)
                          attaagCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGCTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCCAAACCCACGGCAATGGTGATGCCTTGACCTAACCGTCCGGTCGTCGTTTc                                                                                                                             >  1:70629/6‑149 (MQ=255)
                            catGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCca                                                                                                                           >  2:423495/1‑149 (MQ=255)
                              agCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCgacg                                                                                                                         <  1:1004918/148‑2 (MQ=255)
                               gCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCacac                                                                                                                        >  1:553909/1‑149 (MQ=255)
                                 aGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAg                                                                                                                     >  1:881390/1‑150 (MQ=255)
                                  gCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGc                                                                                                                    <  1:239374/150‑1 (MQ=255)
                                  gCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccgg                                                                                                                    <  1:127615/150‑4 (MQ=255)
                                   ccGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCg                                                                                                                   <  2:53700/150‑1 (MQ=255)
                                   ccGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggc                                                                                                                   <  1:745289/150‑5 (MQ=255)
                                    cGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCGCGCGTTCGGCTATCGCGAAACCCACGGCATTGGTGATGCCTTGACCTAACGGTCCGGTCGTCGTTTCGACGcccggcg                                                                                                                  <  1:850370/150‑6 (MQ=255)
                                      tCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGgtg                                                                                                                <  2:1012426/150‑1 (MQ=255)
                                           cATGTAGTCGTAGTTGTATTTGTCGACTATGTCGTGGCCCTGTCTGTTTTTCTGCGCCGCCACCGTTTTTTCTGCAATCGCCATACCGTCTTCGTTTGCTATTCCCTGTCCCTTCGTCCCTTTTTTGGTTTCCACACCTTCGGTTTAAcc                                                                                                           >  2:950824/1‑150 (MQ=255)
                                            aTGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAAcc                                                                                                           >  2:591595/1‑149 (MQ=255)
                                             tGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCTTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCa                                                                                                         <  2:739181/150‑1 (MQ=255)
                                             tGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCa                                                                                                         >  2:280112/1‑150 (MQ=255)
                                                   cGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCg                                                                                                   >  1:195671/1‑150 (MQ=255)
                                                    gTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCgg                                                                                                  >  2:629869/1‑150 (MQ=255)
                                                      aGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                                                >  1:220076/1‑150 (MQ=255)
                                                       ggTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGt                                                                                                >  1:476584/1‑149 (MQ=255)
                                                        gtgtAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGa                                                                                              >  2:864373/1‑150 (MQ=255)
                                                           tAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                                                                            <  2:262391/149‑1 (MQ=255)
                                                           tAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGAcc                                                                                            <  2:643387/149‑1 (MQ=255)
                                                              tGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGa                                                                                        <  2:344278/150‑1 (MQ=255)
                                                                      aaTGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGCGTTTTa                                                                                  >  1:1013063/1‑148 (MQ=255)
                                                                      aaTGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTa                                                                                  >  2:1051860/1‑148 (MQ=255)
                                                                       aTGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTagag                                                                               <  1:12334/150‑1 (MQ=255)
                                                                         gTCGTGGCCAGGACGGTTAAACTGAGACGCCAGCGTTTTTTAAGCAATCGCCATACCGAAAGCGTTGGAAAAAAAAAGAAAAAGCGGACCGGTGGTGGATAACAAACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTg                                                                             <  1:684332/150‑1 (MQ=255)
                                                                          tCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGc                                                                            <  2:12624/150‑1 (MQ=255)
                                                                            gTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAg                                                                          <  1:136053/150‑1 (MQ=255)
                                                                               gCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTg                                                                       <  2:207768/150‑1 (MQ=255)
                                                                               gCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTg                                                                       <  2:376548/150‑1 (MQ=255)
                                                                                 ccGGACGGTTAACCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGAc                                                                     <  2:235526/150‑1 (MQ=255)
                                                                                 ccGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGAc                                                                     >  2:780582/1‑150 (MQ=255)
                                                                                   ggACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACg                                                                    >  2:83950/1‑149 (MQ=255)
                                                                                   ggACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACg                                                                    >  1:450983/1‑149 (MQ=255)
                                                                                    gACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACgg                                                                   >  2:892727/1‑149 (MQ=255)
                                                                                      cGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAg                                                                >  1:582687/1‑150 (MQ=255)
                                                                                       ggTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGt                                                               >  1:645560/1‑150 (MQ=255)
                                                                                        gTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtt                                                              <  2:1003366/150‑1 (MQ=255)
                                                                                        gTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtt                                                              >  1:496371/1‑150 (MQ=255)
                                                                                         ttAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGTAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttt                                                             >  1:829183/1‑150 (MQ=255)
                                                                                         ttAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCGATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttt                                                             >  2:565481/1‑150 (MQ=255)
                                                                                         ttAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttt                                                             <  1:180771/150‑1 (MQ=255)
                                                                                          tAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGtttt                                                            >  2:349286/1‑150 (MQ=255)
                                                                                           aaaCTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGttttt                                                           >  2:757571/1‑150 (MQ=255)
                                                                                               tGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGt                                                       >  1:279998/1‑150 (MQ=255)
                                                                                                 cgccgcCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttc                                                     <  1:506830/150‑1 (MQ=255)
                                                                                                 cgccgcCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttc                                                     >  1:984587/1‑150 (MQ=255)
                                                                                                 cgccgcCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGtt                                                      <  1:790636/149‑1 (MQ=255)
                                                                                                  gccgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttct                                                    <  2:1067366/150‑1 (MQ=255)
                                                                                                    cgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttcttc                                                  >  2:800748/1‑150 (MQ=255)
                                                                                                    cgccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGttcttc                                                  >  2:677458/1‑150 (MQ=255)
                                                                                                      ccAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCa                                                <  1:715040/150‑1 (MQ=255)
                                                                                                       cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                                               >  1:953954/1‑150 (MQ=255)
                                                                                                       cAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCAt                                               >  1:150354/1‑150 (MQ=255)
                                                                                                           gTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCTTCGGc                                           <  2:1024842/150‑1 (MQ=255)
                                                                                                           gTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCgg                                            >  2:50738/1‑149 (MQ=255)
                                                                                                           gTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGc                                           <  1:373900/150‑1 (MQ=255)
                                                                                                            atttttCTGCAATCGCCATACCGACTGCGTTGGCAATACACTGACCCAGCGGACCGATGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGCGACCCGGTGTTTTAGAGTGCAGCTGACGGAAGTATTTCAGTTCTTCCATCgtcg                                          >  2:300492/2‑147 (MQ=255)
                                                                                                               tttCTGCAATCGCCATACCGACTGCGTTGGCAATCCCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGTCCCGGAGTTTTAGAGGGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGa                                        <  1:1024107/149‑1 (MQ=255)
                                                                                                                 tCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATc                                      <  1:86766/149‑1 (MQ=255)
                                                                                                                  cTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCg                                     <  1:610856/149‑1 (MQ=255)
                                                                                                                  cTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCg                                     >  1:433079/1‑149 (MQ=255)
                                                                                                                   tGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGt                                    >  2:324782/1‑149 (MQ=255)
                                                                                                                     cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTaa                                  <  1:816459/149‑1 (MQ=255)
                                                                                                                     cAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTaa                                  <  2:727178/149‑1 (MQ=255)
                                                                                                                      aaTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAAcc                                <  1:731064/150‑1 (MQ=255)
                                                                                                                      aaTCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAAcc                                <  1:801648/150‑1 (MQ=255)
                                                                                                                         cGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGt                             >  2:708100/1‑150 (MQ=255)
                                                                                                                          gCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTg                            <  1:529995/150‑1 (MQ=255)
                                                                                                                                    cTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTgcag                    >  1:632804/1‑148 (MQ=255)
                                                                                                                                       cGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGcagcag                 <  2:697625/148‑1 (MQ=255)
                                                                                                                                        gTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCt              <  2:915285/150‑1 (MQ=255)
                                                                                                                                        gTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCt              <  1:319093/150‑1 (MQ=255)
                                                                                                                                               aTACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGAt        >  1:1059236/1‑149 (MQ=255)
                                                                                                                                                tACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGAt        >  2:657469/1‑148 (MQ=255)
                                                                                                                                                tACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGAt        >  1:668419/1‑148 (MQ=255)
                                                                                                                                                 aCCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAg     >  2:396355/1‑150 (MQ=255)
                                                                                                                                                   ccTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCa   <  2:620851/150‑1 (MQ=255)
                                                                                                                                                    cTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCAt  >  1:901574/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTTGGCAATACCCTGACCCAGCGGACCGGTGGTGGTTTCCACACCAGCGGTGTAACCCACTTCCGGGTGACCCGGAGTTTTAGAGTGCAGCTGACGGAAGTTTTTCAGTTCTTCCATCGGCAGATCGTAACCGGTGAGGTGCAGCAGGCTGTAGATCAGCAT  >  CP009273/3074493‑3074790

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