Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,054 G→T 89.7% S222Y (TCT→TAT)  sthA ← pyridine nucleotide transhydrogenase, soluble

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,150,0540GT89.7% 463.1 / 50.3 184S222Y (TCT→TAT) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (15/4);  new base T (75/89);  total (91/93)
Fisher's exact test for biased strand distribution p-value = 7.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

Only 100 of 195 total aligned reads displayed.

ACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGT  >  CP009273/4149909‑4150198
                                                                                                                                                 |                                                                                                                                                
aCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATaa                                                                                                                                                <  1:2109809/148‑1 (MQ=255)
 ccGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAAt                                                                                                                                               >  2:972718/1‑148 (MQ=255)
 ccGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAAt                                                                                                                                               >  1:1245300/1‑148 (MQ=255)
     tGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATACTCTGa                                                                                                                                           >  2:1798291/1‑148 (MQ=255)
      ggCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGAc                                                                                                                                          <  1:607371/148‑1 (MQ=255)
       gCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACa                                                                                                                                         <  1:783132/148‑1 (MQ=255)
        cATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAt                                                                                                                                        >  1:3265212/1‑148 (MQ=255)
           agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAtctc                                                                                                                                     <  2:3031274/148‑1 (MQ=255)
               caggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGa                                                                                                                                 >  2:1642430/1‑148 (MQ=255)
               caggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGa                                                                                                                                 <  2:408764/148‑1 (MQ=255)
                 ggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATc                                                                                                                               >  2:2852419/1‑148 (MQ=255)
                 ggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATc                                                                                                                               <  2:2938594/148‑1 (MQ=255)
                  gcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCg                                                                                                                              >  1:2972479/1‑148 (MQ=255)
                     gTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCgaga                                                                                                                           <  2:1537485/148‑1 (MQ=255)
                       cAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaaa                                                                                                                         >  1:934396/1‑148 (MQ=255)
                         gCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATg                                                                                                                       >  1:2277796/1‑148 (MQ=255)
                          cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                      <  1:2249970/148‑1 (MQ=255)
                           tttCAGTTTTTTACCCGACTACAGATGCATGATCACACCATCCTCACAGCCTTCGAACTTCTCGTACTCTTCCTTGTGACGAATCACTACGCCACTGTTCCAGAGGTGAGCGGAGAGCTAATCTCACATCTCTTGATCGAGACCTGcc                                                                                                                     >  2:1678958/1‑148 (MQ=255)
                            ttCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCca                                                                                                                    >  2:1122069/1‑148 (MQ=255)
                             tCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCcag                                                                                                                   <  1:1609418/148‑1 (MQ=255)
                                  tttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGc                                                                                                              >  1:816765/1‑148 (MQ=255)
                                    tttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGa                                                                                                            <  2:586021/148‑1 (MQ=255)
                                         ccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgcg                                                                                                       <  1:3296948/148‑1 (MQ=255)
                                          cGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                                      <  1:1549002/148‑1 (MQ=255)
                                          cGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                                      >  1:3275095/1‑148 (MQ=255)
                                             cTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtg                                                                                                   <  2:3022168/148‑1 (MQ=255)
                                              ttCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtgt                                                                                                  >  1:685393/1‑148 (MQ=255)
                                               tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGtt                                                                                                 >  2:1940959/1‑148 (MQ=255)
                                                cAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTg                                                                                                <  2:2153365/148‑1 (MQ=255)
                                                 aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGATAAATGCCAGCAGGCGATCGCGGGTGTTGa                                                                                               <  2:3158554/148‑1 (MQ=255)
                                                  gATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATTCCAGCAGGCGATCGCGGGTGTTGAt                                                                                              <  2:1389393/148‑1 (MQ=255)
                                                   aTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                             >  1:2372002/1‑148 (MQ=255)
                                                   aTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                             <  2:1959798/148‑1 (MQ=255)
                                                     gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAg                                                                                           >  1:2085168/1‑148 (MQ=255)
                                                      cATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                                          <  2:1210038/148‑1 (MQ=255)
                                                        tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATc                                                                                        >  1:1525463/1‑148 (MQ=255)
                                                        tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATc                                                                                        <  2:2140922/148‑1 (MQ=255)
                                                        tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATc                                                                                        >  1:2743271/1‑148 (MQ=255)
                                                        tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATc                                                                                        <  2:476017/148‑1 (MQ=255)
                                                         gATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATcc                                                                                       <  2:2378202/148‑1 (MQ=255)
                                                          aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                      >  2:1869754/1‑148 (MQ=255)
                                                           tCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCAc                                                                                     <  1:1847861/148‑1 (MQ=255)
                                                           tCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCAc                                                                                     <  2:2155765/148‑1 (MQ=255)
                                                            cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                    >  2:3236193/1‑148 (MQ=255)
                                                            cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                    <  1:1161790/148‑1 (MQ=255)
                                                             acacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtt                                                                                   <  2:1025517/148‑1 (MQ=255)
                                                              cacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGCTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACttt                                                                                  >  1:2350932/1‑148 (MQ=255)
                                                               acCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                 <  2:1655201/148‑1 (MQ=255)
                                                               acCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                 <  2:1655191/148‑1 (MQ=255)
                                                                ccATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                                                >  2:839940/1‑148 (MQ=255)
                                                                    cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                            <  1:642248/148‑1 (MQ=255)
                                                                    cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                            >  1:1127989/1‑148 (MQ=255)
                                                                     gTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATcc                                                                           <  1:797485/148‑1 (MQ=255)
                                                                       cacaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCAt                                                                         >  2:1696441/1‑148 (MQ=255)
                                                                         caGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                                       <  1:1539753/148‑1 (MQ=255)
                                                                           gCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcg                                                                     <  1:102030/148‑1 (MQ=255)
                                                                            ccTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcgc                                                                    >  1:3258342/1‑148 (MQ=255)
                                                                            ccTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcgc                                                                    >  1:2271244/1‑148 (MQ=255)
                                                                              ttCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCgg                                                                  >  2:3012341/1‑148 (MQ=255)
                                                                              ttCGATCTTCTCGTACTCTTCGATGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCCGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCgg                                                                  >  2:3185783/1‑148 (MQ=255)
                                                                                 gATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAg                                                               <  1:222358/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              <  2:1320364/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              <  2:2249057/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              <  1:2766765/148‑1 (MQ=255)
                                                                                   tcttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGAt                                                             <  2:334669/148‑1 (MQ=255)
                                                                                    cttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGACGATc                                                            >  1:1758314/1‑148 (MQ=255)
                                                                                    cttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATc                                                            >  2:2071474/1‑148 (MQ=255)
                                                                                     ttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCg                                                           >  2:2526125/1‑148 (MQ=255)
                                                                                      tctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGa                                                          <  2:2415501/148‑1 (MQ=255)
                                                                                         cGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                                       >  1:2487003/1‑148 (MQ=255)
                                                                                            aCTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCata                                                    >  1:2877495/1‑148 (MQ=255)
                                                                                             ctctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCatat                                                   >  2:1816782/1‑148 (MQ=255)
                                                                                                   gTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                                             <  1:3172426/148‑1 (MQ=255)
                                                                                                       tgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGa                                         <  2:362295/148‑1 (MQ=255)
                                                                                                        gACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGAt                                        <  2:1299994/148‑1 (MQ=255)
                                                                                                        gACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTAACAGCCGAt                                        >  2:1190732/1‑148 (MQ=255)
                                                                                                         aCGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                                       >  1:903794/1‑148 (MQ=255)
                                                                                                         aCGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                                       >  1:1277927/1‑148 (MQ=255)
                                                                                                             aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                   >  1:2661784/1‑148 (MQ=255)
                                                                                                             aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                   >  2:2209144/1‑148 (MQ=255)
                                                                                                                 cTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGc                               >  2:2684457/1‑148 (MQ=255)
                                                                                                                   aCGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATAGTCACAGCCGATCACTCCAGCAc                             >  2:1112186/1‑148 (MQ=255)
                                                                                                                       cACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATa                         >  2:3112459/1‑148 (MQ=255)
                                                                                                                       cACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATa                         >  2:3112386/1‑148 (MQ=255)
                                                                                                                            ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATaa                    <  1:617985/148‑1 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               >  2:705134/1‑148 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               >  1:2303099/1‑148 (MQ=255)
                                                                                                                                  aaGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACAt              <  1:1504850/148‑1 (MQ=255)
                                                                                                                                  aaGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACAt              >  2:2280351/1‑148 (MQ=255)
                                                                                                                                  aaGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACAt              >  1:507076/1‑148 (MQ=255)
                                                                                                                                        tAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGcgcg        <  2:2736188/148‑1 (MQ=255)
                                                                                                                                         aGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCgg       <  1:2858592/148‑1 (MQ=255)
                                                                                                                                          ggAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGt      <  1:1248224/148‑1 (MQ=255)
                                                                                                                                          ggAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGt      <  2:2331895/148‑1 (MQ=255)
                                                                                                                                          ggAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGt      <  1:839940/148‑1 (MQ=255)
                                                                                                                                             gagataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCg   >  1:1002500/1‑148 (MQ=255)
                                                                                                                                              agataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgt  <  1:1887182/148‑1 (MQ=255)
                                                                                                                                              agataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgt  >  2:1460102/1‑148 (MQ=255)
                                                                                                                                              agagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgt  >  1:1920305/1‑148 (MQ=255)
                                                                                                                                              agagaATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgt  >  1:2626003/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGT  >  CP009273/4149909‑4150198

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: