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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,858,9630AT59.7% 207.0 / inf 263V224E (GTA→GAA) yidEputative transporter
Reads supporting (aligned to +/- strand):  ref base A (52/53);  new base T (76/81);  total (129/134)
Fisher's exact test for biased strand distribution p-value = 9.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 264 total aligned reads displayed.

ACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATG  >  CP009273/3858816‑3859112
                                                                                                                                                   |                                                                                                                                                     
 ccAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCAt                                                                                                                                                    >  2:1166528/1‑150 (MQ=255)
    gTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCtt                                                                                                                                                 <  1:2166194/150‑1 (MQ=255)
     ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCttt                                                                                                                                                <  2:2301403/150‑1 (MQ=255)
     ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCttt                                                                                                                                                <  1:847439/150‑1 (MQ=255)
     ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTCGAGAATCGGTTCATCttt                                                                                                                                                <  1:4156217/150‑1 (MQ=255)
        cgcaaaaTGAAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGACCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAAt                                                                                                                                             >  1:765808/7‑150 (MQ=255)
        cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAAt                                                                                                                                             <  1:1112933/150‑1 (MQ=255)
        cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAAt                                                                                                                                             <  2:4075543/150‑1 (MQ=255)
        cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAAt                                                                                                                                             >  1:4390580/1‑150 (MQ=255)
        cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAAt                                                                                                                                             >  1:4262282/1‑150 (MQ=255)
        cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAAt                                                                                                                                             <  1:4497186/150‑1 (MQ=255)
          cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTTGATAATGGTATCTGGCGAATGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATgg                                                                                                                                           <  2:2309867/150‑1 (MQ=255)
          cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATgg                                                                                                                                           >  1:1582534/1‑150 (MQ=255)
          cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATgg                                                                                                                                           >  1:1582124/1‑150 (MQ=255)
          cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATgg                                                                                                                                           >  2:2898571/1‑150 (MQ=255)
               tGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCNTTGAGAATCGGTACATCTTTAATGGCTaaa                                                                                                                                      <  1:4286864/150‑1 (MQ=255)
               tGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaaa                                                                                                                                      >  2:4507500/1‑150 (MQ=255)
                 tAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATc                                                                                                                                    >  1:1171232/1‑150 (MQ=255)
                   aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                                                                                                                  >  1:1662119/1‑150 (MQ=255)
                          tggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGtt                                                                                                                           <  1:3335280/150‑1 (MQ=255)
                              tgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAggg                                                                                                                       >  1:4150063/1‑150 (MQ=255)
                              tgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAggg                                                                                                                       >  2:2235844/1‑150 (MQ=255)
                              tgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAggg                                                                                                                       >  1:2827434/1‑150 (MQ=255)
                                aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGtt                                                                                                                     <  1:554681/150‑1 (MQ=255)
                                aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGtt                                                                                                                     <  2:2930201/150‑1 (MQ=255)
                                  ccaccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                                                                   <  2:1803079/150‑1 (MQ=255)
                                   caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                                                                  >  2:3948369/1‑150 (MQ=255)
                                    accaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTCGGGTTCTCa                                                                                                                 >  2:3315897/1‑150 (MQ=255)
                                     ccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCaa                                                                                                                <  2:4305309/150‑1 (MQ=255)
                                      caGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                               <  2:2129950/150‑1 (MQ=255)
                                      caGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                               <  2:2684007/150‑1 (MQ=255)
                                        gATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcgc                                                                                                             >  2:1357772/1‑150 (MQ=255)
                                        gATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcgc                                                                                                             <  2:3811437/150‑1 (MQ=255)
                                                aaaTCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTg                                                                                                     >  1:3513166/1‑150 (MQ=255)
                                                      cccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTc                                                                                               >  1:273399/1‑150 (MQ=255)
                                                       ccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCt                                                                                              >  1:4026385/1‑150 (MQ=255)
                                                       ccAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCt                                                                                              <  1:4577080/150‑1 (MQ=255)
                                                             ggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCa                                                                                        <  1:2397187/150‑1 (MQ=255)
                                                             ggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCa                                                                                        >  2:782644/1‑150 (MQ=255)
                                                              gATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAg                                                                                       <  1:4078723/150‑1 (MQ=255)
                                                               aTAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgc                                                                                      >  1:2850936/1‑150 (MQ=255)
                                                               aTAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgc                                                                                      >  1:441535/1‑150 (MQ=255)
                                                                tAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcg                                                                                     >  2:3071334/1‑150 (MQ=255)
                                                                tAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcg                                                                                     <  2:2451341/150‑1 (MQ=255)
                                                                tAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcg                                                                                     <  1:584486/150‑1 (MQ=255)
                                                                    ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                 <  2:1235682/150‑1 (MQ=255)
                                                                    ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                 <  2:51857/150‑1 (MQ=255)
                                                                        tCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCa                                                                             <  2:1582478/150‑1 (MQ=255)
                                                                        tCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCa                                                                             >  2:2844581/1‑150 (MQ=255)
                                                                        tCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCa                                                                             >  1:3561852/1‑150 (MQ=255)
                                                                           ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTg                                                                          >  2:898974/1‑150 (MQ=255)
                                                                           ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTg                                                                          <  1:4219030/150‑1 (MQ=255)
                                                                           ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCAATTg                                                                          <  1:1357772/150‑1 (MQ=255)
                                                                             cGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGt                                                                        <  1:132984/150‑1 (MQ=255)
                                                                               aaGGAACTTTAAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTAGAGAATCGGTACATCATTAATGGCTAAAACATGCATGATAGGGTACTCAACCCGAATAATGATAGTCTTGATCAGCGCGCCGCCAttcgaga                                                                      >  1:2503835/1‑145 (MQ=255)
                                                                                   aaCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGa                                                                  >  2:272177/1‑150 (MQ=255)
                                                                                    aCTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGaa                                                                 <  2:4172116/150‑1 (MQ=255)
                                                                                     cTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAg                                                                >  2:2395297/1‑150 (MQ=255)
                                                                                     cTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAAGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAg                                                                >  2:4146384/1‑150 (MQ=255)
                                                                                     cTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAg                                                                >  1:2540199/1‑150 (MQ=255)
                                                                                         tAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTc                                                            <  2:442782/150‑1 (MQ=255)
                                                                                         tAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTc                                                            <  2:3194304/150‑1 (MQ=255)
                                                                                          aGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCg                                                           >  2:2973410/1‑150 (MQ=255)
                                                                                          aGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCg                                                           >  2:4390907/1‑150 (MQ=255)
                                                                                              tttCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtgct                                                       <  1:4581910/150‑1 (MQ=255)
                                                                                              tttCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtgct                                                       <  1:1413208/150‑1 (MQ=255)
                                                                                                  ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCATCGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   <  1:4419307/150‑1 (MQ=255)
                                                                                                  ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   <  1:3088536/150‑1 (MQ=255)
                                                                                                  ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   <  1:4643519/150‑1 (MQ=255)
                                                                                                  ttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGctgctg                                                   <  2:1113371/150‑1 (MQ=255)
                                                                                                     gcgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGc                                                >  2:1441656/1‑150 (MQ=255)
                                                                                                     gcgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGc                                                <  2:3567910/150‑1 (MQ=255)
                                                                                                      cgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                               <  1:1724350/150‑1 (MQ=255)
                                                                                                      cgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                               <  1:3858650/150‑1 (MQ=255)
                                                                                                      cgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                               <  1:3858762/150‑1 (MQ=255)
                                                                                                      cgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                               >  2:3700657/1‑150 (MQ=255)
                                                                                                            aGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACGTCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtc                                         >  2:3342258/1‑150 (MQ=255)
                                                                                                            aGACGCGAGCAGATAATTTTGTAGCCGTTGAGAATCGGTTCATCATTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGtc                                         >  1:3447392/1‑150 (MQ=255)
                                                                                                                gcgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGa                                     >  2:2488401/1‑150 (MQ=255)
                                                                                                                 cgAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGAc                                    >  1:3359599/1‑150 (MQ=255)
                                                                                                                    gCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACAtt                                 <  1:2308351/150‑1 (MQ=255)
                                                                                                                     cAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                                >  1:1623699/1‑150 (MQ=255)
                                                                                                                      aGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGa                               >  1:1008457/1‑150 (MQ=255)
                                                                                                                      aGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGa                               <  1:2429802/150‑1 (MQ=255)
                                                                                                                        aTAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAcg                             >  2:3553472/1‑150 (MQ=255)
                                                                                                                          aaTTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTAGTGCTGCTGAGCTTCTGTCTCGACATTGAcgcg                           >  1:1542310/1‑150 (MQ=255)
                                                                                                                          aaTTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAcgcg                           <  2:3513166/150‑1 (MQ=255)
                                                                                                                           aTTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCgg                          >  1:1060955/1‑150 (MQ=255)
                                                                                                                                   gCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAAc                  >  1:2415546/1‑150 (MQ=255)
                                                                                                                                     cGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAAccc                >  2:600587/1‑150 (MQ=255)
                                                                                                                                       ttGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGc              >  2:3399856/1‑150 (MQ=255)
                                                                                                                                        tGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCa             >  2:372240/1‑150 (MQ=255)
                                                                                                                                          agaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAAc           <  2:203942/150‑1 (MQ=255)
                                                                                                                                           gaATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa          >  1:2994297/1‑150 (MQ=255)
                                                                                                                                           gaATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa          <  2:78304/150‑1 (MQ=255)
                                                                                                                                           gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa          <  1:2898571/150‑1 (MQ=255)
                                                                                                                                           gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa          <  1:3725818/150‑1 (MQ=255)
                                                                                                                                            aaTCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACAt         >  1:395038/1‑150 (MQ=255)
                                                                                                                                             aTCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATc        >  1:2156681/1‑150 (MQ=255)
                                                                                                                                                ggTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCcac     >  2:504147/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACATG  >  CP009273/3858816‑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: