Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 555,927 G→T 18.5% D259Y (GAT→TAT)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273555,9270GT18.5% 755.4 / 106.3 268D259Y (GAT→TAT) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (94/124);  new base T (27/23);  total (121/147)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 273 total aligned reads displayed.

TGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTA  >  CP009273/555780‑556074
                                                                                                                                                   |                                                                                                                                                   
 ggCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGa                                                                                                                                                    <  2:2026715/148‑1 (MQ=255)
   cGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGata                                                                                                                                                  <  1:3590725/148‑1 (MQ=255)
     tttACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATAtt                                                                                                                                                <  1:632015/148‑1 (MQ=255)
      ttACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATAttt                                                                                                                                               >  1:2636602/1‑148 (MQ=255)
      ttACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATAttt                                                                                                                                               >  1:984205/1‑148 (MQ=255)
       tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATAtttt                                                                                                                                              >  1:3656682/1‑148 (MQ=255)
         cGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCg                                                                                                                                            <  2:1137151/148‑1 (MQ=255)
          gCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGa                                                                                                                                           <  1:3148383/148‑1 (MQ=255)
          gCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGAGATTTTCGa                                                                                                                                           >  2:2815283/1‑148 (MQ=255)
           caacaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGAt                                                                                                                                          <  2:339114/148‑1 (MQ=255)
             acaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATgg                                                                                                                                        >  1:2497718/1‑148 (MQ=255)
              caaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGc                                                                                                                                       >  1:531046/1‑148 (MQ=255)
                 tGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGtt                                                                                                                                    >  1:812070/1‑148 (MQ=255)
                 tGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGtt                                                                                                                                    <  2:651774/148‑1 (MQ=255)
                   gTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTgg                                                                                                                                  <  1:2809258/148‑1 (MQ=255)
                    gtCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGa                                                                                                                                 <  2:2858769/147‑1 (MQ=255)
                          gAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCgt                                                                                                                           >  2:440222/1‑148 (MQ=255)
                               aTTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGgtgt                                                                                                                      <  2:950567/148‑1 (MQ=255)
                               aTTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGgtgt                                                                                                                      >  2:143263/1‑148 (MQ=255)
                                ttttCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTa                                                                                                                     <  1:3373324/148‑1 (MQ=255)
                                ttttCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTa                                                                                                                     >  1:2126605/1‑148 (MQ=255)
                                ttttCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTa                                                                                                                     <  2:1620978/148‑1 (MQ=255)
                                 tttCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTACTGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTAc                                                                                                                    >  2:3463524/1‑148 (MQ=255)
                                 tttCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTAc                                                                                                                    >  2:224495/1‑148 (MQ=255)
                                   tCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACgg                                                                                                                  <  2:10029/148‑1 (MQ=255)
                                    cgcgGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGc                                                                                                                 <  2:3167827/148‑1 (MQ=255)
                                     gcgGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCt                                                                                                                <  2:503359/148‑1 (MQ=255)
                                     gcgGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCt                                                                                                                <  1:3024967/148‑1 (MQ=255)
                                        gAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTa                                                                                                             <  1:440222/148‑1 (MQ=255)
                                           aTGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTAttc                                                                                                          >  2:1406676/1‑148 (MQ=255)
                                             ggATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTAttctt                                                                                                        >  1:2841678/1‑148 (MQ=255)
                                               atatCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATtcttct                                                                                                      <  2:110258/148‑1 (MQ=255)
                                               atatCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATtcttct                                                                                                      >  1:1000994/1‑148 (MQ=255)
                                                 atCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGa                                                                                                    >  2:2417161/1‑148 (MQ=255)
                                                   cATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGata                                                                                                  >  2:850900/1‑148 (MQ=255)
                                                    aTTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGataa                                                                                                 >  2:3125106/1‑148 (MQ=255)
                                                     ttCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGataat                                                                                                >  2:2514663/1‑148 (MQ=255)
                                                     ttCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGataat                                                                                                <  1:140194/148‑1 (MQ=255)
                                                         gAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTCCGGCTTTATTCTTCTGATAATATGt                                                                                            >  2:2351955/1‑148 (MQ=255)
                                                         gAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGt                                                                                            <  1:1270479/148‑1 (MQ=255)
                                                            cAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATc                                                                                         <  1:3404792/148‑1 (MQ=255)
                                                            cAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATc                                                                                         >  1:1977732/1‑148 (MQ=255)
                                                             aGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATcc                                                                                        <  1:2173399/148‑1 (MQ=255)
                                                             aGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATAAAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATcc                                                                                        >  1:2379035/1‑148 (MQ=255)
                                                                ggAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGa                                                                                     >  2:2567840/1‑148 (MQ=255)
                                                                 gAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGAt                                                                                    <  1:1016602/148‑1 (MQ=255)
                                                                  aataatATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATa                                                                                   >  2:741979/1‑148 (MQ=255)
                                                                   ataataTTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAg                                                                                  <  2:2691258/148‑1 (MQ=255)
                                                                    taataTTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGc                                                                                 >  1:1914453/1‑148 (MQ=255)
                                                                     aataTTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGcc                                                                                >  1:848408/1‑148 (MQ=255)
                                                                          tGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCaa                                                                           <  1:133041/148‑1 (MQ=255)
                                                                            gCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAgg                                                                         >  1:3779683/1‑148 (MQ=255)
                                                                            gCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAgg                                                                         >  1:2155480/1‑148 (MQ=255)
                                                                             cACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAggg                                                                        <  1:1476568/148‑1 (MQ=255)
                                                                             cACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAggg                                                                        <  1:3076431/148‑1 (MQ=255)
                                                                              aCCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGt                                                                       <  1:2227858/148‑1 (MQ=255)
                                                                              aCCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATAATTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGt                                                                       >  2:295184/1‑148 (MQ=255)
                                                                               ccTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGtt                                                                      <  2:1714246/148‑1 (MQ=255)
                                                                                 tGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTg                                                                    <  1:3702138/148‑1 (MQ=255)
                                                                                  gggTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGc                                                                   >  2:69889/1‑148 (MQ=255)
                                                                                   ggTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGcc                                                                  >  2:2196422/1‑148 (MQ=255)
                                                                                    gTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGccc                                                                 <  1:145140/148‑1 (MQ=255)
                                                                                    gTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGccc                                                                 <  2:1273357/148‑1 (MQ=255)
                                                                                    gTACAGCGCGCCACTATTCCGCTGAAAAGTGCACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTATGATAATATGTATCCTGATAGCCAGCAAGGGTTTGccc                                                                 <  1:2073537/148‑1 (MQ=255)
                                                                                       cAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCaa                                                              >  1:1730032/1‑148 (MQ=255)
                                                                                       cAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCaa                                                              <  2:811053/148‑1 (MQ=255)
                                                                                         gcgcgcCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACg                                                            >  1:2240240/1‑148 (MQ=255)
                                                                                          cgcgcCATTATTCCGCTGAAAAGTGAACGGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACgg                                                           >  1:945738/1‑148 (MQ=255)
                                                                                           gcgcCATTAGTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGt                                                          <  2:2001110/148‑1 (MQ=255)
                                                                                              ccATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACg                                                       >  1:2773100/1‑148 (MQ=255)
                                                                                              ccATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACg                                                       <  1:2208126/148‑1 (MQ=255)
                                                                                                attattCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTg                                                     >  2:3106196/1‑148 (MQ=255)
                                                                                                 ttattCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTgg                                                    >  2:3049856/1‑148 (MQ=255)
                                                                                                    ttCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGAt                                                 >  2:2836331/1‑148 (MQ=255)
                                                                                                     tCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGAtt                                                >  1:371784/1‑148 (MQ=255)
                                                                                                         cTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTTATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGccc                                            <  1:2933909/148‑1 (MQ=255)
                                                                                                            aaaaGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTa                                         <  2:1158845/148‑1 (MQ=255)
                                                                                                             aaaGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTAc                                        >  1:2902640/1‑148 (MQ=255)
                                                                                                              aaGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTAcg                                       >  1:183725/1‑148 (MQ=255)
                                                                                                                    aCTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGcc                                 <  2:1487467/148‑1 (MQ=255)
                                                                                                                          aaTGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTa                           <  1:1411116/148‑1 (MQ=255)
                                                                                                                          aaTGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTa                           <  1:3313517/148‑1 (MQ=255)
                                                                                                                           aTGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTaa                          <  2:1933061/148‑1 (MQ=255)
                                                                                                                            tGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAAc                         >  2:1774489/1‑148 (MQ=255)
                                                                                                                              ggAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGa                       <  2:2190296/148‑1 (MQ=255)
                                                                                                                               gagaCGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGAt                      <  2:2794191/148‑1 (MQ=255)
                                                                                                                               gagaCGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGAt                      <  2:1369499/148‑1 (MQ=255)
                                                                                                                                agaCGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGAtt                     <  1:850900/148‑1 (MQ=255)
                                                                                                                                 gaCGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTc                    <  2:1198057/148‑1 (MQ=255)
                                                                                                                                  aCGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCg                   <  2:3776062/148‑1 (MQ=255)
                                                                                                                                     gCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCa                >  2:2139491/1‑148 (MQ=255)
                                                                                                                                       aaTACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCaaa              <  1:143198/148‑1 (MQ=255)
                                                                                                                                         tACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAAt            <  1:252616/148‑1 (MQ=255)
                                                                                                                                             ggAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGtt        >  2:1875846/1‑148 (MQ=255)
                                                                                                                                             ggAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGtt        >  1:1188989/1‑148 (MQ=255)
                                                                                                                                              gAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGttt       >  1:537599/1‑148 (MQ=255)
                                                                                                                                                  tGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatt   <  2:3634557/148‑1 (MQ=255)
                                                                                                                                                  tGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatt   <  2:2604211/148‑1 (MQ=255)
                                                                                                                                                  tGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatt   <  1:3605138/148‑1 (MQ=255)
                                                                                                                                                   gATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatta  >  1:3106979/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTATTCCGCTGAAAAGTGAACTGGTAATGGGAGACGGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTA  >  CP009273/555780‑556074

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: