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,9480GT99.5% 696.4 / ‑4.3 213D503Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (99/113);  total (99/114)
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.89e-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 213 total aligned reads displayed.

GGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAA  >  CP009273/4458803‑4459092
                                                                                                                                                 |                                                                                                                                                
ggtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGAAAAGGTTCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGAGGGTTACTTAt                                                                                                                                                >  1:1010559/1‑148 (MQ=255)
   aGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACg                                                                                                                                             <  2:3783896/148‑1 (MQ=255)
       gCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGa                                                                                                                                         <  1:797006/148‑1 (MQ=255)
        cAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGaa                                                                                                                                        <  1:2818967/148‑1 (MQ=255)
         aGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAAt                                                                                                                                       <  1:2225183/148‑1 (MQ=255)
           aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                     >  1:3769195/1‑148 (MQ=255)
           aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTAATTATACGCTGAATCg                                                                                                                                     <  1:3113969/148‑1 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    >  1:3345941/1‑148 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    >  2:2771712/1‑148 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    >  2:2774902/1‑148 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    >  2:3017970/1‑148 (MQ=255)
            aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGt                                                                                                                                    <  1:392508/148‑1 (MQ=255)
             aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTc                                                                                                                                   >  2:2247949/1‑148 (MQ=255)
             aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTc                                                                                                                                   >  2:3629746/1‑148 (MQ=255)
               tACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAg                                                                                                                                 <  1:2523721/148‑1 (MQ=255)
                aCCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgg                                                                                                                                >  1:2918353/1‑148 (MQ=255)
                  cGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgggg                                                                                                                              <  2:312316/148‑1 (MQ=255)
                   gAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGc                                                                                                                             <  2:2109445/148‑1 (MQ=255)
                        ccATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTgcgc                                                                                                                        >  1:2085332/1‑148 (MQ=255)
                           tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTg                                                                                                                     <  1:1375610/148‑1 (MQ=255)
                            cAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTgg                                                                                                                    >  2:470945/1‑148 (MQ=255)
                              gCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGtt                                                                                                                  >  2:2992874/1‑148 (MQ=255)
                               cTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTg                                                                                                                 >  2:2866845/1‑148 (MQ=255)
                                    ttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                            >  2:34434/1‑148 (MQ=255)
                                         aaGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcgac                                                                                                       >  1:1074846/1‑148 (MQ=255)
                                         aaGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcgac                                                                                                       >  2:1313121/1‑148 (MQ=255)
                                              gCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAAt                                                                                                  <  1:1768373/148‑1 (MQ=255)
                                               cATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATa                                                                                                 <  1:1454666/148‑1 (MQ=255)
                                               cATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATa                                                                                                 <  2:3553596/148‑1 (MQ=255)
                                               cATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATa                                                                                                 >  1:3407410/1‑148 (MQ=255)
                                                aTTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAc                                                                                                <  2:3462016/148‑1 (MQ=255)
                                                 ttGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAcc                                                                                               <  2:2085332/148‑1 (MQ=255)
                                                    cAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGc                                                                                            <  2:1600184/148‑1 (MQ=255)
                                                         aaTCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAgcg                                                                                       <  1:1682267/148‑1 (MQ=255)
                                                           tCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGcgcg                                                                                     >  2:1899098/1‑148 (MQ=255)
                                                            ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCTTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                    >  2:996685/1‑148 (MQ=255)
                                                              tCAATGTGTGTTCGCAGGTGCGTCACACTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                  >  2:1047023/1‑148 (MQ=255)
                                                              tCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                  >  1:3659289/1‑148 (MQ=255)
                                                               cAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGaa                                                                                 <  2:1453820/148‑1 (MQ=255)
                                                                 aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAgg                                                                               >  2:687892/1‑148 (MQ=255)
                                                                 aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAgg                                                                               <  1:898962/148‑1 (MQ=255)
                                                                 aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAgg                                                                               <  1:2058470/148‑1 (MQ=255)
                                                                 aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAgg                                                                               <  2:882663/148‑1 (MQ=255)
                                                                   gtgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGgag                                                                             >  1:98406/1‑148 (MQ=255)
                                                                     gtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                           >  1:48487/1‑148 (MQ=255)
                                                                      tgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAtt                                                                          >  1:1950514/1‑148 (MQ=255)
                                                                       gtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTa                                                                         <  2:2309045/148‑1 (MQ=255)
                                                                          cGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCa                                                                      <  2:2614017/148‑1 (MQ=255)
                                                                              ggTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgg                                                                  >  2:1762242/1‑148 (MQ=255)
                                                                              ggTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgg                                                                  >  1:1955923/1‑148 (MQ=255)
                                                                              ggTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgg                                                                  <  1:2548547/148‑1 (MQ=255)
                                                                               gTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGACGGAGATTACCAgcggg                                                                 >  2:2688767/1‑148 (MQ=255)
                                                                               gTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcggg                                                                 >  1:1764822/1‑148 (MQ=255)
                                                                                  cGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgggcgg                                                              <  2:3550315/148‑1 (MQ=255)
                                                                                     cacaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATg                                                           >  1:3578255/1‑148 (MQ=255)
                                                                                      acaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACTCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGa                                                          <  2:1437196/148‑1 (MQ=255)
                                                                                      acaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGa                                                          <  2:243219/148‑1 (MQ=255)
                                                                                       caATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGaa                                                         <  1:105996/148‑1 (MQ=255)
                                                                                       caATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGaa                                                         <  1:3456224/148‑1 (MQ=255)
                                                                                         aTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATc                                                       >  2:1658762/1‑148 (MQ=255)
                                                                                                 aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                               <  1:2686761/148‑1 (MQ=255)
                                                                                                   tgtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                             >  1:1705857/1‑148 (MQ=255)
                                                                                                   tgtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                             >  1:1090369/1‑148 (MQ=255)
                                                                                                     tgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCt                                           <  1:34434/148‑1 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                       <  2:1225560/148‑1 (MQ=255)
                                                                                                           tCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAgg                                     <  1:2135392/148‑1 (MQ=255)
                                                                                                            cGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGa                                    <  2:2034803/148‑1 (MQ=255)
                                                                                                             gATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGat                                   >  1:2713805/1‑148 (MQ=255)
                                                                                                              aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGata                                  >  1:552020/1‑148 (MQ=255)
                                                                                                               tGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatat                                 >  2:2885396/1‑148 (MQ=255)
                                                                                                                gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatata                                <  1:1244172/148‑1 (MQ=255)
                                                                                                                gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatata                                <  1:91843/148‑1 (MQ=255)
                                                                                                                 aCATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatatat                               >  2:1346909/1‑148 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                            <  1:1539028/148‑1 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                            <  1:3017963/148‑1 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                            <  1:2412978/148‑1 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                            <  1:2413324/148‑1 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                            <  1:2774902/148‑1 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                            <  1:2771715/148‑1 (MQ=255)
                                                                                                                     catGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCt                           >  2:270612/1‑148 (MQ=255)
                                                                                                                      atGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtt                          <  1:1184531/148‑1 (MQ=255)
                                                                                                                       tgctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCttt                         >  1:531763/1‑148 (MQ=255)
                                                                                                                        gctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCCGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtttt                        >  2:2135933/1‑148 (MQ=255)
                                                                                                                        gctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtttt                        >  1:1550276/1‑148 (MQ=255)
                                                                                                                        gctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtttt                        >  1:1549138/1‑148 (MQ=255)
                                                                                                                         ctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTc                       <  2:1173168/148‑1 (MQ=255)
                                                                                                                         ctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTc                       <  2:442267/148‑1 (MQ=255)
                                                                                                                         ctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTc                       <  2:982444/148‑1 (MQ=255)
                                                                                                                            cGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTg                    <  2:808825/148‑1 (MQ=255)
                                                                                                                            cGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTg                    >  1:777643/1‑148 (MQ=255)
                                                                                                                            cGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTg                    <  2:3659332/148‑1 (MQ=255)
                                                                                                                              tAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                  >  1:1610274/1‑148 (MQ=255)
                                                                                                                                aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                >  2:292569/1‑148 (MQ=255)
                                                                                                                                  aTTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGc              >  2:577546/1‑148 (MQ=255)
                                                                                                                                   ttAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGcc             >  2:1663704/1‑148 (MQ=255)
                                                                                                                                    tAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa            >  1:3286223/1‑148 (MQ=255)
                                                                                                                                    tAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa            >  1:3686089/1‑148 (MQ=255)
                                                                                                                                        cGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAag        <  2:1572102/148‑1 (MQ=255)
                                                                                                                                            ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGAc    >  2:2203396/1‑148 (MQ=255)
                                                                                                                                              aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa  <  1:601085/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: