Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 337,558 C→G 6.5% Y241* (TAC→TAG)  yahJ → putative metallo‑dependent hydrolase domain deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273337,5580CG 6.5% 239.5 / 7.2 76Y241* (TAC→TAG) yahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base C (34/37);  new base G (3/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

CTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGG  >  CP009273/337412‑337706
                                                                                                                                                  |                                                                                                                                                    
cTGAAAAATCTGCAAAATTTGCAGGCGGTACTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAATCGATGCAGGCGGGGGCGCATTACGt                                                                                                                                                    >  1:754644/1‑149 (MQ=255)
 tGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGAGGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGGACCTTTAATGCGTGCAGCGATGCAGGCGGGGGCGCATTACGTc                                                                                                                                                   >  2:320301/1‑149 (MQ=255)
 tGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAAACGTCGACATCCCGCAGCACGAGTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCCGGCGGGGGCGCATTACGTc                                                                                                                                                   >  1:732885/1‑149 (MQ=255)
  gAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTcg                                                                                                                                                  >  2:498108/1‑149 (MQ=255)
      aaTCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTAGGTcggcgg                                                                                                                                              <  2:79260/149‑1 (MQ=255)
      aaTCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTAGGTcggcgg                                                                                                                                              <  1:1235964/149‑1 (MQ=255)
         cTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCt                                                                                                                                           >  2:94407/1‑149 (MQ=255)
         cTGCAAAATTGGCAGGCGGTCCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGGCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCt                                                                                                                                           <  2:796024/149‑1 (MQ=255)
         cTGCAAAATCTGCCGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAATTTTGAACCTTTAATGCGTGAAGCGAGGCAGGCGGGGGCGCATTACGTCGGCGGCCt                                                                                                                                           <  2:1063112/149‑1 (MQ=255)
           gCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTgg                                                                                                                                         >  2:587801/1‑149 (MQ=255)
            cAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGa                                                                                                                                        >  1:953544/1‑149 (MQ=255)
            cAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCGTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGa                                                                                                                                        >  1:733580/1‑149 (MQ=255)
            cAAAATTCGCAGGCCGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGACCCTTTAATGCGTGAAGCGATGCAGGCGGGTGCGCATTACGTCGGCGGCCTGGa                                                                                                                                        <  2:642099/149‑1 (MQ=255)
             aaaaTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGACATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCCGCGGCGCGGAc                                                                                                                                       >  2:1163947/1‑149 (MQ=255)
             aaaaTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCGGGCGGGGGCGCATTACGTCGGCGGCCTGGAc                                                                                                                                       >  1:500504/1‑149 (MQ=255)
             aaaaTTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGAc                                                                                                                                       >  1:970626/1‑149 (MQ=255)
             aaaaTTGGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGCTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGAc                                                                                                                                       <  2:612045/149‑1 (MQ=255)
                       ggCGGTGCTGGCGCGGCGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAgtg                                                                                                                             <  1:159850/149‑1 (MQ=255)
                         cGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGtt                                                                                                                           <  1:1160340/149‑1 (MQ=255)
                               tGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTAGGTCGGCGGCCTGGACCCGGCCAGTGTTGATGgc                                                                                                                     >  2:936479/1‑149 (MQ=255)
                               tGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTAGGTCGGCGGCCTGGACCCGACCCGTGTTGATGgc                                                                                                                     >  2:835947/1‑149 (MQ=255)
                                   tggACGTCCGGCGGGCTTTTAGTGTGAAATCGGCGCCCTCCCGCAGCACGGTGTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGa                                                                                                                 <  1:534749/147‑1 (MQ=255)
                                   gcgACGTCCGGCGGGGTTTGCGTGTGGAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGa                                                                                                                 <  2:21018/149‑1 (MQ=255)
                                   gcgACGTCAGGCGGGTTTTGAGTGTGAAAGCGGCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGa                                                                                                                 <  1:579695/149‑1 (MQ=255)
                                    ggACGTCAGCCGGGCTTTGAGTGTGAAATCGTCGCATTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGAt                                                                                                                <  2:413837/148‑1 (MQ=255)
                                    cgACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGAt                                                                                                                <  2:413882/149‑1 (MQ=255)
                                    cgACGCCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGAt                                                                                                                <  1:360436/149‑1 (MQ=255)
                                       cGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCGTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGGGCATTACGTCGGCGGCCTGGACCCGACCACGGTTGATGCGGCGATGGa                                                                                                             >  1:1043787/1‑149 (MQ=255)
                                             gCGGGCTTTGAGTTTGAAATCTTCGCCGTCCCGCATCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCTATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATc                                                                                                       <  1:1118220/149‑1 (MQ=255)
                                              cGGGTTTGGCGTTTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGACATCTGAACCTTTAATGCGTGAAGCGAGGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGTCCCGACCAGTGTTGATGGCGCGATGGAAAAATcc                                                                                                      <  2:1070034/149‑1 (MQ=255)
                                               gggCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATccc                                                                                                     >  2:320708/1‑149 (MQ=255)
                                                 gCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTc                                                                                                   <  1:441339/149‑1 (MQ=255)
                                                   tttGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGa                                                                                                 >  2:646651/1‑149 (MQ=255)
                                                   tttGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGa                                                                                                 >  2:292959/1‑149 (MQ=255)
                                                    ttGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCCTTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCCCGATGGCAAAATCCCCCgcc                                                                                                >  2:850642/1‑147 (MQ=255)
                                                      gAGTGTGAAAGCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGacac                                                                                              <  1:703477/149‑1 (MQ=255)
                                                      gAGTGCGAAATCGTCGCCATCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCCTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGacac                                                                                              <  1:1062825/149‑1 (MQ=255)
                                                         tgtgAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGCTGGCGCGATGGAAAAATCCCTCGACACCAt                                                                                           >  1:588212/1‑149 (MQ=255)
                                                          gtgAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATg                                                                                          <  1:645622/149‑1 (MQ=255)
                                                           tgAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACCTCGGCGGCCTGGACCCGACCAGGGTTGATGGCGCGATGGAAAAATCCCGCGACACCATGt                                                                                         >  2:442506/1‑149 (MQ=255)
                                                             aaaTCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTc                                                                                       >  1:292373/1‑149 (MQ=255)
                                                                tcgtcgCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAg                                                                                    <  1:683936/149‑1 (MQ=255)
                                                                      ccTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTAGGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTgcg                                                                              >  1:1053217/1‑149 (MQ=255)
                                                                         ttccGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACCTCGGCGCCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTg                                                                           <  1:124699/147‑1 (MQ=255)
                                                                         tCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGGAAAATCCCTCGACACCATGTTCCAGAGTGCGGTg                                                                           >  2:427776/1‑149 (MQ=255)
                                                                           ccGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGATCGGATCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGa                                                                         <  2:380624/149‑1 (MQ=255)
                                                                           ccGCAGCACGGTTTGCCGCCGGCGAAATCTGAACCTTTACTGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGa                                                                         <  1:1105453/149‑1 (MQ=255)
                                                                              cagaaCTGTTTGCGGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACTCCATGTTCCAGATTGCGCTGGACTa                                                                      <  1:247519/149‑1 (MQ=255)
                                                                               agcaCGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGGGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTAc                                                                     <  1:1168691/149‑1 (MQ=255)
                                                                                 caCGGTTTGCTGCTGTCGAAAGCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGAAGAGAAAAAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGa                                                                   <  2:797717/149‑1 (MQ=255)
                                                                                           tgctgTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTc                                                         >  2:840175/1‑149 (MQ=255)
                                                                                           tgctgTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACAAGTGTTAGTGGCGCTATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAACGGCGTc                                                         >  2:443913/1‑149 (MQ=255)
                                                                                             ctgTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGcgacga                                                       >  1:599626/1‑149 (MQ=255)
                                                                                             ctgTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGa                                                       >  1:893372/1‑149 (MQ=255)
                                                                                             ctgTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGa                                                       >  1:662880/1‑149 (MQ=255)
                                                                                             ctgTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGGCGa                                                       >  1:1033694/1‑149 (MQ=255)
                                                                                                tCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGCCAAGAGCGTCGatgt                                                    >  2:22622/1‑147 (MQ=255)
                                                                                                      tCTGAACCTTTGATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCt                                              <  2:110061/149‑1 (MQ=255)
                                                                                                          aaCCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCAc                                          <  1:498108/149‑1 (MQ=255)
                                                                                                             cggtAAGGCGTGAAGCGGTGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATGGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGaa                                       <  2:109124/146‑1 (MQ=255)
                                                                                                              tttAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGaaa                                      <  2:292373/149‑1 (MQ=255)
                                                                                                              tttAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGaaa                                      >  1:101428/1‑149 (MQ=255)
                                                                                                                 aaTGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCa                                   <  1:476755/149‑1 (MQ=255)
                                                                                                                    gCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACCCCATGTTCTAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAGACCACTc                                <  1:1163947/149‑1 (MQ=255)
                                                                                                                        gAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGc                            >  2:404198/1‑149 (MQ=255)
                                                                                                                        aaagggcggCAGGCCGGGGCGCACTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGc                            <  1:1199553/141‑1 (MQ=255)
                                                                                                                            cGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATGGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACGCCGGCAGGc                        >  2:139715/1‑149 (MQ=255)
                                                                                                                                 cAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGAGGc                   >  2:1235279/1‑149 (MQ=255)
                                                                                                                                  aGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCa                  >  1:437317/1‑149 (MQ=255)
                                                                                                                                   ggCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCCCCTGCACGAAACCACTCCGGCAGGCGTGGCAg                 >  1:690356/1‑149 (MQ=255)
                                                                                                                                   ggCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAg                 >  1:638340/1‑149 (MQ=255)
                                                                                                                                    gCGGGGTCGCATTCCGTCCGCGGCCTGGACCCCACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCACTTTGCGCTGCACTACGACAAAGGCGTCTACATCCACCTGCCCGAAACCACTCCCGCATGCGTGGCAGc                <  2:1154221/149‑1 (MQ=255)
                                                                                                                                    gCGGGGGGGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGc                <  2:266171/149‑1 (MQ=255)
                                                                                                                                       ggggCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCTAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCAt             <  1:320708/149‑1 (MQ=255)
                                                                                                                                         ggCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCa           <  1:1006254/149‑1 (MQ=255)
                                                                                                                                             cATTACGTCGGCGGCCTGGACCCGACCAGTGTGGCTGGCGCGATGGAAAAAGCCCTCGACACCATGTTCCAGATGGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATta       <  1:614925/149‑1 (MQ=255)
                                                                                                                                                  cGTCGGCGGCCTGGACCCGACCAGTGTTGCTGGCGCGATGGAAAAATCCCTCGACACCATGTTACCGAGTTCGCTGGACCAACACAAAGTCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATgg  <  1:292959/149‑1 (MQ=255)
                                                                                                                                                  cGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGCCAGGCGTGGCAGCCATCGATTATATga  >  1:508545/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGAAAAATCTGCAAAATTTGCAGGCGGTGCTGGCGCGACGTCAGGCGGGCTTTGAGTGTGAAATCGTCGCCTTCCCGCAGCACGGTTTGCTGCTGTCGAAATCTGAACCTTTAATGCGTGAAGCGATGCAGGCGGGGGCGCATTACGTCGGCGGCCTGGACCCGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGG  >  CP009273/337412‑337706

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