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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273800,7790GT99.7% 1003.8 / ‑4.7 300E607D (GAG→GAT) ybhJputative hydratase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (158/141);  minor base C (0/1);  total (158/142)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 302 total aligned reads displayed.

GGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTC  >  CP009273/800633‑800924
                                                                                                                                                  |                                                                                                                                                 
ggTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATgt                                                                                                                                                 >  2:4534607/1‑149 (MQ=255)
   gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGATTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGtt                                                                                                                                              <  2:4417443/149‑1 (MQ=255)
   gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGtt                                                                                                                                              >  2:3497435/1‑149 (MQ=255)
   gAAACCTCTTCTTATCGTTCAAATCCGATGGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGtt                                                                                                                                              <  2:4584526/149‑1 (MQ=255)
   gAAACCTCTTCTTATCGTACAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGtt                                                                                                                                              <  1:4285918/149‑1 (MQ=255)
    aaaCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGttt                                                                                                                                             >  2:3407884/1‑149 (MQ=255)
      aCCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgc                                                                                                                                           >  1:1435117/1‑149 (MQ=255)
        ctctTGTTATCGTTAAAATCCGTTTGGTCTGGCGGATTTTACCCTGTCTCGCCGCGATCCCGTTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTgcgc                                                                                                                                         <  2:1619004/149‑1 (MQ=255)
           tttttATCTTTAAAATCCGTTTGGTGTGGGGTTTTTTACCCTGTCGCGCGGCGATCCCGTTTTTTTTAGCACAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCa                                                                                                                                      <  1:463711/146‑1 (MQ=255)
           ttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCa                                                                                                                                      >  1:3359371/1‑149 (MQ=255)
           ttcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCa                                                                                                                                      >  2:2504922/1‑149 (MQ=255)
            tcttATCGTTCAAATCCGTTTGGTCTGGGGGTGTTTACCCTGTCTCGCCGCGATCCGGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                     <  1:3075521/149‑1 (MQ=255)
            tcttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                     >  2:528216/1‑149 (MQ=255)
            acttAACGTAAAAATCCTTTTGGTGTGGCGTAGTTTACCCTGTCTCGCCGCTAACCCGTTTATTTTAGCAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAt                                                                                                                                     <  1:1847578/148‑1 (MQ=255)
               tATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTaa                                                                                                                                  >  2:1722076/1‑149 (MQ=255)
                 tCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGc                                                                                                                                >  2:4659537/1‑149 (MQ=255)
                   gTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCTATCCCGTTTATGTTAGAAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAg                                                                                                                              <  1:1171771/149‑1 (MQ=255)
                   gTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAg                                                                                                                              >  1:2384627/1‑149 (MQ=255)
                      cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGTTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGAtt                                                                                                                           <  1:1200581/149‑1 (MQ=255)
                       aattCCGATTGGTCTGGGTGTGTTTACCGTGTCTCGCGGCGTTCCCGTTTTTGTTAGCAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTg                                                                                                                          <  1:3873069/146‑1 (MQ=255)
                        aaTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGc                                                                                                                         >  1:3684762/1‑149 (MQ=255)
                         aTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCt                                                                                                                        >  2:1111647/1‑149 (MQ=255)
                           ccGATTGGTCTGGCGGTGTTTACCCTGTCTCGCCGCTATCCCGTTTATTTTAGCATATGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTgg                                                                                                                      <  1:1837101/149‑1 (MQ=255)
                            cGATTGGTCTGGGGGAGTTTACCCTGTCTCGCCGCGATCCGGGTTTTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                                     <  1:312257/149‑1 (MQ=255)
                            cGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                                     >  1:3934658/1‑149 (MQ=255)
                             gATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTc                                                                                                                    >  2:3852069/1‑149 (MQ=255)
                               ttGGTCGGGGGGTGTTTACCCTGTCTCGCCGCGATCCCGGTTTTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                                  <  2:2848793/149‑1 (MQ=255)
                                 ggTGTGGGGGTGTTTACCCTGTCTCGCCGCGATCCGGTTTATTTTAGCATATGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGa                                                                                                                <  1:4623572/149‑1 (MQ=255)
                                  gtggggCGGTGTTTACCCGGCCTCGCCGCTAGCCGGTTTTTGTTAGCACAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                               <  1:4408582/145‑1 (MQ=255)
                                       gCGGAGTTTACCCGGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                                                                          <  1:4366515/149‑1 (MQ=255)
                                         ggAGTTTACCCTGTCTCGCCGCGATCCGGTTTTTGTTAGCAGATGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTg                                                                                                        <  1:4709432/149‑1 (MQ=255)
                                         ggAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTg                                                                                                        >  1:3462544/1‑149 (MQ=255)
                                           aGTTTCCCCTGTCCCGCCGCGATCCGGTTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                      <  1:148925/149‑1 (MQ=255)
                                            ttttCCCCTGTCTCGCCGCTATCCGGTTTATGTTAGCTTAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                                     <  1:251119/148‑1 (MQ=255)
                                            ttttACCCTGACTCGCCGCGATCCCGTTTATGTTAGCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                                     <  1:2273736/148‑1 (MQ=255)
                                             tttACCGTGTCTCGCCGCGAGCCGGGTTTTTTTAGCAGAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    <  1:4546633/149‑1 (MQ=255)
                                             tttACCCTGTCTCGCCGCGTGCCCGGTTATGTTAGCAGAAGAAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    <  1:2154853/149‑1 (MQ=255)
                                             tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCATAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    <  2:3776232/149‑1 (MQ=255)
                                             tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    >  1:1319641/1‑149 (MQ=255)
                                             tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    >  2:1510804/1‑149 (MQ=255)
                                             tttACCCTGGCTCGCCGCGATCCCGTTTATTTTAGCATAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATcc                                                                                                    <  1:4477609/149‑1 (MQ=255)
                                               tACCCTGTCTCGCCGCGAGCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGc                                                                                                  <  1:2597105/149‑1 (MQ=255)
                                                ccccGGTCGCGCCGCGATCCGGGTTTGGTTAGCAGACGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                                                                                 <  1:1268112/148‑1 (MQ=255)
                                                aCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGACGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCt                                                                                                 <  2:1557029/149‑1 (MQ=255)
                                                  ccTGTCTCGCCGCGATCCGGGTTTTTTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGc                                                                                               <  1:2254702/149‑1 (MQ=255)
                                                    tGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCaa                                                                                             >  1:1049649/1‑149 (MQ=255)
                                                      tctcGCCGCTTTCCCGGTTTTGTTAGCATAATTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAAc                                                                                           <  1:1584484/149‑1 (MQ=255)
                                                          gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaa                                                                                       <  1:1722076/149‑1 (MQ=255)
                                                          gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaa                                                                                       >  2:1295387/1‑149 (MQ=255)
                                                              cgATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTg                                                                                   <  2:4718541/149‑1 (MQ=255)
                                                               gATCCCGGTTATGTTAGCAGAAGTAGAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                                                                                  <  2:1258043/149‑1 (MQ=255)
                                                               gATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                                                                                  <  2:36106/149‑1 (MQ=255)
                                                                aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                                                                                 >  1:269979/1‑149 (MQ=255)
                                                                   ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                              >  2:3006451/1‑149 (MQ=255)
                                                                   ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                              >  1:2145263/1‑149 (MQ=255)
                                                                    cGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagcc                                                                             >  2:2424056/1‑148 (MQ=255)
                                                                    cGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagca                                                                             >  2:1165914/1‑149 (MQ=255)
                                                                     ggTTATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                            <  1:958418/149‑1 (MQ=255)
                                                                      gTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATg                                                                           >  2:92169/1‑149 (MQ=255)
                                                                         aTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                        >  1:2994242/1‑149 (MQ=255)
                                                                         aTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTc                                                                        <  2:2244508/149‑1 (MQ=255)
                                                                          tGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCt                                                                       >  1:581103/1‑149 (MQ=255)
                                                                              aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGc                                                                   >  2:4247566/1‑149 (MQ=255)
                                                                               gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                                  >  2:704394/1‑149 (MQ=255)
                                                                                  gAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                               >  1:799577/1‑149 (MQ=255)
                                                                                   aaGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa                                                              <  2:1909652/149‑1 (MQ=255)
                                                                                    aGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                             <  2:1681651/149‑1 (MQ=255)
                                                                                    aGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                             >  1:3643318/1‑149 (MQ=255)
                                                                                     ttAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                            <  1:877828/148‑1 (MQ=255)
                                                                                     gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                            >  2:3477562/1‑149 (MQ=255)
                                                                                     gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                            >  2:4135059/1‑149 (MQ=255)
                                                                                     gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaaa                                                            >  2:1955979/1‑149 (MQ=255)
                                                                                      tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTGTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                           <  1:4783273/149‑1 (MQ=255)
                                                                                       aaaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc                                                          >  1:2055377/1‑149 (MQ=255)
                                                                                        aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                         >  1:694263/1‑149 (MQ=255)
                                                                                          gCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAgg                                                       <  1:1020712/149‑1 (MQ=255)
                                                                                          gCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAgg                                                       >  1:3284228/1‑149 (MQ=255)
                                                                                             aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                                    <  1:3832989/149‑1 (MQ=255)
                                                                                              ctgctgAGCTGGAAAATAAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGAt                                                   >  2:3257400/1‑149 (MQ=255)
                                                                                                 ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTGTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGt                                                <  1:2297129/149‑1 (MQ=255)
                                                                                                  tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                               >  2:3729271/1‑149 (MQ=255)
                                                                                                  tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCAGTAAGCAGAGTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                               <  2:578433/149‑1 (MQ=255)
                                                                                                     gCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTg                                            >  1:4784842/1‑149 (MQ=255)
                                                                                                       tGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgcg                                          >  1:4603791/1‑149 (MQ=255)
                                                                                                        ggAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTgcgc                                         <  1:3047176/149‑1 (MQ=255)
                                                                                                           aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTg                                      >  2:50605/1‑149 (MQ=255)
                                                                                                             aTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                                    <  2:3625754/149‑1 (MQ=255)
                                                                                                               cAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACa                                  >  2:2129584/1‑149 (MQ=255)
                                                                                                                   gTCTGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAggcg                              >  2:4112045/1‑149 (MQ=255)
                                                                                                                      tGGCGGGGAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                           >  2:3690212/1‑149 (MQ=255)
                                                                                                                           gggAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTg                      >  1:791817/1‑149 (MQ=255)
                                                                                                                            ggAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc                     >  1:4518299/1‑149 (MQ=255)
                                                                                                                             gAATGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGcc                    <  2:3024318/149‑1 (MQ=255)
                                                                                                                              aaTGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCa                   >  2:1636810/1‑149 (MQ=255)
                                                                                                                                tGTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGc                 >  2:2703463/1‑149 (MQ=255)
                                                                                                                                 gTCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCg                <  1:3039115/149‑1 (MQ=255)
                                                                                                                                  tCAGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgt               <  1:4013048/149‑1 (MQ=255)
                                                                                                                                    aGCGAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCgtgt             >  2:2021793/1‑149 (MQ=255)
                                                                                                                                      cGAGCTGACAGATGTGTTTGCGCGCATTACGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGa           <  2:1860045/149‑1 (MQ=255)
                                                                                                                                       gAGCTGACAGATGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt          >  1:4213842/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGGCGGTC  >  CP009273/800633‑800924

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