Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,457,499 C→A 92.8% R20S (CGT→AGT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,457,4990CA92.8% 656.5 / 42.4 236R20S (CGT→AGT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base A (118/101);  total (126/110)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 236 total aligned reads displayed.

CCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGG  >  CP009273/4457354‑4457643
                                                                                                                                                 |                                                                                                                                                
cccTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGt                                                                                                                                                <  2:641880/148‑1 (MQ=255)
   tGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGAt                                                                                                                                             <  1:932915/148‑1 (MQ=255)
    gCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATc                                                                                                                                            >  1:249689/1‑148 (MQ=255)
      tCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCa                                                                                                                                          <  1:642982/148‑1 (MQ=255)
         gCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTg                                                                                                                                       >  1:1268185/1‑148 (MQ=255)
         gCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTg                                                                                                                                       <  2:1880385/148‑1 (MQ=255)
          cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGc                                                                                                                                      <  2:56025/148‑1 (MQ=255)
           tttATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGcc                                                                                                                                     <  1:3144852/148‑1 (MQ=255)
              attaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTgg                                                                                                                                  >  2:1745487/1‑148 (MQ=255)
                taCCTTCAGGTAAGGCTTCGCCACGCCTGAAGAAATATCTGTACTGTTTCAGACAGAGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCAAAGTGATCCATTGCCTGgcg                                                                                                                                >  2:46445/1‑148 (MQ=255)
                 aCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACAGTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGgcgc                                                                                                                               >  1:429909/1‑148 (MQ=255)
                    ttCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCgca                                                                                                                            >  1:1135671/1‑148 (MQ=255)
                        ggTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGcagcag                                                                                                                        >  1:899365/1‑148 (MQ=255)
                         gTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGa                                                                                                                       <  1:2684081/148‑1 (MQ=255)
                              gCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTg                                                                                                                  >  1:1506837/1‑148 (MQ=255)
                               cTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                 <  1:1052399/148‑1 (MQ=255)
                               cTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                 >  1:2004129/1‑148 (MQ=255)
                               cTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                 <  1:1414942/148‑1 (MQ=255)
                                ttCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGaa                                                                                                                >  1:987626/1‑148 (MQ=255)
                                 tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               <  1:2014390/148‑1 (MQ=255)
                                 tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAAt                                                                                                               >  2:3068773/1‑148 (MQ=255)
                                       cGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtg                                                                                                         >  2:502549/1‑148 (MQ=255)
                                       cGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtg                                                                                                         >  2:2864399/1‑148 (MQ=255)
                                         ccTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggtggt                                                                                                       >  1:68974/1‑148 (MQ=255)
                                              aGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgc                                                                                                  <  1:2254404/148‑1 (MQ=255)
                                                aCATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgccg                                                                                                >  1:392286/1‑148 (MQ=255)
                                                 cATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGt                                                                                               <  2:244750/148‑1 (MQ=255)
                                                  aTTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTc                                                                                              <  1:965400/148‑1 (MQ=255)
                                                  aTTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTc                                                                                              <  2:2026780/148‑1 (MQ=255)
                                                   tttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTcc                                                                                             >  1:3193211/1‑148 (MQ=255)
                                                    ttCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCt                                                                                            >  2:241898/1‑148 (MQ=255)
                                                     tCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCtt                                                                                           <  2:1912330/148‑1 (MQ=255)
                                                       tGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTaa                                                                                         >  1:1034051/1‑148 (MQ=255)
                                                       tGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTaa                                                                                         >  1:238821/1‑148 (MQ=255)
                                                        gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        <  1:1793909/148‑1 (MQ=255)
                                                        gTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        <  2:1543836/148‑1 (MQ=255)
                                                        gTACTGGTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                        <  2:3183865/148‑1 (MQ=255)
                                                          aCTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTg                                                                                      <  2:2675245/148‑1 (MQ=255)
                                                             gTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTgcgc                                                                                   <  1:2466639/148‑1 (MQ=255)
                                                              tttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                  <  1:2785991/148‑1 (MQ=255)
                                                              tttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                  <  1:211808/148‑1 (MQ=255)
                                                              tttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                  <  1:211644/148‑1 (MQ=255)
                                                               ttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAt                                                                                 >  2:1509379/1‑148 (MQ=255)
                                                                   gACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGcc                                                                             >  1:1243371/1‑148 (MQ=255)
                                                                    aCAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATGGCCt                                                                            <  2:2928910/148‑1 (MQ=255)
                                                                       gTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaa                                                                         >  1:1960330/1‑148 (MQ=255)
                                                                       gTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaa                                                                         >  2:2514982/1‑148 (MQ=255)
                                                                        tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                        >  1:175390/1‑148 (MQ=255)
                                                                        tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                        >  2:1682916/1‑148 (MQ=255)
                                                                            gAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATgg                                                                    <  1:3068773/148‑1 (MQ=255)
                                                                               ggACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCgg                                                                 >  2:2729354/1‑148 (MQ=255)
                                                                               ggACTCCTTATGTTTAAAGAAATTTTTACAGGGCTAATTCGCCATTTAACTTCCCGTCTGGTTCATAGTGATCCATTTCCTGGCGCGCAGCAGACAGTGAATACGGAGGTCCCGCCGTCCTTAAGTGCGCATTGACTTAAAATGGCgg                                                                 >  2:1108965/1‑148 (MQ=255)
                                                                                  cTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGa                                                              >  2:849320/1‑148 (MQ=255)
                                                                                  cTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGa                                                              <  1:1770727/148‑1 (MQ=255)
                                                                                   tCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTACATAGGGATCCATTGCATGGCGCGCAGCAGACAGAGAATACGGTGGTCCAGCCGTCCTTAAGTGCGCATTGCAGGAACAAGGCGGTGAt                                                             >  2:2732180/1‑148 (MQ=255)
                                                                                    ccTTATGTTTAAAGAAATTTTTACCCGGCTCATACGCCATTTACATTCCCGTCTGGTTCATAGTGATACATTTCCTGGCGCGCAGCAGACATTGAATACGGTGGTCCCGACGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATg                                                            >  1:2418850/1‑148 (MQ=255)
                                                                                     cTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGc                                                           <  2:374308/148‑1 (MQ=255)
                                                                                      ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGcc                                                          >  1:3237986/1‑148 (MQ=255)
                                                                                      ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGcc                                                          <  1:1512039/148‑1 (MQ=255)
                                                                                      ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGcc                                                          <  2:1452822/148‑1 (MQ=255)
                                                                                       tATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGccc                                                         >  1:2492305/1‑148 (MQ=255)
                                                                                         tgattAAAGAAATTTGTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                       <  2:2518465/145‑1 (MQ=255)
                                                                                         tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCga                                                       <  2:2979990/148‑1 (MQ=255)
                                                                                            ttAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCACGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaaga                                                    >  2:1307724/1‑148 (MQ=255)
                                                                                             tAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaa                                                   >  1:2706214/1‑148 (MQ=255)
                                                                                             tAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaa                                                   <  1:863919/148‑1 (MQ=255)
                                                                                               aaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaaga                                                 <  1:157111/148‑1 (MQ=255)
                                                                                                aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                >  1:3102532/1‑148 (MQ=255)
                                                                                                aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                >  2:2477198/1‑148 (MQ=255)
                                                                                                 gAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                               >  1:44102/1‑148 (MQ=255)
                                                                                                 gAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                               <  2:1554013/148‑1 (MQ=255)
                                                                                                 gAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                               >  1:881233/1‑148 (MQ=255)
                                                                                                    aTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtgt                                            <  1:3162358/148‑1 (MQ=255)
                                                                                                      ttttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACANTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTgg                                          >  1:1680914/1‑148 (MQ=255)
                                                                                                             cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGt                                   <  1:570518/148‑1 (MQ=255)
                                                                                                                cTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCg                                >  1:1931229/1‑148 (MQ=255)
                                                                                                                cTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCg                                >  1:2526890/1‑148 (MQ=255)
                                                                                                                cTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTAAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCg                                >  1:2527103/1‑148 (MQ=255)
                                                                                                                   aTTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGaca                             <  2:733625/148‑1 (MQ=255)
                                                                                                                    ttCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGACAACGTTCGacac                            >  1:2601938/1‑148 (MQ=255)
                                                                                                                    ttCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGacac                            >  2:3274335/1‑148 (MQ=255)
                                                                                                                     tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                           <  2:2406382/148‑1 (MQ=255)
                                                                                                                     tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                           <  2:1200376/148‑1 (MQ=255)
                                                                                                                     tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAcc                           >  2:1492057/1‑148 (MQ=255)
                                                                                                                        ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                        >  2:2301740/1‑148 (MQ=255)
                                                                                                                        ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                        >  1:1443120/1‑148 (MQ=255)
                                                                                                                        ccATTTACCTTCCCGTCTGGTTCATAGTAATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCAt                        >  1:1656462/1‑148 (MQ=255)
                                                                                                                         cATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATc                       <  1:422500/148‑1 (MQ=255)
                                                                                                                             tACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGa                   <  1:325508/148‑1 (MQ=255)
                                                                                                                                 ttCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg               <  1:680867/148‑1 (MQ=255)
                                                                                                                                 ttCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg               <  2:68732/148‑1 (MQ=255)
                                                                                                                                    ccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa            >  1:1108732/1‑148 (MQ=255)
                                                                                                                                     cGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTaa           >  2:2193515/1‑148 (MQ=255)
                                                                                                                                      gTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTaaa          <  2:1955134/148‑1 (MQ=255)
                                                                                                                                      gTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTaaa          <  2:823082/148‑1 (MQ=255)
                                                                                                                                       tCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAAt         <  2:429909/148‑1 (MQ=255)
                                                                                                                                        cTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATc        <  1:1884559/148‑1 (MQ=255)
                                                                                                                                        cTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCAGGAAGGGTTAAATc        >  1:1417755/1‑148 (MQ=255)
                                                                                                                                          ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAg      >  1:3277712/1‑148 (MQ=255)
                                                                                                                                          ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAg      >  1:1843444/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGG  >  CP009273/4457354‑4457643

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: