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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,6090GT99.5% 690.3 / ‑4.1 208D390Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (123/84);  total (123/85)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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 210 total aligned reads displayed.

GGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC  >  CP009273/4458464‑4458753
                                                                                                                                                 |                                                                                                                                                
 gcgcgcGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCtata                                                                                                                                               <  2:2672516/148‑1 (MQ=255)
    cgcgGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtct                                                                                                                                            >  1:1047679/1‑148 (MQ=255)
     gcggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtctc                                                                                                                                           <  2:1282493/148‑1 (MQ=255)
     gcggcgCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATAtctc                                                                                                                                           <  1:2542987/148‑1 (MQ=255)
        gcgcAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCgg                                                                                                                                        >  2:2944975/1‑148 (MQ=255)
             gTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTaaaa                                                                                                                                   <  2:3571776/148‑1 (MQ=255)
              tAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAc                                                                                                                                  >  1:2154508/1‑148 (MQ=255)
               aaaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAAcc                                                                                                                                 <  2:324947/148‑1 (MQ=255)
                aaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCa                                                                                                                                >  2:2042212/1‑148 (MQ=255)
                aaGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCa                                                                                                                                <  1:2863560/148‑1 (MQ=255)
                aaGCTGTCGAAACAGAAAGTGATCGTCAAAAATCTGGATGCTATTCAGAACTTTGGCGCATTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGTATAAAATTGTGCTGGAGAATCATACCTATATCTCCGTTAAAACCa                                                                                                                                <  1:1946275/148‑1 (MQ=255)
                 aGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAg                                                                                                                               >  1:2838104/1‑148 (MQ=255)
                  gCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGc                                                                                                                              >  1:72936/1‑148 (MQ=255)
                    tGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGa                                                                                                                            <  1:3328987/148‑1 (MQ=255)
                     gTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGaa                                                                                                                           <  2:2695474/148‑1 (MQ=255)
                       cGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAAcg                                                                                                                         >  2:1961159/1‑148 (MQ=255)
                           aCAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtg                                                                                                                     <  1:2127936/148‑1 (MQ=255)
                               aaaGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGtgctgc                                                                                                                 >  2:837765/1‑148 (MQ=255)
                                 aGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACTCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACAAGCGAACGCGTGCAGCAt                                                                                                               >  2:2963906/1‑148 (MQ=255)
                                 aGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCAt                                                                                                               >  2:3692323/1‑148 (MQ=255)
                                  gTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATa                                                                                                              >  2:2117971/1‑148 (MQ=255)
                                   tGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAg                                                                                                             >  2:1363081/1‑148 (MQ=255)
                                   tGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAg                                                                                                             >  2:1772396/1‑148 (MQ=255)
                                    gATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                            >  1:2403449/1‑148 (MQ=255)
                                    gATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgc                                                                                                            >  1:1417529/1‑148 (MQ=255)
                                     aTCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgcg                                                                                                           >  1:720023/1‑148 (MQ=255)
                                      tCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAgcgc                                                                                                          >  1:320698/1‑148 (MQ=255)
                                       cGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGcgcg                                                                                                         >  1:3135841/1‑148 (MQ=255)
                                        gTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGt                                                                                                        <  1:2283221/148‑1 (MQ=255)
                                         tCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTg                                                                                                       <  1:2025788/148‑1 (MQ=255)
                                           aaaCATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGt                                                                                                     >  2:839652/1‑148 (MQ=255)
                                             aCATCAGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACACTGACGCAGTATAAAATAGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGAATAGCGCGTGGTTg                                                                                                   >  2:3038855/1‑148 (MQ=255)
                                              cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGa                                                                                                  >  1:3005212/1‑148 (MQ=255)
                                              cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGa                                                                                                  >  2:2388943/1‑148 (MQ=255)
                                              cATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGa                                                                                                  >  2:1831367/1‑148 (MQ=255)
                                               aTCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGAACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGaa                                                                                                 >  2:540121/1‑148 (MQ=255)
                                                tCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAAc                                                                                                >  1:3602052/1‑148 (MQ=255)
                                                 cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACa                                                                                               >  1:1380001/1‑148 (MQ=255)
                                                 cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACa                                                                                               >  2:1027967/1‑148 (MQ=255)
                                                 cTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACa                                                                                               <  1:814530/148‑1 (MQ=255)
                                                  tGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAg                                                                                              >  1:1171381/1‑148 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                             >  1:1073381/1‑148 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                             <  1:41389/148‑1 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                             >  1:1721391/1‑148 (MQ=255)
                                                   ggATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGACACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGt                                                                                             >  2:148729/1‑148 (MQ=255)
                                                     aTGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGtga                                                                                           <  2:2926475/148‑3 (MQ=255)
                                                       gCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCAtt                                                                                         >  2:1408487/1‑148 (MQ=255)
                                                        cTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTa                                                                                        <  1:612491/148‑1 (MQ=255)
                                                         tATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTAt                                                                                       <  2:1915121/148‑1 (MQ=255)
                                                          aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATc                                                                                      >  1:1697644/1‑148 (MQ=255)
                                                          aTTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATc                                                                                      >  1:3404768/1‑148 (MQ=255)
                                                             cAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGa                                                                                   <  1:1117005/148‑1 (MQ=255)
                                                               gAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGAcc                                                                                 >  1:1192271/1‑148 (MQ=255)
                                                                  cTTTGTCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGa                                                                              <  2:3577945/148‑1 (MQ=255)
                                                                   tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGAc                                                                             >  2:3495375/1‑148 (MQ=255)
                                                                   tttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGAc                                                                             >  2:287610/1‑148 (MQ=255)
                                                                    ttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACt                                                                            >  2:3431987/1‑148 (MQ=255)
                                                                    ttGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACt                                                                            >  2:423266/1‑148 (MQ=255)
                                                                          cAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaaa                                                                      >  2:3658401/1‑148 (MQ=255)
                                                                          cAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTaaaaa                                                                      >  2:2416112/1‑148 (MQ=255)
                                                                           aaTGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAAc                                                                     <  1:1772396/148‑1 (MQ=255)
                                                                              ggATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTg                                                                  <  2:1870768/148‑1 (MQ=255)
                                                                              ggATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTg                                                                  >  2:2147091/1‑148 (MQ=255)
                                                                                  atTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCg                                                              >  2:1218926/1‑148 (MQ=255)
                                                                                   ttCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGa                                                             >  1:462817/1‑148 (MQ=255)
                                                                                      tgtgCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATAc                                                          <  1:2337910/148‑1 (MQ=255)
                                                                                          cACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCg                                                      <  2:1094329/148‑1 (MQ=255)
                                                                                          cACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCg                                                      <  2:1096749/148‑1 (MQ=255)
                                                                                              gATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGc                                                  <  2:2207356/148‑1 (MQ=255)
                                                                                               aTAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCt                                                 >  2:1104205/1‑148 (MQ=255)
                                                                                                  aaaCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGa                                              <  2:3514129/148‑1 (MQ=255)
                                                                                                    aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                            >  1:780997/1‑148 (MQ=255)
                                                                                                    aCCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAg                                            >  1:46500/1‑148 (MQ=255)
                                                                                                         cACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACg                                       >  2:3213302/1‑148 (MQ=255)
                                                                                                           cccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGa                                     >  1:1243900/1‑148 (MQ=255)
                                                                                                            ccTGACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGAt                                    <  2:748953/148‑1 (MQ=255)
                                                                                                               gACGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGaa                                 <  1:1374957/148‑1 (MQ=255)
                                                                                                                aCGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAag                                <  2:719433/148‑1 (MQ=255)
                                                                                                                aCGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAag                                <  2:720023/148‑1 (MQ=255)
                                                                                                                 cGCAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAaga                               >  2:2271758/1‑148 (MQ=255)
                                                                                                                   cAGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGt                             <  2:1073381/148‑1 (MQ=255)
                                                                                                                    aGGATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTc                            <  2:1617271/148‑1 (MQ=255)
                                                                                                                      gATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAg                          >  2:13219/1‑148 (MQ=255)
                                                                                                                      gATAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAg                          >  2:3562264/1‑148 (MQ=255)
                                                                                                                       aTAAAATTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgc                         >  1:505173/1‑148 (MQ=255)
                                                                                                                           aaTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAgcgcgc                     >  1:853152/1‑148 (MQ=255)
                                                                                                                            aTTGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgct                    >  1:3553596/1‑148 (MQ=255)
                                                                                                                             ttGTGCTGGAGAATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGcgctc                   >  1:3107178/1‑148 (MQ=255)
                                                                                                                                      agaATCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAg          >  2:797006/1‑148 (MQ=255)
                                                                                                                                        aaTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTc        >  2:2648141/1‑148 (MQ=255)
                                                                                                                                         aTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCg       <  1:2042151/148‑1 (MQ=255)
                                                                                                                                         aTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCg       <  1:2042259/148‑1 (MQ=255)
                                                                                                                                         aTCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCg       <  1:2861914/148‑1 (MQ=255)
                                                                                                                                          tCATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGt      >  2:2187525/1‑148 (MQ=255)
                                                                                                                                           cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt     >  2:2134895/1‑148 (MQ=255)
                                                                                                                                           cATACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGtt     >  2:2027571/1‑148 (MQ=255)
                                                                                                                                             tACCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTgg   <  2:699480/148‑1 (MQ=255)
                                                                                                                                              aCCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc  >  2:2523721/1‑148 (MQ=255)
                                                                                                                                              aCCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc  <  1:3010330/148‑1 (MQ=255)
                                                                                                                                              aCCTATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGc  <  1:1103460/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGCGCGGCGCAGTAAAGCTGTCGAAACAGAAAGTGATCGTCAAACATCTGGATGCTATTCAGAACTTTGGCGCAATGGATATTCTGTGCACTGATAAAACCGGCACCCTGACGCAGGATAAAATTGTGCTGGAGAATCATACCGATATCTCCGGTAAAACCAGCGAACGCGTGCTGCATAGCGCGTGGTTGAACAGTCATTATCAGACCGGACTTAAAAACCTGCTCGATACAGCGGTGCTCGAAGGTACGGATGAAGAGTCAGCGCGCTCGCTGGCCAGTCGTTGGC  >  CP009273/4458464‑4458753

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: