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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,9480GT88.7% 483.4 / 60.3 195D503Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base T (82/91);  total (93/102)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 198 total aligned reads displayed.

GGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAC  >  CP009273/4458803‑4459093
                                                                                                                                                 |                                                                                                                                                 
ggtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTAt                                                                                                                                                 <  2:2306200/148‑1 (MQ=255)
 gtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTtata                                                                                                                                                <  2:356065/148‑1 (MQ=255)
 gtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTtata                                                                                                                                                <  1:769486/148‑1 (MQ=255)
   aGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACg                                                                                                                                              <  2:2134110/148‑1 (MQ=255)
     tGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCt                                                                                                                                            >  2:1029401/1‑148 (MQ=255)
      ggCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTg                                                                                                                                           <  1:330504/148‑1 (MQ=255)
       gCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGa                                                                                                                                          >  1:2054080/1‑148 (MQ=255)
        cAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGaa                                                                                                                                         <  1:2230756/148‑1 (MQ=255)
           aaaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                      <  2:926097/148‑1 (MQ=255)
             aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTc                                                                                                                                    >  2:2087375/1‑148 (MQ=255)
             aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTc                                                                                                                                    <  1:3196134/148‑1 (MQ=255)
                 ccGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgg                                                                                                                                 <  2:2579962/147‑1 (MQ=255)
                  cGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgggg                                                                                                                               <  1:1825581/148‑1 (MQ=255)
                  cGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgggg                                                                                                                               <  2:518810/148‑1 (MQ=255)
                   gAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGc                                                                                                                              <  1:215770/148‑1 (MQ=255)
                    aGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCt                                                                                                                             >  1:1658729/1‑148 (MQ=255)
                     gCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTg                                                                                                                            <  1:1507489/148‑1 (MQ=255)
                     gCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTg                                                                                                                            <  1:2397225/148‑1 (MQ=255)
                       aCCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTgcg                                                                                                                          <  1:1177665/148‑1 (MQ=255)
                           tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTg                                                                                                                      <  1:2853047/148‑1 (MQ=255)
                           tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTg                                                                                                                      >  2:1708400/1‑148 (MQ=255)
                           tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTg                                                                                                                      >  2:1708185/1‑148 (MQ=255)
                             aGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGt                                                                                                                    <  1:494026/148‑1 (MQ=255)
                             aGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGAGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGt                                                                                                                    >  1:185617/1‑148 (MQ=255)
                                    ttGCAAAGGTGCATTTCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             <  1:493604/148‑1 (MQ=255)
                                    ttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             >  2:29040/1‑148 (MQ=255)
                                    ttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             <  2:907359/148‑1 (MQ=255)
                                       cAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcg                                                                                                          <  1:3304037/148‑1 (MQ=255)
                                            gTGCATTGCAGGAAATCCTCACTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaa                                                                                                     <  2:609207/148‑1 (MQ=255)
                                            gTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaa                                                                                                     <  2:390294/148‑1 (MQ=255)
                                            gTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaa                                                                                                     >  1:2715126/1‑148 (MQ=255)
                                             tGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaaa                                                                                                    >  2:1555525/1‑148 (MQ=255)
                                               cATTGCGGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATa                                                                                                  >  2:725435/1‑148 (MQ=255)
                                                aTTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAc                                                                                                 <  1:3003494/148‑1 (MQ=255)
                                                 ttGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAcc                                                                                                >  2:3158139/1‑148 (MQ=255)
                                                  tGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCt                                                                                               >  2:511472/1‑148 (MQ=255)
                                                     aGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGcc                                                                                            <  2:164206/148‑1 (MQ=255)
                                                      ggAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCa                                                                                           >  1:766643/1‑148 (MQ=255)
                                                         aaTCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAgcg                                                                                        >  1:744184/1‑148 (MQ=255)
                                                            ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                     >  1:2516877/1‑148 (MQ=255)
                                                            ccTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                     <  1:464110/148‑1 (MQ=255)
                                                              tCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGa                                                                                   <  1:1119008/148‑1 (MQ=255)
                                                                   gtgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGgag                                                                              >  2:3154783/1‑148 (MQ=255)
                                                                     gtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                            >  2:1051085/1‑148 (MQ=255)
                                                                         tCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTAcc                                                                        <  1:2194224/148‑1 (MQ=255)
                                                                             aGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcg                                                                    <  1:1555525/148‑1 (MQ=255)
                                                                               gTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcggg                                                                  >  1:1600456/1‑148 (MQ=255)
                                                                                    tCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGAt                                                             >  1:1299293/1‑148 (MQ=255)
                                                                                      acaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGa                                                           <  1:2457710/148‑1 (MQ=255)
                                                                                      acaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGa                                                           <  2:1050656/148‑1 (MQ=255)
                                                                                       caATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGaa                                                          >  1:1070806/1‑148 (MQ=255)
                                                                                       caATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGaa                                                          >  1:1658622/1‑148 (MQ=255)
                                                                                         aTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATc                                                        <  1:2787122/148‑1 (MQ=255)
                                                                                         aTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATc                                                        >  2:516902/1‑148 (MQ=255)
                                                                                          tGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATcc                                                       >  1:1887732/1‑148 (MQ=255)
                                                                                          tGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATcc                                                       >  2:3193859/1‑148 (MQ=255)
                                                                                            gCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGa                                                     <  1:2217155/148‑1 (MQ=255)
                                                                                             cGAGTTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGACTCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCGAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAc                                                    <  2:195342/148‑1 (MQ=255)
                                                                                             cGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAc                                                    <  1:940792/148‑1 (MQ=255)
                                                                                             cGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAc                                                    <  1:2911262/148‑1 (MQ=255)
                                                                                                gaTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTg                                                 <  2:2873726/148‑1 (MQ=255)
                                                                                                 aTTGTGCCGCTCGATGACATCATTCTGAGTAAGATTAAGAGGGTTACTTATACGCTGAATCGTCATGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAACTGCCAGCGCTTTAAGGAGATTACCAGCGGGCTGATGAATCCGACCTGa                                                <  2:987805/148‑1 (MQ=255)
                                                                                                 aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                                <  2:1907250/148‑1 (MQ=255)
                                                                                                 aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGa                                                <  1:550101/148‑1 (MQ=255)
                                                                                                  ttGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGAt                                               <  2:1917675/148‑1 (MQ=255)
                                                                                                   tgtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                              >  2:3025368/1‑148 (MQ=255)
                                                                                                   tgtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                              <  2:2170190/148‑1 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                        >  1:3300852/1‑148 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                        <  1:634366/148‑1 (MQ=255)
                                                                                                         gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                        <  2:2256552/148‑1 (MQ=255)
                                                                                                            cGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGa                                     <  1:21654/148‑1 (MQ=255)
                                                                                                                gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGATGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatcta                                 >  2:2437845/1‑148 (MQ=255)
                                                                                                                 aCATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatatat                                >  1:692957/1‑148 (MQ=255)
                                                                                                                 aCATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatatat                                >  2:305360/1‑148 (MQ=255)
                                                                                                                    tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGc                             >  2:1673067/1‑148 (MQ=255)
                                                                                                                      atGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCtt                           <  1:3154867/148‑1 (MQ=255)
                                                                                                                         ctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTc                        <  2:1406739/148‑1 (MQ=255)
                                                                                                                          tgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCt                       <  1:472756/148‑1 (MQ=255)
                                                                                                                          tgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCt                       <  1:194334/148‑1 (MQ=255)
                                                                                                                           gcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCtt                      <  2:82845/148‑1 (MQ=255)
                                                                                                                            cGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTg                     >  2:1929673/1‑148 (MQ=255)
                                                                                                                              tAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   <  1:997207/148‑1 (MQ=255)
                                                                                                                              tAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   <  2:2935525/148‑1 (MQ=255)
                                                                                                                              tAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   <  2:1382691/148‑1 (MQ=255)
                                                                                                                                aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                 >  1:3028987/1‑148 (MQ=255)
                                                                                                                                aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                 >  1:3028951/1‑148 (MQ=255)
                                                                                                                                  aTTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGc               <  1:2633739/148‑1 (MQ=255)
                                                                                                                                   ttAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGcc              >  2:2731438/1‑148 (MQ=255)
                                                                                                                                   ttAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGcc              >  1:769045/1‑148 (MQ=255)
                                                                                                                                   ttAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGcc              >  2:1270796/1‑148 (MQ=255)
                                                                                                                                    tAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa             >  2:1167983/1‑148 (MQ=255)
                                                                                                                                      aGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCaaa           >  2:2970027/1‑148 (MQ=255)
                                                                                                                                        cGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAag         <  1:2087375/148‑1 (MQ=255)
                                                                                                                                          ggTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAagag       >  2:3109128/1‑148 (MQ=255)
                                                                                                                                            ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAACAGAGAc     >  2:3119675/1‑148 (MQ=255)
                                                                                                                                            ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATAACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGAc     >  1:1400225/1‑148 (MQ=255)
                                                                                                                                             tACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACa    <  1:326875/148‑1 (MQ=255)
                                                                                                                                              aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa   >  1:2313792/1‑148 (MQ=255)
                                                                                                                                               cTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAc  <  1:1708400/148‑1 (MQ=255)
                                                                                                                                               cTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAc  <  1:1708185/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAC  >  CP009273/4458803‑4459093

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