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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273800,7790GT93.8% 823.6 / 42.4 276E607D (GAG→GAT) ybhJputative hydratase
Reads supporting (aligned to +/- strand):  ref base G (11/6);  new base T (128/129);  total (140/136)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 286 total aligned reads displayed.

CGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGG  >  CP009273/800632‑800927
                                                                                                                                                   |                                                                                                                                                    
cGGTGAAACCTCTTCTTATCGTTCAAATCCGATGGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATgt                                                                                                                                                    <  2:2444571/150‑1 (MQ=255)
 ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAtgtg                                                                                                                                                   <  1:2702472/150‑1 (MQ=255)
 ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAtgtg                                                                                                                                                   <  1:2702269/150‑1 (MQ=255)
    gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGttt                                                                                                                                                >  2:1072298/1‑150 (MQ=255)
     aaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTg                                                                                                                                               >  1:955703/1‑150 (MQ=255)
     aaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTg                                                                                                                                               >  1:3232016/1‑150 (MQ=255)
     aaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTg                                                                                                                                               <  2:2853222/150‑1 (MQ=255)
         ctctTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgcg                                                                                                                                           >  2:3450808/1‑150 (MQ=255)
           cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCa                                                                                                                                         >  1:2829968/1‑150 (MQ=255)
           cttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCa                                                                                                                                         >  1:1346194/1‑150 (MQ=255)
            ttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                        >  1:916233/1‑150 (MQ=255)
            ttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                        >  1:4322460/1‑150 (MQ=255)
            ttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                        <  1:2431260/150‑1 (MQ=255)
            ttcttATCGTTAAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                        <  1:1926239/150‑1 (MQ=255)
             tcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAtt                                                                                                                                       <  1:4240167/150‑1 (MQ=255)
                tATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAg                                                                                                                                    >  1:3540801/1‑150 (MQ=255)
                 aTCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGc                                                                                                                                   >  1:2113194/1‑150 (MQ=255)
                   cGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAg                                                                                                                                 >  1:3840128/1‑150 (MQ=255)
                   cGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAg                                                                                                                                 >  1:3093802/1‑150 (MQ=255)
                   cGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAg                                                                                                                                 >  2:3402615/1‑150 (MQ=255)
                     ttCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGAt                                                                                                                               >  2:2198408/1‑150 (MQ=255)
                        aaaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGc                                                                                                                            >  2:295626/1‑150 (MQ=255)
                         aaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCt                                                                                                                           <  2:2250862/150‑1 (MQ=255)
                          ctCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTg                                                                                                                          <  2:1592773/149‑1 (MQ=255)
                          aTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTg                                                                                                                          <  2:1592790/150‑1 (MQ=255)
                          aTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTg                                                                                                                          <  1:2993290/150‑1 (MQ=255)
                            ccGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                                        >  2:1894479/1‑150 (MQ=255)
                               aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                                     <  1:501381/150‑1 (MQ=255)
                                  ggTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                                  >  1:4033267/1‑150 (MQ=255)
                                     cgggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                               <  2:3860043/148‑1 (MQ=255)
                                     cTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                               >  1:3620216/1‑150 (MQ=255)
                                     cTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                               <  1:2829042/150‑1 (MQ=255)
                                     cTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                               <  2:272392/150‑1 (MQ=255)
                                      tGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCata                                                                                                              >  2:530846/1‑150 (MQ=255)
                                      tGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCata                                                                                                              <  2:2241450/150‑1 (MQ=255)
                                       ggCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTAGGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                                                                             >  1:636466/1‑150 (MQ=255)
                                        gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATAtt                                                                                                            <  2:4277851/150‑1 (MQ=255)
                                         cGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTg                                                                                                           <  2:4423330/150‑1 (MQ=255)
                                          ggAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGa                                                                                                          >  2:2974398/1‑150 (MQ=255)
                                          ggAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGa                                                                                                          <  2:2268696/150‑1 (MQ=255)
                                           gAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                         <  1:3982367/150‑1 (MQ=255)
                                                 aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTg                                                                                                   <  1:202396/150‑1 (MQ=255)
                                                  cccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                                  >  2:4000334/1‑150 (MQ=255)
                                                   ccTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCa                                                                                                 <  1:4476041/150‑1 (MQ=255)
                                                       tctcGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACt                                                                                             <  1:792917/150‑1 (MQ=255)
                                                           gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaaa                                                                                         <  2:2882327/150‑1 (MQ=255)
                                                            ccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAt                                                                                        <  1:3881578/150‑1 (MQ=255)
                                                             cgcgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAtt                                                                                       >  2:4263571/1‑150 (MQ=255)
                                                               cgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                                                                                     <  2:3287944/150‑1 (MQ=255)
                                                                 aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgca                                                                                   <  1:4145072/150‑1 (MQ=255)
                                                                 aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgca                                                                                   >  2:1194326/1‑150 (MQ=255)
                                                                       gTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATgg                                                                             <  1:1681351/150‑1 (MQ=255)
                                                                       gTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATgg                                                                             <  1:2691494/150‑1 (MQ=255)
                                                                        ttATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGt                                                                            >  1:3719872/1‑150 (MQ=255)
                                                                         tATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                           >  1:1606016/1‑150 (MQ=255)
                                                                         tATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                           <  2:3327921/150‑1 (MQ=255)
                                                                         aaTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                           <  2:4396707/149‑1 (MQ=255)
                                                                           tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                                         >  2:3188845/1‑150 (MQ=255)
                                                                            gTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTAt                                                                        <  2:1973614/150‑1 (MQ=255)
                                                                             ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATg                                                                       >  2:3368920/1‑150 (MQ=255)
                                                                              tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc                                                                      >  2:687661/1‑150 (MQ=255)
                                                                              tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc                                                                      <  2:3455131/150‑1 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                                   >  2:364142/1‑150 (MQ=255)
                                                                                 cAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGt                                                                   >  2:3598274/1‑150 (MQ=255)
                                                                                    aaGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                                <  1:3602284/150‑1 (MQ=255)
                                                                                        aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAAGGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                            <  2:148232/150‑1 (MQ=255)
                                                                                          aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAgg                                                          >  2:1751213/1‑150 (MQ=255)
                                                                                          aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAgg                                                          >  1:1200132/1‑150 (MQ=255)
                                                                                             gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                                       >  2:732790/1‑150 (MQ=255)
                                                                                             gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                                       >  1:2236328/1‑150 (MQ=255)
                                                                                              aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGAt                                                      >  2:475000/1‑150 (MQ=255)
                                                                                              aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGAt                                                      >  2:3050530/1‑150 (MQ=255)
                                                                                                tgctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATgg                                                    >  1:1349081/1‑150 (MQ=255)
                                                                                                 gctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGt                                                   <  2:4174640/150‑1 (MQ=255)
                                                                                                  ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                                  >  2:2949230/1‑150 (MQ=255)
                                                                                                   tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTc                                                 <  1:1330652/150‑1 (MQ=255)
                                                                                                   tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTc                                                 <  2:102198/150‑1 (MQ=255)
                                                                                                      gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgc                                              <  1:1151908/150‑1 (MQ=255)
                                                                                                           aaaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTg                                         <  1:1224975/150‑1 (MQ=255)
                                                                                                            aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGa                                        >  1:3675204/1‑150 (MQ=255)
                                                                                                            aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGa                                        >  1:1408772/1‑150 (MQ=255)
                                                                                                             aaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                                       >  1:1287235/1‑150 (MQ=255)
                                                                                                             aaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                                       >  1:1885995/1‑150 (MQ=255)
                                                                                                              aTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAAc                                      <  1:692388/150‑1 (MQ=255)
                                                                                                                 aGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                                   >  1:1023950/1‑150 (MQ=255)
                                                                                                                       tGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGa                             >  2:1771461/1‑150 (MQ=255)
                                                                                                                          cGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCt                          <  1:3450808/150‑1 (MQ=255)
                                                                                                                           ggggAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                         <  1:1072298/150‑1 (MQ=255)
                                                                                                                              gAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa                      <  2:1388403/150‑1 (MQ=255)
                                                                                                                                 tGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                   <  2:820672/150‑1 (MQ=255)
                                                                                                                                         aGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg           <  2:3197474/150‑1 (MQ=255)
                                                                                                                                         aGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg           >  1:61889/1‑150 (MQ=255)
                                                                                                                                         aGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTg           >  1:61729/1‑150 (MQ=255)
                                                                                                                                           cTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGc         <  2:618886/150‑1 (MQ=255)
                                                                                                                                            tGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCg        >  2:3705466/1‑150 (MQ=255)
                                                                                                                                            tGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCg        <  1:3585571/150‑1 (MQ=255)
                                                                                                                                             gACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCgg       >  2:3744557/1‑150 (MQ=255)
                                                                                                                                               cAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc     <  2:3359559/150‑1 (MQ=255)
                                                                                                                                               cAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTc     >  2:1808387/1‑150 (MQ=255)
                                                                                                                                                agatGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCt    >  2:1331447/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTCTGG  >  CP009273/800632‑800927

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: ‑