Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,009,389 G→T 8.7% A47A (GCC→GCA ysgA ← putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,009,3890GT 8.7% 1179.4 / 54.2 321A47A (GCC→GCAysgAputative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base G (145/148);  new base T (17/11);  total (162/159)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 329 total aligned reads displayed.

TCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTC  >  CP009273/4009240‑4009538
                                                                                                                                                     |                                                                                                                                                     
  gCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGc                                                                                                                                                     >  1:439077/1‑150 (MQ=255)
  gCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGc                                                                                                                                                     <  1:572767/150‑1 (MQ=255)
   ccTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAGTTACCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCa                                                                                                                                                    >  2:4185433/1‑150 (MQ=255)
    cTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCAt                                                                                                                                                   <  2:3918231/150‑1 (MQ=255)
         cgGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTaa                                                                                                                                              <  1:3082756/150‑1 (MQ=255)
             aGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAg                                                                                                                                          >  1:4542568/1‑150 (MQ=255)
             aGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAg                                                                                                                                          <  2:286160/150‑1 (MQ=255)
             aGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGAAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAg                                                                                                                                          <  1:2926188/150‑1 (MQ=255)
              gTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAgg                                                                                                                                         <  2:2058779/150‑1 (MQ=255)
                aaaGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCa                                                                                                                                       >  1:1918023/1‑150 (MQ=255)
                 aaGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCAt                                                                                                                                      <  2:4184933/150‑1 (MQ=255)
                  aGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATg                                                                                                                                     >  1:4231293/1‑150 (MQ=255)
                     tCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTa                                                                                                                                  >  1:1496753/1‑150 (MQ=255)
                      cAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTAt                                                                                                                                 <  2:887029/150‑1 (MQ=255)
                       aGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATc                                                                                                                                <  2:1302525/150‑1 (MQ=255)
                       aGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATc                                                                                                                                <  2:1590872/150‑1 (MQ=255)
                         gTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATccc                                                                                                                              >  1:665739/1‑150 (MQ=255)
                         gTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATccc                                                                                                                              <  2:1888854/150‑1 (MQ=255)
                          tGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATcccc                                                                                                                             <  1:1990479/150‑1 (MQ=255)
                           gCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCt                                                                                                                            <  2:2915083/150‑1 (MQ=255)
                           gCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCt                                                                                                                            >  1:1754804/1‑150 (MQ=255)
                            cGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTCAGCAGGCATGTTATCCCCtt                                                                                                                           >  2:3730676/1‑150 (MQ=255)
                             gATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTtg                                                                                                                          >  2:4144996/1‑150 (MQ=255)
                                aGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGa                                                                                                                       <  1:1003623/150‑1 (MQ=255)
                                 gCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGaa                                                                                                                      >  1:130423/1‑150 (MQ=255)
                                 gCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGaa                                                                                                                      >  1:2788704/1‑150 (MQ=255)
                                  ccAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAg                                                                                                                     <  2:2803132/150‑1 (MQ=255)
                                   cAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAgg                                                                                                                    >  1:4356339/1‑150 (MQ=255)
                                       tACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATa                                                                                                                <  2:1625571/150‑1 (MQ=255)
                                        aCCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATag                                                                                                               <  2:2035242/150‑1 (MQ=255)
                                        aCCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATag                                                                                                               <  2:1156896/150‑1 (MQ=255)
                                          cccTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATagag                                                                                                             <  1:504821/150‑1 (MQ=255)
                                              ccAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGa                                                                                                         <  1:2465462/150‑1 (MQ=255)
                                                aGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAg                                                                                                       >  1:3755345/1‑150 (MQ=255)
                                                   gcCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGccc                                                                                                    <  2:1020050/150‑1 (MQ=255)
                                                     cAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGc                                                                                                  >  2:1509951/1‑150 (MQ=255)
                                                      tgACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGcc                                                                                                 <  1:2106427/149‑1 (MQ=255)
                                                      aGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGcc                                                                                                 <  2:4193547/150‑1 (MQ=255)
                                                         cGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCcaca                                                                                              >  1:482731/1‑150 (MQ=255)
                                                          ggCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACaa                                                                                             >  1:2571255/1‑150 (MQ=255)
                                                                 aaTGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGcgtc                                                                                      >  1:4472739/1‑150 (MQ=255)
                                                                  aTGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGcgtcg                                                                                     <  2:4029623/150‑1 (MQ=255)
                                                                         ggATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCgggg                                                                              >  2:1470946/1‑150 (MQ=255)
                                                                         ggATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCgggg                                                                              >  1:4050498/1‑150 (MQ=255)
                                                                            tatGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCt                                                                           >  1:4166444/1‑150 (MQ=255)
                                                                            tatGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCt                                                                           <  1:2245019/150‑1 (MQ=255)
                                                                               gTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGaa                                                                        >  2:2253746/1‑150 (MQ=255)
                                                                                 tCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACg                                                                      <  1:877776/150‑1 (MQ=255)
                                                                                       aCGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTc                                                                <  2:4075122/150‑1 (MQ=255)
                                                                                        cGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCg                                                               <  2:4071149/150‑1 (MQ=255)
                                                                                        cGCCAAAAATTTCCTGCACTACAATGACCACAGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCg                                                               >  1:3572593/1‑150 (MQ=255)
                                                                                             aaaaTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCg                                                          >  1:3245444/1‑150 (MQ=255)
                                                                                             aaaaTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCg                                                          >  1:3937819/1‑150 (MQ=255)
                                                                                             aaaaTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCg                                                          <  1:3282703/150‑1 (MQ=255)
                                                                                             aaaaTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCg                                                          <  1:2382110/150‑1 (MQ=255)
                                                                                               aaTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCgag                                                        <  2:1031674/150‑1 (MQ=255)
                                                                                                aTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCgaga                                                       <  1:1074370/150‑1 (MQ=255)
                                                                                                aTTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCgaga                                                       <  1:1244384/150‑1 (MQ=255)
                                                                                                 tttCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGagag                                                      >  2:331526/1‑150 (MQ=255)
                                                                                                 tttCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGagag                                                      >  1:1855638/1‑150 (MQ=255)
                                                                                                   tCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGc                                                    <  2:595273/150‑1 (MQ=255)
                                                                                                   accTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGc                                                    <  1:4185433/149‑1 (MQ=255)
                                                                                                     cTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACACTCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGa                                                  >  2:3497355/1‑150 (MQ=255)
                                                                                                      tGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGAt                                                 >  2:3574376/1‑150 (MQ=255)
                                                                                                          cTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcag                                             <  1:669604/150‑1 (MQ=255)
                                                                                                            aCAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATgcagca                                           >  2:3618233/1‑150 (MQ=255)
                                                                                                             cAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                          <  2:4144061/150‑1 (MQ=255)
                                                                                                             cAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGcagcag                                          <  2:1891713/150‑1 (MQ=255)
                                                                                                              aaTGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAgg                                         >  2:502646/1‑150 (MQ=255)
                                                                                                               aTGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGt                                        <  1:2586504/150‑1 (MQ=255)
                                                                                                                 gACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGc                                      >  1:3139106/1‑150 (MQ=255)
                                                                                                                     aCTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAAt                                  >  2:255078/1‑150 (MQ=255)
                                                                                                                       tGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATcc                                <  1:121686/150‑1 (MQ=255)
                                                                                                                        ggCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                               >  2:245554/1‑150 (MQ=255)
                                                                                                                        cgaaGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCACCTTGTGAAGGGATAGAGGTGAAGCCCTCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCa                               >  2:3388592/4‑150 (MQ=255)
                                                                                                                             tGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttg                          <  2:26752/150‑1 (MQ=255)
                                                                                                                             tGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttg                          >  1:1472209/1‑150 (MQ=255)
                                                                                                                              gggCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgt                         <  2:1833201/150‑1 (MQ=255)
                                                                                                                                gCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGAttgttg                       >  1:3630016/1‑150 (MQ=255)
                                                                                                                                  cATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtg                     <  2:4362496/150‑1 (MQ=255)
                                                                                                                                  cATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtg                     <  1:2081697/150‑1 (MQ=255)
                                                                                                                                  cATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtg                     <  2:2309899/150‑1 (MQ=255)
                                                                                                                                   aTCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                    <  2:856976/150‑1 (MQ=255)
                                                                                                                                   aTCGCTTTGCTTTGGTCTTGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                    <  2:130423/150‑1 (MQ=255)
                                                                                                                                   aTCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgt                    >  2:2612055/1‑150 (MQ=255)
                                                                                                                                    tCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATtgttgtgtt                   <  1:3882719/150‑1 (MQ=255)
                                                                                                                                       ctttgctttgGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTgttgtt                <  1:3395174/150‑1 (MQ=255)
                                                                                                                                        tttgctttgGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGttgttg               >  2:3415431/1‑150 (MQ=255)
                                                                                                                                         ttgctttgGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc              >  2:4231229/1‑150 (MQ=255)
                                                                                                                                         ttgctttgGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGc              >  2:4231292/1‑150 (MQ=255)
                                                                                                                                           gctttgGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCa            >  2:4261799/1‑150 (MQ=255)
                                                                                                                                               tgGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGt        >  1:423190/1‑150 (MQ=255)
                                                                                                                                                ggTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTa       <  2:3784514/150‑1 (MQ=255)
                                                                                                                                                ggTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTa       <  2:4530808/150‑1 (MQ=255)
                                                                                                                                                 gTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt      <  1:747880/150‑1 (MQ=255)
                                                                                                                                                 gTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt      <  2:1211923/150‑1 (MQ=255)
                                                                                                                                                 gTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt      <  1:3027421/150‑1 (MQ=255)
                                                                                                                                                 gTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTAt      >  2:3079934/1‑150 (MQ=255)
                                                                                                                                                   cTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtc    <  2:676929/150‑1 (MQ=255)
                                                                                                                                                     ggCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATtctc  >  1:4101824/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATCCCCTTGTGAAGGGATAGAGGTGAAGCCCGCCACAATCGCGTCGTCCGGGGTCTGAACGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTC  >  CP009273/4009240‑4009538

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: