Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,858,963 A→T 19.2% V224E (GTA→GAA)  yidE ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,858,9630AT19.2% 610.1 / 137.3 250V224E (GTA→GAA) yidEputative transporter
Reads supporting (aligned to +/- strand):  ref base A (116/85);  new base T (22/26);  total (139/111)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 253 total aligned reads displayed.

TCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATG  >  CP009273/3858814‑3859112
                                                                                                                                                     |                                                                                                                                                     
tCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTa                                                                                                                                                       <  1:1484429/150‑1 (MQ=255)
 cACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTAc                                                                                                                                                      <  2:2526507/150‑1 (MQ=255)
  aCCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACa                                                                                                                                                     >  1:2045318/1‑150 (MQ=255)
   ccAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCAt                                                                                                                                                    >  2:3829803/1‑150 (MQ=255)
          cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAAt                                                                                                                                             <  1:625541/150‑1 (MQ=255)
          cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAAt                                                                                                                                             >  2:4312198/1‑150 (MQ=255)
           gcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATg                                                                                                                                            >  2:2609909/1‑150 (MQ=255)
                aTGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaa                                                                                                                                       >  1:976430/1‑150 (MQ=255)
                  gTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAAt                                                                                                                                     <  1:564271/150‑1 (MQ=255)
                     aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                                                                                                                  <  1:986586/150‑1 (MQ=255)
                     aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                                                                                                                  >  1:2416294/1‑150 (MQ=255)
                      aTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATg                                                                                                                                 >  1:3379113/1‑150 (MQ=255)
                        ccGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCa                                                                                                                               <  2:1412134/150‑1 (MQ=255)
                           ctggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGt                                                                                                                            >  1:277709/1‑150 (MQ=255)
                            tggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGtt                                                                                                                           >  2:3363290/1‑150 (MQ=255)
                            tggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGtt                                                                                                                           >  1:1385004/1‑150 (MQ=255)
                               ctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAgg                                                                                                                        <  2:3004902/150‑1 (MQ=255)
                                tgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAggg                                                                                                                       >  2:1323294/1‑150 (MQ=255)
                                 gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                      >  2:18122/1‑150 (MQ=255)
                                 gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                      >  1:3661669/1‑150 (MQ=255)
                                 gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                      <  2:1303271/150‑1 (MQ=255)
                                  aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGtt                                                                                                                     <  2:4030032/150‑1 (MQ=255)
                                  aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGtt                                                                                                                     <  2:1138178/150‑1 (MQ=255)
                                   cccACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtc                                                                                                                    >  2:2858444/1‑150 (MQ=255)
                                   cccACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtc                                                                                                                    >  2:1778215/1‑150 (MQ=255)
                                    ccaccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                                                                   >  1:2607558/1‑150 (MQ=255)
                                     caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                                                                  <  1:3648116/150‑1 (MQ=255)
                                     caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                                                                  <  1:231887/150‑1 (MQ=255)
                                     caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                                                                  <  1:1681822/150‑1 (MQ=255)
                                           aTGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcgcg                                                                                                            <  2:4362094/150‑1 (MQ=255)
                                             gcagcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGaa                                                                                                          >  2:4240648/1‑150 (MQ=255)
                                               agcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                        >  2:450721/1‑150 (MQ=255)
                                                gcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAatat                                                                                                       <  1:2150337/150‑1 (MQ=255)
                                                 caAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATAtt                                                                                                      <  2:2271383/150‑1 (MQ=255)
                                                  aaaTCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTg                                                                                                     >  2:3758274/1‑150 (MQ=255)
                                                     tCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATa                                                                                                  <  2:899433/150‑1 (MQ=255)
                                                     tCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATa                                                                                                  >  1:3958217/1‑150 (MQ=255)
                                                      cGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCACCGCGAATATTGATAg                                                                                                 >  2:269542/1‑150 (MQ=255)
                                                      cGCCCAGTTGGATAATGGAATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCCCCGTAGAGAATCGGTACAACTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGTATatagatag                                                                                                 >  2:3472620/1‑150 (MQ=255)
                                                       gCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGt                                                                                                >  1:444992/1‑150 (MQ=255)
                                                        cccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTc                                                                                               >  1:3167906/1‑150 (MQ=255)
                                                         ccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCt                                                                                              >  1:82010/1‑150 (MQ=255)
                                                         ccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCt                                                                                              >  2:4187767/1‑150 (MQ=255)
                                                          cAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt                                                                                             <  2:4544252/150‑1 (MQ=255)
                                                          cAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt                                                                                             <  2:4544120/150‑1 (MQ=255)
                                                           aGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTg                                                                                            <  1:1401672/150‑1 (MQ=255)
                                                           aGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTg                                                                                            <  1:3045727/150‑1 (MQ=255)
                                                              tGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAAGTTTGTCGCCGTTGAGTATCGGTACATCTTTAATGGCTAATTCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATc                                                                                         >  2:4268976/1‑150 (MQ=255)
                                                                   aaTGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcgc                                                                                    <  1:1034139/150‑1 (MQ=255)
                                                                           cTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCAt                                                                            >  1:70849/1‑150 (MQ=255)
                                                                            tGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCTGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCAtt                                                                           <  2:3650140/150‑1 (MQ=255)
                                                                               cGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGTTAGTCTTGATCAGCGCGCCGCCATTGGt                                                                        >  2:2877616/1‑150 (MQ=255)
                                                                                gAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTg                                                                       <  1:878550/150‑1 (MQ=255)
                                                                                gAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTg                                                                       <  2:1754460/150‑1 (MQ=255)
                                                                                 aaGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGc                                                                      <  2:1350961/150‑1 (MQ=255)
                                                                                   ggAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGt                                                                    <  2:2288235/150‑1 (MQ=255)
                                                                                     aaCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGa                                                                  <  1:715300/150‑1 (MQ=255)
                                                                                         tttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTGAGGGTTCTCAACGCGACTATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGAc                                                              >  2:3399393/1‑150 (MQ=255)
                                                                                         tttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGAc                                                              >  1:1922138/1‑150 (MQ=255)
                                                                                           tAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTc                                                            <  1:1711284/150‑1 (MQ=255)
                                                                                           tAGGGTTTCTTCGCGTTTAAGACGCGATCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCTGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATAAGCGCGCCGCCATTGGTGCGTGAAGACTc                                                            <  2:2136443/150‑1 (MQ=255)
                                                                                            aGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCg                                                           >  1:4307137/1‑150 (MQ=255)
                                                                                            aGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCg                                                           >  1:4173572/1‑150 (MQ=255)
                                                                                                   cttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTgctgct                                                    <  1:467361/150‑1 (MQ=255)
                                                                                                    ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   >  2:1554153/1‑150 (MQ=255)
                                                                                                     tcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGa                                                  >  2:1371047/1‑150 (MQ=255)
                                                                                                      cgcgTTTCAGACGCGATCAGATAATTTTTTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGAAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAg                                                 >  2:4359925/1‑150 (MQ=255)
                                                                                                      cgcgTTTCAGACGCGAGCAGATTATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAg                                                 <  1:2728592/150‑1 (MQ=255)
                                                                                                      cgcgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAg                                                 <  2:2118896/150‑1 (MQ=255)
                                                                                                        cgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                               >  2:1524775/1‑150 (MQ=255)
                                                                                                        cgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                               >  2:3072180/1‑150 (MQ=255)
                                                                                                         gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTCCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCtt                                              >  1:3295197/1‑150 (MQ=255)
                                                                                                         gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCtt                                              >  1:2496696/1‑150 (MQ=255)
                                                                                                              aGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtc                                         >  2:1223070/1‑150 (MQ=255)
                                                                                                               gACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtct                                        <  1:3140649/150‑1 (MQ=255)
                                                                                                                aCGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtctc                                       <  2:2062026/150‑1 (MQ=255)
                                                                                                                 cgcgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCg                                      <  2:3339886/150‑1 (MQ=255)
                                                                                                                   cgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGAc                                    <  2:4486225/150‑1 (MQ=255)
                                                                                                                   cgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGAc                                    >  2:3571830/1‑150 (MQ=255)
                                                                                                                    gAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACa                                   <  2:617819/150‑1 (MQ=255)
                                                                                                                       cAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                                <  1:3710453/150‑1 (MQ=255)
                                                                                                                       cAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                                >  2:3885740/1‑150 (MQ=255)
                                                                                                                       cAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                                >  2:797196/1‑150 (MQ=255)
                                                                                                                           tAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAcgc                            >  2:242804/1‑150 (MQ=255)
                                                                                                                               tttGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGaa                        >  2:823665/1‑150 (MQ=255)
                                                                                                                                   tCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGACAATa                    >  2:4637039/1‑150 (MQ=255)
                                                                                                                                    cGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATaa                   <  1:4045162/150‑1 (MQ=255)
                                                                                                                                     gCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAAc                  <  2:2965665/150‑1 (MQ=255)
                                                                                                                                       cGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAAccc                <  1:2139866/150‑1 (MQ=255)
                                                                                                                                          tGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCa             >  1:2477933/1‑150 (MQ=255)
                                                                                                                                           gagaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCaa            >  1:1205066/1‑150 (MQ=255)
                                                                                                                                             gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa          >  2:2507525/1‑150 (MQ=255)
                                                                                                                                              aaTCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACAt         >  1:2715584/1‑150 (MQ=255)
                                                                                                                                               aTCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATc        >  2:1722583/1‑150 (MQ=255)
                                                                                                                                               aTCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATc        >  1:1101568/1‑150 (MQ=255)
                                                                                                                                                 cGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCca      >  2:1944874/1‑150 (MQ=255)
                                                                                                                                                  ggTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCcac     >  1:2831749/1‑150 (MQ=255)
                                                                                                                                                   gTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCcaca    >  2:3265653/1‑150 (MQ=255)
                                                                                                                                                    tACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAt   >  2:3084656/1‑150 (MQ=255)
                                                                                                                                                     aCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATg  >  1:2790550/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATG  >  CP009273/3858814‑3859112

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