Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,458,948 G→T 100% D503Y (GAT→TAT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,9480GT98.4% 602.1 / 0.8 188D503Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (91/94);  total (92/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 189 total aligned reads displayed.

GGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAA  >  CP009273/4458803‑4459092
                                                                                                                                                 |                                                                                                                                                
ggtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTAt                                                                                                                                                <  1:1036959/148‑1 (MQ=255)
  tAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATAc                                                                                                                                              >  1:260189/1‑148 (MQ=255)
   aGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACg                                                                                                                                             <  1:1264248/148‑1 (MQ=255)
   aGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACg                                                                                                                                             <  1:2544170/148‑1 (MQ=255)
       gCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGa                                                                                                                                         >  2:429119/1‑148 (MQ=255)
        cAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGaa                                                                                                                                        <  1:2406096/148‑1 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    <  2:2434364/148‑1 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    >  2:747509/1‑148 (MQ=255)
             aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTc                                                                                                                                   <  2:540327/148‑1 (MQ=255)
             aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTc                                                                                                                                   <  1:2122632/148‑1 (MQ=255)
              aTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCa                                                                                                                                  >  2:1546173/1‑148 (MQ=255)
              aTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCa                                                                                                                                  >  2:919059/1‑148 (MQ=255)
                aCCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgg                                                                                                                                <  1:3341324/148‑1 (MQ=255)
                 ccGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAggg                                                                                                                               >  1:2193513/1‑148 (MQ=255)
                 ccGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAggg                                                                                                                               <  2:1724505/148‑1 (MQ=255)
                 ccGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGAGGGTTACTTATACGCTGAATCGTCAggg                                                                                                                               >  1:403537/1‑148 (MQ=255)
                  cGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgggg                                                                                                                              <  2:2868266/148‑1 (MQ=255)
                   gAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGc                                                                                                                             <  1:3103017/148‑1 (MQ=255)
                         cATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGcgcg                                                                                                                       >  2:2087304/1‑148 (MQ=255)
                         cATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGcgcg                                                                                                                       <  2:2354354/148‑1 (MQ=255)
                         cATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGcgcg                                                                                                                       >  2:2388144/1‑148 (MQ=255)
                              gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  >  2:304223/1‑148 (MQ=255)
                              gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  <  1:2265526/148‑1 (MQ=255)
                              gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  >  1:477686/1‑148 (MQ=255)
                              gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  <  1:3784/148‑1 (MQ=255)
                              gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGATGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  >  2:2200702/1‑148 (MQ=255)
                               cTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTg                                                                                                                 >  2:56545/1‑148 (MQ=255)
                               cTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTg                                                                                                                 >  2:507829/1‑148 (MQ=255)
                               cTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTg                                                                                                                 >  1:3472227/1‑148 (MQ=255)
                                tGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGc                                                                                                                >  2:2278462/1‑148 (MQ=255)
                                tGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGc                                                                                                                >  2:2278874/1‑148 (MQ=255)
                                  gTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCgg                                                                                                              >  2:1685798/1‑148 (MQ=255)
                                  gTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCgg                                                                                                              >  2:2619839/1‑148 (MQ=255)
                                   tttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGt                                                                                                             >  1:1998256/1‑148 (MQ=255)
                                         aaGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcgac                                                                                                       <  1:49181/148‑1 (MQ=255)
                                           ggTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcgacga                                                                                                     >  2:1096138/1‑148 (MQ=255)
                                            gTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaa                                                                                                    >  2:416242/1‑148 (MQ=255)
                                                aTTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAc                                                                                                >  1:1831498/1‑148 (MQ=255)
                                                aTTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAc                                                                                                <  1:1323124/148‑1 (MQ=255)
                                                 ttGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAcc                                                                                               <  2:719845/148‑1 (MQ=255)
                                                 ttGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTAGACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAcc                                                                                               <  1:932737/148‑1 (MQ=255)
                                                  tGCAGGAAATCCACAATGTGTGATCGCAGGTGCGTCACTATGGCTAGATTGTCCCGCTCGATGACATCATGATGCGTAAGATTAAGCGGGTTACTTATAAGCTGAATCGGCAGGGGCTTCGCGTGTTTGCGGTGGCGACGAAATACCt                                                                                              >  2:448867/1‑148 (MQ=255)
                                                   gCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTg                                                                                             >  1:2915743/1‑148 (MQ=255)
                                                    cAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGc                                                                                            <  2:3458081/148‑1 (MQ=255)
                                                     aGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGcc                                                                                           >  2:3243/1‑148 (MQ=255)
                                                      ggAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCa                                                                                          >  1:712652/1‑148 (MQ=255)
                                                      ggAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCa                                                                                          >  1:2495798/1‑148 (MQ=255)
                                                         aaTCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAgcg                                                                                       <  2:2796344/148‑1 (MQ=255)
                                                            ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCTTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                    >  1:1694936/1‑148 (MQ=255)
                                                            ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                    <  1:2885009/148‑1 (MQ=255)
                                                            ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                    >  1:2963763/1‑148 (MQ=255)
                                                              tCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                  >  2:119424/1‑148 (MQ=255)
                                                              tCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAACCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                  <  1:2671855/148‑1 (MQ=255)
                                                               cAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGaa                                                                                 <  2:2633638/148‑1 (MQ=255)
                                                                 aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGACATACCTGCCAGCGCGTGAAgg                                                                               >  2:2365396/1‑148 (MQ=255)
                                                                     gtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                           >  1:39633/1‑148 (MQ=255)
                                                                              ggTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgg                                                                  <  2:1016762/148‑1 (MQ=255)
                                                                               gTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGAGCGTGAAGGAGATTACCAgcggg                                                                 >  2:1827467/1‑148 (MQ=255)
                                                                                  cGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgggcgg                                                              <  2:375615/148‑1 (MQ=255)
                                                                                    tCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGAt                                                            >  2:1508331/1‑148 (MQ=255)
                                                                                          tGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATcc                                                      <  1:2514034/148‑1 (MQ=255)
                                                                                            gCGAGTTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGa                                                    <  2:619135/148‑1 (MQ=255)
                                                                                            gCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGa                                                    >  1:72584/1‑148 (MQ=255)
                                                                                             cGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAc                                                   <  1:2552574/148‑1 (MQ=255)
                                                                                               agaTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCt                                                 >  2:547710/1‑148 (MQ=255)
                                                                                                 aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                               >  1:3232379/1‑148 (MQ=255)
                                                                                                 aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                               <  2:117622/148‑1 (MQ=255)
                                                                                                   tgtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                             <  1:1207159/148‑1 (MQ=255)
                                                                                                    gtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATcc                                            <  1:1045292/148‑1 (MQ=255)
                                                                                                     tgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCt                                           >  1:2887216/1‑148 (MQ=255)
                                                                                                        cGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGa                                        <  2:1320774/148‑1 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                       >  2:2065446/1‑148 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                       <  2:745201/148‑1 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                       <  2:3162460/148‑1 (MQ=255)
                                                                                                           tCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAgg                                     >  1:801454/1‑148 (MQ=255)
                                                                                                             gATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGat                                   <  1:441297/148‑1 (MQ=255)
                                                                                                              aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  >  2:1980879/1‑148 (MQ=255)
                                                                                                              aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  <  2:625677/148‑1 (MQ=255)
                                                                                                              aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  <  2:2917540/148‑1 (MQ=255)
                                                                                                              aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  <  2:3016386/148‑1 (MQ=255)
                                                                                                              aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  >  2:1881722/1‑148 (MQ=255)
                                                                                                                gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatata                                >  1:1923052/1‑148 (MQ=255)
                                                                                                                gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatata                                <  2:2660190/148‑1 (MQ=255)
                                                                                                                  catcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATAtt                              >  1:2783777/1‑148 (MQ=255)
                                                                                                                  catcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATAtt                              >  2:1009749/1‑148 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGTATCGTCAGGGGCTGCGCGTGGGTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGTCCTGATCCTCGAAGGATATATTGc                            <  2:2224386/148‑1 (MQ=255)
                                                                                                                       tgctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCttt                         >  2:1728692/1‑148 (MQ=255)
                                                                                                                       tgctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCttt                         >  1:2868380/1‑148 (MQ=255)
                                                                                                                                aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                >  2:2538233/1‑148 (MQ=255)
                                                                                                                                  aTTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGc              <  1:1996040/148‑1 (MQ=255)
                                                                                                                                  aTTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGc              >  1:1464721/1‑148 (MQ=255)
                                                                                                                                   ttAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGcc             >  1:574552/1‑148 (MQ=255)
                                                                                                                                    tAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa            >  1:671325/1‑148 (MQ=255)
                                                                                                                                    tAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa            >  2:1940411/1‑148 (MQ=255)
                                                                                                                                          ggTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGTAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGTTCCTAGAAGGATATATTGCTTTTCTTGATCCGCCAAAagag      <  2:2653158/148‑1 (MQ=255)
                                                                                                                                           gTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAgaga     <  2:936031/148‑1 (MQ=255)
                                                                                                                                            ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGAc    <  2:3411022/148‑1 (MQ=255)
                                                                                                                                             tACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACa   <  1:747509/148‑1 (MQ=255)
                                                                                                                                              aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa  <  2:176950/148‑1 (MQ=255)
                                                                                                                                              aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa  <  1:2556062/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAA  >  CP009273/4458803‑4459092

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: ‑