Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 800,779 G→T 81.6% E607D (GAG→GAT ybhJ → putative hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273800,7790GT81.6% 824.1 / 218.1 358E607D (GAG→GATybhJputative hydratase
Reads supporting (aligned to +/- strand):  ref base G (23/43);  new base T (146/146);  total (169/189)
Fisher's exact test for biased strand distribution p-value = 2.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

Only 100 of 365 total aligned reads displayed.

CGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTC  >  CP009273/800632‑800924
                                                                                                                                                   |                                                                                                                                                 
 ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATgt                                                                                                                                                 >  2:4502578/1‑149 (MQ=255)
 ggTCAAACTTCTTCTTATCGTACAAATCCGTTTGGTCGGGGGGTTTTTACCGTGGCTCGCCGCGACCCGGTTTTTGTTAGCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgt                                                                                                                                                 <  1:582086/149‑1 (MQ=255)
  tttaaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCGGTTTATGTTAGCAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGgtg                                                                                                                                                <  1:3731942/146‑1 (MQ=255)
  gTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAtgtg                                                                                                                                                >  1:4321561/1‑149 (MQ=255)
     aaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGttt                                                                                                                                             >  2:2692207/1‑149 (MQ=255)
           cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcgc                                                                                                                                       >  1:3157077/1‑149 (MQ=255)
           cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcgc                                                                                                                                       >  1:2464853/1‑149 (MQ=255)
           cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTGATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcgc                                                                                                                                       >  2:4665083/1‑149 (MQ=255)
             tcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                     <  2:752428/149‑1 (MQ=255)
              cttACCGTTCAAATCCGATTGGTCTGGCGGAGTTGACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAtt                                                                                                                                    <  2:2036536/149‑1 (MQ=255)
                  tCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGc                                                                                                                                >  2:1401380/1‑149 (MQ=255)
                  tCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGc                                                                                                                                <  2:2559213/149‑1 (MQ=255)
                   cGTTAAAATCCGATTGGTCTGTCGGAGTTTACCTTGTCTCGCCGCGATCCCTGTTATTTTAGCAGAAGTAAAGCGACTGCTGAGATGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCa                                                                                                                               <  2:3299594/149‑1 (MQ=255)
                     ttCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGa                                                                                                                             <  2:2496773/149‑1 (MQ=255)
                     ttCAAATCCGATTGGTCTGGCGGAGTTTACCCGGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGa                                                                                                                             <  2:3479026/149‑1 (MQ=255)
                      taaaaTCCGTTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGTTTTTGTTAGCTGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGAt                                                                                                                            <  1:4537203/147‑1 (MQ=255)
                      tCAAATCCTATTGTTCTGGCGTATTTTACCCTGTCTCGCCGCGATCCCGGTTTTGTTAGCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGAt                                                                                                                            <  1:2195918/149‑1 (MQ=255)
                      tCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGAt                                                                                                                            >  2:1990817/1‑149 (MQ=255)
                         aaTCCTTTTGGTGTGGCTGTGTTTACCGTGTCGCGCCGCGATCCGGTTTTTTTTCTCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGc                                                                                                                         <  1:3528483/149‑1 (MQ=255)
                             cGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCGGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                                     <  2:3589074/149‑1 (MQ=255)
                               aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                                   >  2:985530/1‑149 (MQ=255)
                                ttGGTCTGGCGGTTTTTACCCTGTCTCGCCGCGAGCCGGGTTTTGTTTGCTTAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                                  <  1:3318885/149‑1 (MQ=255)
                                 tGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAgg                                                                                                                 >  2:3958120/1‑149 (MQ=255)
                                      tGGCGGTGTTTACCCTGTCTCGCCGCGATCCGGGTTTTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                            <  1:1770961/149‑1 (MQ=255)
                                       ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCata                                                                                                           >  2:2791659/1‑149 (MQ=255)
                                        gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                                                                          <  2:3541853/149‑1 (MQ=255)
                                          ggATTTTACCCTGTCTCGCCGCGATCCCGGTTTTGTTAGCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTg                                                                                                        <  1:741728/149‑1 (MQ=255)
                                           gAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGa                                                                                                       >  2:613402/1‑149 (MQ=255)
                                              tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    >  2:3713170/1‑149 (MQ=255)
                                                tACCCTGTCTCGCCGCGATCCCGGTTTTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                                                                                  <  1:714955/149‑1 (MQ=255)
                                                tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                                                                                  <  2:3653745/149‑1 (MQ=255)
                                                 aCCCTGTCTCGCCGCGATCCCGGTTTTGTTAGCAGACTTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                                                                                 <  1:1118075/149‑1 (MQ=255)
                                                   ccTGTCTCGCCGCGATCCCGTTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                               <  1:257098/149‑1 (MQ=255)
                                                    cTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCa                                                                                              >  2:498124/1‑149 (MQ=255)
                                                    cTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCa                                                                                              >  2:2270427/1‑149 (MQ=255)
                                                       tctcGCCGCGATCCCGGTTATGTTAGCAGTATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                           <  2:162576/149‑1 (MQ=255)
                                                       tctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                           <  1:766314/149‑1 (MQ=255)
                                                         tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTg                                                                                         >  2:3848065/1‑149 (MQ=255)
                                                         tcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTg                                                                                         <  2:4126607/149‑1 (MQ=255)
                                                           gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaa                                                                                       >  1:2357174/1‑149 (MQ=255)
                                                            ccgcGATCCCGGTTTTGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaaa                                                                                      <  1:320441/149‑1 (MQ=255)
                                                              gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                    >  1:1908210/1‑149 (MQ=255)
                                                              gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                    >  1:3665161/1‑149 (MQ=255)
                                                              gcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                    <  2:2190491/149‑1 (MQ=255)
                                                                gagcccGGTTTTTTTAGCATAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                                                                                  <  1:2866741/146‑1 (MQ=255)
                                                                 aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                                                                                 >  2:885057/1‑149 (MQ=255)
                                                                  tCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgca                                                                                >  2:2645405/1‑149 (MQ=255)
                                                                    ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                              >  1:836967/1‑149 (MQ=255)
                                                                      ggTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                            >  1:1205595/1‑149 (MQ=255)
                                                                          aTGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                        <  2:3198069/149‑1 (MQ=255)
                                                                          aTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                        <  2:861435/149‑1 (MQ=255)
                                                                          aTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTGTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                        >  1:3113353/1‑149 (MQ=255)
                                                                           tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCt                                                                       >  2:2199179/1‑149 (MQ=255)
                                                                            gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                                      <  2:3026429/149‑1 (MQ=255)
                                                                             ttAGCATAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTAt                                                                     <  1:953646/149‑1 (MQ=255)
                                                                                gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                                  >  1:4681123/1‑149 (MQ=255)
                                                                                 cAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                                 <  1:4063949/149‑1 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                                 <  1:3470133/149‑1 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                                 <  1:1387369/149‑1 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                                 >  1:4204494/1‑149 (MQ=255)
                                                                                  aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                                >  2:2234611/1‑149 (MQ=255)
                                                                                  aGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                                >  2:3455336/1‑149 (MQ=255)
                                                                                     aGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                             >  2:2977870/1‑149 (MQ=255)
                                                                                       tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                           <  2:3378898/149‑1 (MQ=255)
                                                                                       tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                           <  1:528687/149‑1 (MQ=255)
                                                                                        aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc                                                          <  1:4577824/149‑1 (MQ=255)
                                                                                        aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc                                                          <  2:4219454/149‑1 (MQ=255)
                                                                                         acgCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                         <  1:414317/147‑1 (MQ=255)
                                                                                         aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                         >  2:4603680/1‑149 (MQ=255)
                                                                                          aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAg                                                        <  2:820735/149‑1 (MQ=255)
                                                                                             gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCg                                                     <  1:2259407/149‑1 (MQ=255)
                                                                                              aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                                    <  2:2091085/149‑1 (MQ=255)
                                                                                                tgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATg                                                  >  2:4040339/1‑149 (MQ=255)
                                                                                                tgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATg                                                  <  1:1025214/149‑1 (MQ=255)
                                                                                                    gAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTc                                              <  1:3779697/149‑1 (MQ=255)
                                                                                                         ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgcgc                                         >  2:1801672/1‑149 (MQ=255)
                                                                                                           aaaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATGAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGt                                       <  1:884726/149‑1 (MQ=255)
                                                                                                            aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTg                                      >  1:3475842/1‑149 (MQ=255)
                                                                                                              aTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                                    <  1:1538563/149‑1 (MQ=255)
                                                                                                               tCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAAc                                   <  1:4112778/149‑1 (MQ=255)
                                                                                                                cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACa                                  <  2:2624696/149‑1 (MQ=255)
                                                                                                                cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACa                                  <  1:3124202/149‑1 (MQ=255)
                                                                                                                cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACa                                  <  2:2152368/149‑1 (MQ=255)
                                                                                                                 aGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                                 >  2:1120095/1‑149 (MQ=255)
                                                                                                                  gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                                >  1:3946354/1‑149 (MQ=255)
                                                                                                                    gTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcg                              <  1:1836071/149‑1 (MQ=255)
                                                                                                                     tCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcgg                             >  1:924146/1‑149 (MQ=255)
                                                                                                                       tGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                           <  2:3115536/149‑1 (MQ=255)
                                                                                                                       tGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                           <  2:2875388/149‑1 (MQ=255)
                                                                                                                        ggCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGa                          <  2:2786931/149‑1 (MQ=255)
                                                                                                                              gAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGcc                    <  1:3011010/149‑1 (MQ=255)
                                                                                                                               aaTGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa                   >  1:3859439/1‑149 (MQ=255)
                                                                                                                                aTGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAg                  >  2:12390/1‑149 (MQ=255)
                                                                                                                                    cAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtg              <  2:2438346/149‑1 (MQ=255)
                                                                                                                                       cGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa           <  1:2403626/149‑1 (MQ=255)
                                                                                                                                       cGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa           <  1:2714395/149‑1 (MQ=255)
                                                                                                                                         aGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAtt         >  2:772062/1‑149 (MQ=255)
                                                                                                                                               cAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGt   <  2:450789/149‑1 (MQ=255)
                                                                                                                                                agagGTGTTTGCGCGCATTTAGCAGATTGCTTGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTTTGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc  <  1:3000688/149‑1 (MQ=255)
                                                                                                                                                agagGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc  <  1:1707199/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTC  >  CP009273/800632‑800924

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