Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,594,273 C→T 5.9% A529A (GCG→GCA tmcA ← elongator methionine tRNA (ac4C34) acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,594,2730CT 5.9% 217.8 / 2.2 69A529A (GCG→GCAtmcAelongator methionine tRNA (ac4C34) acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CGTCGTCAGAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAGGTCTTGCGTATATTGCAAAGC  >  NC_000913/2594131‑2594422
                                                                                                                                              |                                                                                                                                                     
ggtcgtcAGAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCg                                                                                                                                                <  1:545272/149‑1 (MQ=255)
 gtcgtcAGAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGTGCCATCGc                                                                                                                                               >  1:644454/1‑150 (MQ=255)
 gtcgtcAGAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGTGCCATCGc                                                                                                                                               >  1:594314/1‑150 (MQ=255)
  tcgtcAGAAAGGACGGCATCGTTTAGTGGAGCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGcc                                                                                                                                              <  1:263511/150‑1 (MQ=255)
   cgtcAGAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCg                                                                                                                                             >  1:149160/1‑150 (MQ=255)
        gAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGCCGTAACCGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAg                                                                                                                                        >  2:104271/1‑150 (MQ=255)
         aaaGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGc                                                                                                                                       <  1:1151171/150‑1 (MQ=255)
           aGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCaa                                                                                                                                     >  2:1331142/1‑150 (MQ=255)
             gACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAAcc                                                                                                                                   <  2:889176/150‑1 (MQ=255)
                ggCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTCCCCGCATCACTCATCGGTAACAGCGCCCTCGCCGTATAGCAACCgct                                                                                                                                >  2:99135/1‑150 (MQ=255)
                 gCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCgctg                                                                                                                               >  1:123597/1‑150 (MQ=255)
                 gCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCgctg                                                                                                                               >  1:843363/1‑150 (MQ=255)
                      gTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGc                                                                                                                          <  2:888898/150‑1 (MQ=255)
                      gTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGc                                                                                                                          <  2:1020398/150‑1 (MQ=255)
                       tttAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCt                                                                                                                         >  2:972753/1‑150 (MQ=255)
                          aGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCAGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTcc                                                                                                                      >  1:664016/1‑150 (MQ=255)
                           gTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTccc                                                                                                                     >  2:514519/1‑150 (MQ=255)
                            tGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGGCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCg                                                                                                                    <  2:234625/150‑1 (MQ=255)
                             ggATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCga                                                                                                                   >  1:1281842/1‑150 (MQ=255)
                                  aaCAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGGAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGAtt                                                                                                              <  1:836066/150‑1 (MQ=255)
                                   aCAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTa                                                                                                             <  2:1083644/150‑1 (MQ=255)
                                   aCAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTa                                                                                                             <  2:1089020/150‑1 (MQ=255)
                                           cGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCg                                                                                                     <  1:523232/150‑1 (MQ=255)
                                             tttCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCa                                                                                                   <  2:115588/150‑1 (MQ=255)
                                                cGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCa                                                                                                >  2:1213529/1‑150 (MQ=255)
                                                          cTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAccc                                                                                      >  1:150103/1‑150 (MQ=255)
                                                            gcgcGAGAGCTTGCGCATCGCGACGTAAACGGTACTGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGc                                                                                    <  2:1288491/150‑1 (MQ=255)
                                                            gcgcGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGCTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGc                                                                                    <  2:58032/150‑1 (MQ=255)
                                                              gcgAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAg                                                                                  >  1:1301481/1‑150 (MQ=255)
                                                               cgAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGc                                                                                 >  2:1170597/1‑150 (MQ=255)
                                                                 agagCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGt                                                                               >  1:57478/1‑150 (MQ=255)
                                                                    gCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGc                                                                            >  2:435040/1‑150 (MQ=255)
                                                                    gCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGc                                                                            >  2:343/1‑150 (MQ=255)
                                                                    gCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGACAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCAGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGc                                                                            >  1:794407/1‑150 (MQ=255)
                                                                      ttGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCa                                                                          <  2:772709/150‑1 (MQ=255)
                                                                       tGCGCATCGCGACGTAAACTGTAATTCTCACGTTCAGCCAGCTGTTTACTCGCATCACTCATCGGTAACAGGGCCATCGACGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAg                                                                         <  2:529923/150‑1 (MQ=255)
                                                                         cgcATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGaa                                                                       >  2:493150/1‑150 (MQ=255)
                                                                            aTCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGTGCCATCGCCGTATAGCTACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACg                                                                    <  2:644454/150‑1 (MQ=255)
                                                                            aTCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGTGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACg                                                                    <  2:594314/150‑1 (MQ=255)
                                                                             tCGCGACGTAAACTGTAATGCTCACGTTAAGCCAGCTGTTTACCCGCATCACTCATCTGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTTCCATAAACGc                                                                   <  2:800795/150‑1 (MQ=255)
                                                                             tCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGc                                                                   >  1:771416/1‑150 (MQ=255)
                                                                             tCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGc                                                                   >  1:776660/1‑150 (MQ=255)
                                                                                cgACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTTCCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCAt                                                                <  2:792030/150‑1 (MQ=255)
                                                                                cgACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCAt                                                                <  2:500144/150‑1 (MQ=255)
                                                                                       aaCGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTccc                                                         <  2:1067798/150‑1 (MQ=255)
                                                                                         cGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCCCCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATACCTCCCCg                                                       >  2:1220697/1‑150 (MQ=255)
                                                                                          ggTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCGGCACAACACCGCAGCGTTGCCAGAAACGCCATACCTCCCCgg                                                      >  1:240666/1‑150 (MQ=255)
                                                                                          ggTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCgg                                                      <  2:776660/150‑1 (MQ=255)
                                                                                          ggTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCgg                                                      <  1:972753/150‑1 (MQ=255)
                                                                                          ggTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCgg                                                      <  2:771416/150‑1 (MQ=255)
                                                                                            tAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGgtg                                                    <  1:937571/150‑1 (MQ=255)
                                                                                            tAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGGTGCCAGAAACGCCATAACTCCCCGgtg                                                    <  2:732536/150‑1 (MQ=255)
                                                                                             aaTGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGgtgt                                                   <  1:978325/150‑1 (MQ=255)
                                                                                              aTGCTCACGTTCAGTCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTa                                                  >  2:547342/1‑150 (MQ=255)
                                                                                               tGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTaa                                                 <  1:99135/150‑1 (MQ=255)
                                                                                                gCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAAc                                                <  2:1186483/150‑1 (MQ=255)
                                                                                                  tCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCa                                              <  2:1043782/150‑1 (MQ=255)
                                                                                                    aCGTTCAGCCAGCTGTTTACCCGCCTCACACCTCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCGGGCTTCCCGATGATTACCCATCCGCACCCGGCCAAACCCGCCGCGTTGCCAGGAACGCCAGTACTCCCCGGTGTAACCaaa                                            >  2:239452/1‑150 (MQ=255)
                                                                                                      gTTCAGCCAGCGGTTTACCCGGATCACTCATCGGTATCAGCGCCATCGCCGTTTAGCAACCGCTGCTGGATTCCAGTTGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAAc                                          <  2:1186121/150‑1 (MQ=255)
                                                                                                          agccagcTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTACCCAAAACTCAc                                      >  1:889837/1‑150 (MQ=255)
                                                                                                          agccagcTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCAc                                      >  2:155886/1‑150 (MQ=255)
                                                                                                                   tttACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATa                             <  1:1213529/150‑1 (MQ=255)
                                                                                                                    ttACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAg                            <  2:843363/150‑1 (MQ=255)
                                                                                                                       cccGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGATAGATAGTCg                         <  2:1012991/150‑1 (MQ=255)
                                                                                                                        ccGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCCCCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCCCCGAAAGATAGTCGa                        >  1:13598/1‑150 (MQ=255)
                                                                                                                         cGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCCCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAg                       <  2:123597/150‑1 (MQ=255)
                                                                                                                                   aTCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAGGTCTTGCGtt             >  2:1018955/1‑149 (MQ=255)
                                                                                                                                      ggTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAGGTCTTGCGTatat          >  2:661894/1‑150 (MQ=255)
                                                                                                                                       gTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAGGTCTTGCGTATAtt         <  1:1296691/150‑1 (MQ=255)
                                                                                                                                              cgcCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAGGTCTTGCGTATATTGCAAAGc  <  1:1331142/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CGTCGTCAGAAAGGACGGCATCGTTTAGTGGATCAACAGGAAGCGTTTCGCCATTCCACTGCGCGAGAGCTTGCGCATCGCGACGTAAACGGTAATGCTCACGTTCAGCCAGCTGTTTACCCGCATCACTCATCGGTAACAGCGCCATCGCCGTATAGCAACCGCTGCTGGCTTCCCGATGATTACCCATCCGCACCAGCACAAAACCGCAGCGTTGCCAGAAACGCCATAACTCCCCGGTGTAACCAAAACTCACCGAAAGATAGTCGAGGTCTTGCGTATATTGCAAAGC  >  NC_000913/2594131‑2594422

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