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,0540GT95.9% 665.7 / 17.9 220S222Y (TCT→TAT) sthApyridine nucleotide transhydrogenase, soluble
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (111/100);  total (115/105)
Fisher's exact test for biased strand distribution p-value = 7.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

ACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGT  >  CP009273/4149909‑4150198
                                                                                                                                                 |                                                                                                                                                
aCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATaa                                                                                                                                                >  2:1038084/1‑148 (MQ=255)
aCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATaa                                                                                                                                                >  2:505139/1‑148 (MQ=255)
 ccGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAAt                                                                                                                                               >  1:1548214/1‑148 (MQ=255)
    ttGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTg                                                                                                                                            >  1:117564/1‑148 (MQ=255)
         aTAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAtc                                                                                                                                       <  2:2617783/148‑1 (MQ=255)
         aTAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCGTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAtc                                                                                                                                       <  2:2635491/148‑1 (MQ=255)
          tAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACAtct                                                                                                                                      <  1:2963095/148‑1 (MQ=255)
          tAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCACACAtct                                                                                                                                      >  1:2500712/1‑148 (MQ=255)
          tAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACAtct                                                                                                                                      <  1:1151036/148‑1 (MQ=255)
            gagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATctct                                                                                                                                    <  1:1284636/148‑1 (MQ=255)
            gagCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATctct                                                                                                                                    <  1:3025074/148‑1 (MQ=255)
                 ggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATc                                                                                                                               >  2:1650773/1‑148 (MQ=255)
                  gcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCg                                                                                                                              >  1:1399217/1‑148 (MQ=255)
                       cAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGaaa                                                                                                                         <  1:2974067/148‑1 (MQ=255)
                         gCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATg                                                                                                                       >  1:1686055/1‑148 (MQ=255)
                          cTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGc                                                                                                                      >  2:1039068/1‑148 (MQ=255)
                                   ttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCg                                                                                                             <  1:866153/148‑1 (MQ=255)
                                    tttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGa                                                                                                            >  1:812478/1‑148 (MQ=255)
                                     ttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGAt                                                                                                           <  1:2911763/148‑1 (MQ=255)
                                        cccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgc                                                                                                        <  2:604339/148‑1 (MQ=255)
                                         ccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATcgcg                                                                                                       >  2:2818830/1‑148 (MQ=255)
                                           gACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCggg                                                                                                     <  1:760051/148‑1 (MQ=255)
                                           gACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCggg                                                                                                     <  2:3090560/148‑1 (MQ=255)
                                            aCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgt                                                                                                    >  2:1203742/1‑148 (MQ=255)
                                            aCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgt                                                                                                    >  2:2584017/1‑148 (MQ=255)
                                              ttCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGgtgt                                                                                                  <  2:1456663/148‑1 (MQ=255)
                                               tCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGtt                                                                                                 >  2:1544123/1‑148 (MQ=255)
                                                cAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTg                                                                                                <  1:3218573/148‑1 (MQ=255)
                                                 aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGa                                                                                               >  1:2132422/1‑148 (MQ=255)
                                                  gATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGAt                                                                                              <  1:1038499/148‑1 (MQ=255)
                                                   aTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATc                                                                                             >  2:1358396/1‑148 (MQ=255)
                                                    tGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCa                                                                                            <  1:884240/148‑1 (MQ=255)
                                                     gCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAg                                                                                           >  2:2083246/1‑148 (MQ=255)
                                                      cATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGa                                                                                          <  2:1991246/148‑1 (MQ=255)
                                                       aTGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGAt                                                                                         >  2:2665384/1‑148 (MQ=255)
                                                        tGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATc                                                                                        >  1:46927/1‑148 (MQ=255)
                                                          aTCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                      >  2:1326753/1‑148 (MQ=255)
                                                           tCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCAc                                                                                     >  2:3055441/1‑148 (MQ=255)
                                                            cacaCCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                    <  1:1976625/148‑1 (MQ=255)
                                                             acacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtt                                                                                   <  2:1537802/148‑1 (MQ=255)
                                                             acacCATCGTCACAGCCTTCGATCTTCTCGGACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTGGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtt                                                                                   <  2:3216695/148‑1 (MQ=255)
                                                               acCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACtttt                                                                                 >  1:1524996/1‑148 (MQ=255)
                                                                ccATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTa                                                                                >  2:756822/1‑148 (MQ=255)
                                                                  aTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACa                                                                              >  1:1856902/1‑148 (MQ=255)
                                                                   tCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTGTACAt                                                                             <  2:2305797/148‑1 (MQ=255)
                                                                    cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                            <  2:11596/148‑1 (MQ=255)
                                                                    cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                            >  2:1282850/1‑148 (MQ=255)
                                                                    cGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                            >  1:345492/1‑148 (MQ=255)
                                                                     gTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATcc                                                                           <  2:564859/148‑1 (MQ=255)
                                                                       cacaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCAt                                                                         >  1:1531162/1‑148 (MQ=255)
                                                                        acaGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATa                                                                        >  2:3427404/1‑148 (MQ=255)
                                                                         caGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAc                                                                       >  2:314886/1‑148 (MQ=255)
                                                                          aGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATAcc                                                                      >  1:3039874/1‑148 (MQ=255)
                                                                           gCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcg                                                                     >  1:54816/1‑148 (MQ=255)
                                                                                cGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATACTCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGaa                                                                >  2:536845/1‑148 (MQ=255)
                                                                                 gATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAg                                                               <  2:870523/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              <  2:767753/148‑1 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              >  1:1250351/1‑148 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              >  1:3021326/1‑148 (MQ=255)
                                                                                  aTCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGa                                                              <  1:3191661/148‑1 (MQ=255)
                                                                                    cttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATc                                                            >  1:601687/1‑148 (MQ=255)
                                                                                     ttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCg                                                           <  1:1480287/148‑1 (MQ=255)
                                                                                     ttctCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCg                                                           <  1:2943722/148‑1 (MQ=255)
                                                                                           tACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                                     <  1:1429965/148‑1 (MQ=255)
                                                                                            aCTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCata                                                    <  2:2993206/148‑1 (MQ=255)
                                                                                            aCTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCata                                                    >  2:2738231/1‑148 (MQ=255)
                                                                                             ctctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCatat                                                   >  2:3192861/1‑148 (MQ=255)
                                                                                               ctTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTc                                                 <  2:2389082/148‑1 (MQ=255)
                                                                                                      gtgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCg                                          >  2:3176738/1‑148 (MQ=255)
                                                                                                       tgACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGa                                         <  2:3089662/148‑1 (MQ=255)
                                                                                                        gACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGAt                                        <  1:1388631/148‑1 (MQ=255)
                                                                                                        gACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGAt                                        >  1:332533/1‑148 (MQ=255)
                                                                                                        gACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGAt                                        <  2:1759582/148‑1 (MQ=255)
                                                                                                         aCGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                                       >  2:2016781/1‑148 (MQ=255)
                                                                                                          cGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCa                                      >  1:657619/1‑148 (MQ=255)
                                                                                                            aaTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCAca                                    >  1:2351165/1‑147 (MQ=255)
                                                                                                             aTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                   <  1:2085900/148‑1 (MQ=255)
                                                                                                              tCACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTcc                                  >  2:1704512/1‑148 (MQ=255)
                                                                                                               cACTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCa                                 <  1:1319515/148‑1 (MQ=255)
                                                                                                                 cTACGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGc                               >  1:2629927/1‑148 (MQ=255)
                                                                                                                   aCGCCACTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAc                             <  2:2178256/148‑1 (MQ=255)
                                                                                                                        aCTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAg                        <  1:2567982/148‑1 (MQ=255)
                                                                                                                         cTGTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGa                       >  2:1798565/1‑148 (MQ=255)
                                                                                                                           gTTCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATa                     >  2:36015/1‑148 (MQ=255)
                                                                                                                            ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATaa                    <  1:3328086/148‑1 (MQ=255)
                                                                                                                            ttCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATaa                    <  1:1023155/148‑1 (MQ=255)
                                                                                                                             tCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAg                   >  1:1766804/1‑148 (MQ=255)
                                                                                                                             tCCAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAg                   >  1:3462235/1‑148 (MQ=255)
                                                                                                                              ccAGAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                  <  2:345492/148‑1 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               >  1:2410679/1‑148 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               <  1:1767044/148‑1 (MQ=255)
                                                                                                                                 gAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACa               <  1:571406/148‑1 (MQ=255)
                                                                                                                                   aGTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATg             <  1:1563601/148‑1 (MQ=255)
                                                                                                                                    gTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATgg            <  2:1256172/148‑1 (MQ=255)
                                                                                                                                    gTGATAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATgg            <  1:510322/148‑1 (MQ=255)
                                                                                                                                       aTAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgcgc         >  2:3244148/1‑148 (MQ=255)
                                                                                                                                       aTAGGAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGgcgc         >  2:492043/1‑148 (MQ=255)
                                                                                                                                          ggAGAGATAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGt      <  1:413998/148‑1 (MQ=255)
                                                                                                                                           gagagataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGtt     >  1:1619485/1‑148 (MQ=255)
                                                                                                                                              agataATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgt  <  1:2713470/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: