Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,175,193 A→T 23.7% E36V (GAA→GTA)  cobB → deacetylase of acs and cheY, chemotaxis regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,175,1930AT23.7% 616.9 / 208.9 283E36V (GAA→GTA) cobBdeacetylase of acs and cheY, chemotaxis regulator
Reads supporting (aligned to +/- strand):  ref base A (116/100);  new base T (28/39);  total (144/139)
Fisher's exact test for biased strand distribution p-value = 9.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 288 total aligned reads displayed.

GCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGAT  >  CP009273/1175044‑1175341
                                                                                                                                                     |                                                                                                                                                    
gCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGa                                                                                                                                                      >  2:4031342/1‑150 (MQ=255)
 cGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGaa                                                                                                                                                     <  1:4683568/150‑1 (MQ=255)
  gTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAg                                                                                                                                                    >  2:4182151/1‑150 (MQ=255)
    cAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGAt                                                                                                                                                  >  1:4381764/1‑150 (MQ=255)
      gCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGAtgt                                                                                                                                                <  1:3702510/150‑1 (MQ=255)
            ttttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAAc                                                                                                                                          <  2:2209394/150‑1 (MQ=255)
             tttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACg                                                                                                                                         <  2:1455419/150‑1 (MQ=255)
             tttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACg                                                                                                                                         >  2:668630/1‑150 (MQ=255)
              ttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTCTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGACGAACATCGGGTTGAAGATGTGGCAACGc                                                                                                                                        >  2:4308169/1‑150 (MQ=255)
              ttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGc                                                                                                                                        >  1:1810655/1‑150 (MQ=255)
              ttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGc                                                                                                                                        <  2:1623464/150‑1 (MQ=255)
              ttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGc                                                                                                                                        <  2:2804151/150‑1 (MQ=255)
                ttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCg                                                                                                                                      >  1:3879332/1‑150 (MQ=255)
                   agagATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGaa                                                                                                                                   <  1:3679383/150‑1 (MQ=255)
                    gagaTAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAg                                                                                                                                  >  2:4519438/1‑150 (MQ=255)
                        tAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGttt                                                                                                                              >  1:4909589/1‑150 (MQ=255)
                         aaaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTc                                                                                                                             >  1:3934914/1‑150 (MQ=255)
                          aaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCg                                                                                                                            <  2:5098342/150‑1 (MQ=255)
                          aaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCg                                                                                                                            <  2:2554407/150‑1 (MQ=255)
                             tggtgCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATc                                                                                                                         <  1:3671523/150‑1 (MQ=255)
                              ggtgCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcg                                                                                                                        >  1:3539210/1‑150 (MQ=255)
                               gtgCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATcgc                                                                                                                       <  2:3618982/150‑1 (MQ=255)
                                 gCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGa                                                                                                                     >  2:1156766/1‑150 (MQ=255)
                                    ggAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATcc                                                                                                                  <  1:2199072/150‑1 (MQ=255)
                                       aGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATGGCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGa                                                                                                               <  2:4292763/150‑1 (MQ=255)
                                            tGGAAAAACCAAGAGTACTCGTACTGACAGGGGCATGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTgg                                                                                                          <  2:3070965/150‑1 (MQ=255)
                                            tGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCATGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTgg                                                                                                          >  2:45685/1‑150 (MQ=255)
                                              gAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTg                                                                                                        >  1:4371687/1‑150 (MQ=255)
                                              gAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTg                                                                                                        >  1:2543284/1‑150 (MQ=255)
                                               aaaaaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGc                                                                                                       <  1:618012/150‑1 (MQ=255)
                                                  aaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAg                                                                                                    <  2:1897780/150‑1 (MQ=255)
                                                  aaCCAAGAGTACTAGTACTGACAGGGGCAGGAATTTCTGCGTAATCAGGTATTCGTACATTTCGCGCCGCATATGGCCTGTGGGAAGTACATCGGGTTGAAGATGTGTCAACGCAGTTATTTTTCGATCGCGATCCTTAACTGGTGCAAg                                                                                                    <  2:3460640/150‑1 (MQ=255)
                                                     cAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGt                                                                                                 >  1:1682593/1‑150 (MQ=255)
                                                     cAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGt                                                                                                 >  1:2619365/1‑150 (MQ=255)
                                                      aaGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGtt                                                                                                >  1:3548774/1‑150 (MQ=255)
                                                        gagTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGtttt                                                                                              >  2:2260963/1‑150 (MQ=255)
                                                          gTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTtat                                                                                            <  2:2245496/150‑1 (MQ=255)
                                                            aCTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATaa                                                                                          >  2:2013836/1‑150 (MQ=255)
                                                                 tACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGccc                                                                                     <  1:2943831/150‑1 (MQ=255)
                                                                 tACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGccc                                                                                     >  2:2495206/1‑150 (MQ=255)
                                                                     gACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtcg                                                                                 <  2:2711376/150‑1 (MQ=255)
                                                                     gACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtcg                                                                                 <  2:111308/150‑1 (MQ=255)
                                                                     gACAGGGGCAGGAATGTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtcg                                                                                 <  2:3287841/150‑1 (MQ=255)
                                                                      aCAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtcgt                                                                                >  2:2079002/1‑150 (MQ=255)
                                                                            gCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAg                                                                          <  1:919350/150‑1 (MQ=255)
                                                                            gCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAg                                                                          <  1:34703/150‑1 (MQ=255)
                                                                             cAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGc                                                                         >  2:1642314/1‑150 (MQ=255)
                                                                                  aTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTgcag                                                                    <  1:1013930/150‑1 (MQ=255)
                                                                                   tttCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTgcagc                                                                   <  1:4492070/150‑1 (MQ=255)
                                                                                      cTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCagcagc                                                                >  1:706005/1‑150 (MQ=255)
                                                                                      cTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCagcagc                                                                <  1:2915908/150‑1 (MQ=255)
                                                                                         cGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAg                                                             >  1:1539070/1‑150 (MQ=255)
                                                                                          ggAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGa                                                            <  1:1808924/150‑1 (MQ=255)
                                                                                           gAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGaa                                                           >  2:1016511/1‑150 (MQ=255)
                                                                                            aaTCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGaaa                                                          >  2:2625104/1‑150 (MQ=255)
                                                                                            aaTCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGaaa                                                          >  2:1511569/1‑150 (MQ=255)
                                                                                             aTCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGACGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAAt                                                         >  2:4703265/1‑150 (MQ=255)
                                                                                                aGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCa                                                      >  1:1969523/1‑150 (MQ=255)
                                                                                                aGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCa                                                      >  2:3781863/1‑150 (MQ=255)
                                                                                                   tATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGcc                                                   >  1:158987/1‑150 (MQ=255)
                                                                                                    aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCt                                                  >  2:3029443/1‑150 (MQ=255)
                                                                                                     ttCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTa                                                 >  2:2989100/1‑150 (MQ=255)
                                                                                                       cGTACCTTTCGCGCCGCAGCTGGCCTGGGGGACGAACATCGGGTTGAAGAGGTGGCACCGCCGGAAGGGTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTGATCACGCCCGTCGTCGACAGCTGCCGCCGCCAGACCGTCAGCCTCAc                                               >  2:2886874/1‑150 (MQ=255)
                                                                                                       cGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAc                                               >  1:376849/1‑150 (MQ=255)
                                                                                                        gTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAcg                                              >  2:1684202/1‑150 (MQ=255)
                                                                                                         tACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAcgc                                             >  1:2338000/1‑150 (MQ=255)
                                                                                                          aCCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAcgcc                                            <  2:3901048/150‑1 (MQ=255)
                                                                                                          aCCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAcgcc                                            >  1:597722/1‑150 (MQ=255)
                                                                                                           ccTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAcgccg                                           >  2:1458848/1‑150 (MQ=255)
                                                                                                             tttCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCcgcg                                         >  1:4313371/1‑150 (MQ=255)
                                                                                                             tttCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCcgcg                                         <  2:650948/150‑1 (MQ=255)
                                                                                                             tttCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCcgcg                                         <  2:5025724/150‑1 (MQ=255)
                                                                                                                cgcgCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGACGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCAt                                      >  2:4103756/1‑150 (MQ=255)
                                                                                                                cgcgCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCAt                                      <  2:4226754/150‑1 (MQ=255)
                                                                                                                  cgccgcAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCt                                    >  2:3158437/1‑150 (MQ=255)
                                                                                                                        aGATGGCCTGGGGGAAGCCCATCGGGGTGTAGCTGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATCCCGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCt                              <  2:3502845/150‑1 (MQ=255)
                                                                                                                          aTGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTgg                            <  2:1568779/150‑1 (MQ=255)
                                                                                                                           tGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGc                           >  1:2135013/1‑150 (MQ=255)
                                                                                                                            ggCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCt                          <  1:1115937/150‑1 (MQ=255)
                                                                                                                            ggCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCt                          <  2:1979790/150‑1 (MQ=255)
                                                                                                                              ccTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTaa                        >  2:1296825/1‑150 (MQ=255)
                                                                                                                               cTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTaaa                       >  1:3752752/1‑150 (MQ=255)
                                                                                                                                  tgGGAAGTACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTg                    >  1:4606934/1‑150 (MQ=255)
                                                                                                                                  tgGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTg                    <  1:3058952/150‑1 (MQ=255)
                                                                                                                                   gggAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGc                   >  1:2786060/1‑150 (MQ=255)
                                                                                                                                    ggAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCa                  <  2:4091675/150‑1 (MQ=255)
                                                                                                                                     gaagaaCATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCaa                 >  1:1334294/1‑150 (MQ=255)
                                                                                                                                     gaagaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCaa                 <  1:4651127/150‑1 (MQ=255)
                                                                                                                                      aagaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAg                <  2:3813179/150‑1 (MQ=255)
                                                                                                                                       agaaCATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGa               <  1:337542/150‑1 (MQ=255)
                                                                                                                                       agaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGa               >  2:3630919/1‑150 (MQ=255)
                                                                                                                                       agaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGa               >  1:1125841/1‑150 (MQ=255)
                                                                                                                                         aaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACg             <  1:2687810/150‑1 (MQ=255)
                                                                                                                                         aaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACg             <  1:3260881/150‑1 (MQ=255)
                                                                                                                                            aTCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGccc          >  1:227909/1‑150 (MQ=255)
                                                                                                                                             tCGGGTTGAAGCTGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGATCTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCt         <  1:1366162/150‑1 (MQ=255)
                                                                                                                                              cGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTc        <  1:4031367/150‑1 (MQ=255)
                                                                                                                                                ggTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCgg      <  2:1802691/150‑1 (MQ=255)
                                                                                                                                                   tGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGa   >  1:3918802/1‑150 (MQ=255)
                                                                                                                                                   tGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGa   >  1:4848266/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGAT  >  CP009273/1175044‑1175341

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: