Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,150,054 G→T 100% 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,0540GT98.3% 764.1 / 2.9 239S222Y (TCT→TAT) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base T (124/111);  total (128/111)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
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 245 total aligned reads displayed.

ACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGT  >  CP009273/4149909‑4150198
                                                                                                                                                 |                                                                                                                                                
aCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATaa                                                                                                                                                <  1:2599631/148‑1 (MQ=255)
 ccGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAAt                                                                                                                                               <  1:1860450/148‑1 (MQ=255)
     tGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGa                                                                                                                                           <  2:3379982/148‑1 (MQ=255)
      ggCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGAc                                                                                                                                          <  2:2223637/148‑1 (MQ=255)
      ggCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGAc                                                                                                                                          <  2:2223572/148‑1 (MQ=255)
       gCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACa                                                                                                                                         >  2:1691187/1‑148 (MQ=255)
       gCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACa                                                                                                                                         <  2:2080832/148‑1 (MQ=255)
        cATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGACGTGAGAGGAGAGATAATCTGACAt                                                                                                                                        >  2:2692765/1‑148 (MQ=255)
        cATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAt                                                                                                                                        >  2:2731625/1‑148 (MQ=255)
        cATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAt                                                                                                                                        <  1:1790551/148‑1 (MQ=255)
           agagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAtctc                                                                                                                                     >  1:2622298/1‑148 (MQ=255)
            gagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATctct                                                                                                                                    >  1:2078695/1‑148 (MQ=255)
             agCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATACTCTGACATCTCtt                                                                                                                                   >  2:1823631/1‑148 (MQ=255)
              gcaggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTg                                                                                                                                  >  2:191117/1‑148 (MQ=255)
               caggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGa                                                                                                                                 >  1:2471912/1‑148 (MQ=255)
               caggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGa                                                                                                                                 <  2:2043207/148‑1 (MQ=255)
                aggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGAt                                                                                                                                <  1:3180088/148‑1 (MQ=255)
                 ggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATc                                                                                                                               >  1:3386157/1‑148 (MQ=255)
                  gcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCg                                                                                                                              <  2:2376562/148‑1 (MQ=255)
                  gcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCg                                                                                                                              <  2:2915988/148‑1 (MQ=255)
                    agTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCgag                                                                                                                            <  1:214806/148‑1 (MQ=255)
                    agTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCgag                                                                                                                            <  1:214297/148‑1 (MQ=255)
                    agTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCgag                                                                                                                            <  2:635128/148‑1 (MQ=255)
                      tCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaa                                                                                                                          >  1:3314869/1‑148 (MQ=255)
                       cAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGGCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaaa                                                                                                                         >  1:3334555/1‑148 (MQ=255)
                         gCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATg                                                                                                                       >  1:1852507/1‑148 (MQ=255)
                          cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                      >  1:1594709/1‑148 (MQ=255)
                          cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                      <  2:2927734/148‑1 (MQ=255)
                          cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                      >  1:3036454/1‑148 (MQ=255)
                              cAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCcagc                                                                                                                  <  1:911614/148‑1 (MQ=255)
                               aGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCcagca                                                                                                                 >  2:938115/1‑148 (MQ=255)
                                gTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCcagcag                                                                                                                >  1:971964/1‑148 (MQ=255)
                                 ttttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCcagcagt                                                                                                               >  2:2480779/1‑147 (MQ=255)
                                 ttttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAgg                                                                                                               <  2:3208637/148‑1 (MQ=255)
                                    tttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGa                                                                                                            <  2:17093/148‑1 (MQ=255)
                                    tttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGa                                                                                                            <  1:3221625/148‑1 (MQ=255)
                                    tttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGa                                                                                                            <  2:3158289/148‑1 (MQ=255)
                                       aCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcg                                                                                                         >  1:2292611/1‑148 (MQ=255)
                                       aCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcg                                                                                                         >  1:499818/1‑148 (MQ=255)
                                       aCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcg                                                                                                         >  2:1127907/1‑148 (MQ=255)
                                         cagACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgcg                                                                                                       <  1:2557429/146‑1 (MQ=255)
                                          cGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                                      >  2:1028542/1‑148 (MQ=255)
                                          cGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                                      >  1:519282/1‑148 (MQ=255)
                                          cGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCgg                                                                                                      <  2:2013029/148‑1 (MQ=255)
                                            aCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGTAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgt                                                                                                    <  2:429736/148‑1 (MQ=255)
                                             cTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtg                                                                                                   >  1:3457054/1‑148 (MQ=255)
                                             cTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtg                                                                                                   <  1:536015/148‑1 (MQ=255)
                                              ttCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtgt                                                                                                  >  2:533405/1‑148 (MQ=255)
                                                   aTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                             >  2:3565584/1‑148 (MQ=255)
                                                     gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGcc                                                                                                                     >  1:2827696/1‑122 (MQ=255)
                                                     gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAg                                                                                           <  2:688747/148‑1 (MQ=255)
                                                     gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAg                                                                                           >  2:2014842/1‑148 (MQ=255)
                                                     gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAg                                                                                           >  1:3203918/1‑148 (MQ=255)
                                                      cATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                                          >  2:2648858/1‑148 (MQ=255)
                                                         gATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATca                                                                                       >  1:1657526/1‑147 (MQ=255)
                                                          aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                      >  1:489358/1‑148 (MQ=255)
                                                          aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                      <  1:1398009/148‑1 (MQ=255)
                                                            cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                    >  1:672053/1‑148 (MQ=255)
                                                            cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                    <  2:2592871/148‑1 (MQ=255)
                                                                 cATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTAc                                                                               >  1:1735235/1‑148 (MQ=255)
                                                                      tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCa                                                                          <  2:2318343/148‑1 (MQ=255)
                                                                       cacaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCAt                                                                         >  1:240384/1‑148 (MQ=255)
                                                                        acaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                        <  1:2578515/148‑1 (MQ=255)
                                                                        acaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                        <  1:1828432/148‑1 (MQ=255)
                                                                            ccTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcgc                                                                    >  1:1323296/1‑148 (MQ=255)
                                                                               tCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGa                                                                 <  1:856216/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              <  2:2803648/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              >  2:1047756/1‑148 (MQ=255)
                                                                                    cttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATc                                                            >  2:2627517/1‑148 (MQ=255)
                                                                                        tcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACg                                                        <  2:1642139/148‑1 (MQ=255)
                                                                                        tcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACg                                                        >  1:3503276/1‑148 (MQ=255)
                                                                                         cGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGc                                                       <  1:3071114/148‑1 (MQ=255)
                                                                                           tACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                                     >  2:2548823/1‑148 (MQ=255)
                                                                                           tACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                                     <  2:1381101/148‑1 (MQ=255)
                                                                                                 tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTcac                                               >  2:3324949/1‑148 (MQ=255)
                                                                                                 tCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTcac                                               >  1:1581267/1‑148 (MQ=255)
                                                                                                     tgtgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                                           <  1:93043/148‑1 (MQ=255)
                                                                                                       tgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGa                                         >  1:389514/1‑148 (MQ=255)
                                                                                                          cGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCa                                      >  2:52752/1‑148 (MQ=255)
                                                                                                            aaTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACt                                    >  2:3447007/1‑148 (MQ=255)
                                                                                                             aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                   >  2:2821635/1‑148 (MQ=255)
                                                                                                             aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                   <  1:2047854/148‑1 (MQ=255)
                                                                                                               cACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa                                 >  2:1244197/1‑148 (MQ=255)
                                                                                                                aCTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAg                                >  1:1182479/1‑148 (MQ=255)
                                                                                                                aCTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAg                                >  1:3273713/1‑148 (MQ=255)
                                                                                                                   aCGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAc                             >  2:839790/1‑148 (MQ=255)
                                                                                                                          tGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGAt                      <  1:3101585/148‑1 (MQ=255)
                                                                                                                            ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATaa                    <  2:1521080/148‑1 (MQ=255)
                                                                                                                            ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATaa                    <  2:2249502/148‑1 (MQ=255)
                                                                                                                              ccAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                  <  1:2482485/148‑1 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               <  1:1699326/148‑1 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               <  1:3565584/148‑1 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTAAAGACAGACCGCGGAATCTAGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               <  1:3497566/148‑1 (MQ=255)
                                                                                                                                    gTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATgg            <  1:111363/148‑1 (MQ=255)
                                                                                                                                     tGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgc           >  1:1023362/1‑148 (MQ=255)
                                                                                                                                       aTAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgcgc         >  2:2144945/1‑148 (MQ=255)
                                                                                                                                        tAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGcgcg        >  1:2246149/1‑148 (MQ=255)
                                                                                                                                         aGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCgg       <  2:2620379/148‑1 (MQ=255)
                                                                                                                                          ggAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGt      >  1:2926215/1‑148 (MQ=255)
                                                                                                                                              agataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgt  <  2:2078695/148‑1 (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: