Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,614,837 G→T 64.6% I186I (ATC→ATA eamA ← cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,614,8370GT64.6% 101.5 / 218.1 181I186I (ATC→ATAeamAcysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base G (27/37);  new base T (61/56);  total (88/93)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-01

Only 100 of 185 total aligned reads displayed.

TAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATGCT  >  CP009273/1614690‑1614982
                                                                                                                                                   |                                                                                                                                                 
tAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATg                                                                                                                                                   >  1:2869221/1‑148 (MQ=255)
tAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATg                                                                                                                                                   >  1:2111513/1‑148 (MQ=255)
  aCGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTAt                                                                                                                                                 <  2:587285/148‑1 (MQ=255)
  aCGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGAt                                                                                                                                                 >  2:1525853/1‑148 (MQ=255)
     tCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGCATTATTgg                                                                                                                                              >  2:3560991/1‑148 (MQ=255)
       ccccAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGa                                                                                                                                            <  1:1998279/148‑1 (MQ=255)
           aGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTaa                                                                                                                                        <  1:1058576/148‑1 (MQ=255)
                 cATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCt                                                                                                                                  >  1:1656907/1‑148 (MQ=255)
                 cATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCt                                                                                                                                  >  1:1656897/1‑148 (MQ=255)
                  aTAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTc                                                                                                                                 >  2:1917410/1‑148 (MQ=255)
                  aTAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTc                                                                                                                                 >  1:2298475/1‑148 (MQ=255)
                     aCCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAg                                                                                                                              <  1:2971135/148‑1 (MQ=255)
                     aCCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTACTGGGATTAAAGCGCTCCAg                                                                                                                              <  1:3563255/148‑1 (MQ=255)
                      ccAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGa                                                                                                                             <  2:3083411/148‑1 (MQ=255)
                       caacaaTTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGAt                                                                                                                            >  2:639937/1‑148 (MQ=255)
                        aacaaTTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGAtt                                                                                                                           >  1:3373163/1‑148 (MQ=255)
                          caaTTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTAc                                                                                                                         <  2:720998/148‑1 (MQ=255)
                            aTTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCa                                                                                                                       >  1:2892893/1‑148 (MQ=255)
                             ttGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAg                                                                                                                      >  1:3293956/1‑148 (MQ=255)
                                 cGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGAc                                                                                                                  <  2:3219661/148‑1 (MQ=255)
                                  gCCAAAATCGCCAGATACATCAGAGACAAGATGGTGGTCATATCGTTATTAACCAGAGTTTGAATCATGGTTGCGGAACCATCTAGAATCAGCGAGGCAACAAAGAATGCAATTTTTGGGATTAAAGCGATCCAGATTACCAGCGACa                                                                                                                 <  1:2418478/148‑1 (MQ=255)
                                   ccACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACAt                                                                                                                >  1:3517962/1‑148 (MQ=255)
                                   ccACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACAt                                                                                                                >  1:3518468/1‑148 (MQ=255)
                                   ccACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACAt                                                                                                                >  1:2432074/1‑148 (MQ=255)
                                   ccACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACAt                                                                                                                >  1:2432386/1‑148 (MQ=255)
                                   ccACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACAt                                                                                                                <  2:2760402/148‑1 (MQ=255)
                                    cacaAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATc                                                                                                               <  1:2123873/148‑1 (MQ=255)
                                    cacaAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATc                                                                                                               <  2:1332012/148‑1 (MQ=255)
                                     acaAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCa                                                                                                              <  1:1332192/148‑1 (MQ=255)
                                     acaAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCa                                                                                                              <  1:1112833/148‑1 (MQ=255)
                                          cGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCAccgcc                                                                                                         >  2:816202/1‑148 (MQ=255)
                                          cGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCAccgcc                                                                                                         >  1:862345/1‑148 (MQ=255)
                                               gATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGAcg                                                                                                    >  2:3016146/1‑148 (MQ=255)
                                               gATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGAcg                                                                                                    >  2:904484/1‑148 (MQ=255)
                                                aTACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGAcgc                                                                                                   >  1:3405226/1‑148 (MQ=255)
                                                 tACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGAcgcg                                                                                                  >  1:590086/1‑148 (MQ=255)
                                                 tACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGAcgcg                                                                                                  >  1:187698/1‑148 (MQ=255)
                                                       agagACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGt                                                                                            <  1:434149/148‑1 (MQ=255)
                                                         agaCAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGc                                                                                          >  1:3495949/1‑148 (MQ=255)
                                                         agaCAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGc                                                                                          >  1:834259/1‑148 (MQ=255)
                                                          gaCAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCg                                                                                         <  2:2694084/148‑1 (MQ=255)
                                                          gaCAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCg                                                                                         >  1:553591/1‑148 (MQ=255)
                                                           aCAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGa                                                                                        >  2:1684094/1‑148 (MQ=255)
                                                            cAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGAc                                                                                       <  2:1146416/148‑1 (MQ=255)
                                                            cAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGAc                                                                                       >  2:2883690/1‑148 (MQ=255)
                                                            cAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGAc                                                                                       <  2:800721/148‑1 (MQ=255)
                                                            cAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGAc                                                                                       <  1:627815/148‑1 (MQ=255)
                                                             aaGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACa                                                                                      <  2:76484/148‑1 (MQ=255)
                                                               gATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATg                                                                                    >  1:3143119/1‑148 (MQ=255)
                                                               gATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATg                                                                                    <  1:1566053/148‑1 (MQ=255)
                                                                aTGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGa                                                                                   >  1:1048765/1‑148 (MQ=255)
                                                                  ggtggtCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATc                                                                                 <  1:1474937/148‑1 (MQ=255)
                                                                    tggtCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCtt                                                                               <  1:1646592/148‑1 (MQ=255)
                                                                     ggtCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCttt                                                                              <  1:1036410/148‑1 (MQ=255)
                                                                     ggtCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCttt                                                                              <  1:2081097/148‑1 (MQ=255)
                                                                     ggtCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCttt                                                                              <  1:1983985/148‑1 (MQ=255)
                                                                       tCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCttttt                                                                            <  2:1656907/148‑1 (MQ=255)
                                                                       tCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCttttt                                                                            <  2:1297252/148‑1 (MQ=255)
                                                                        cATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTa                                                                           <  1:3003620/148‑1 (MQ=255)
                                                                             cGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGaa                                                                      >  2:1378313/1‑148 (MQ=255)
                                                                              gATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAg                                                                     <  1:831975/148‑1 (MQ=255)
                                                                               aTAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGa                                                                    >  1:2528697/1‑148 (MQ=255)
                                                                                  gTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGACTTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGt                                                                 >  1:1930495/1‑148 (MQ=255)
                                                                                  gTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGt                                                                 >  2:3493362/1‑148 (MQ=255)
                                                                                     aCCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGc                                                              <  2:1829975/148‑1 (MQ=255)
                                                                                        aGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCcac                                                           <  2:843504/148‑1 (MQ=255)
                                                                                          aCTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCacac                                                         >  1:3045983/1‑148 (MQ=255)
                                                                                          aCTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCacac                                                         >  2:1454391/1‑148 (MQ=255)
                                                                                             gtgAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGcc                                                      <  2:2296000/148‑1 (MQ=255)
                                                                                               gAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCa                                                    <  2:420843/148‑1 (MQ=255)
                                                                                                aaTCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCaa                                                   >  2:3183807/1‑148 (MQ=255)
                                                                                                  tCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACACAGAAGGGACTTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACt                                                 >  2:3303302/1‑148 (MQ=255)
                                                                                                  tCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACt                                                 >  1:3521756/1‑148 (MQ=255)
                                                                                                    aTGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTaa                                               <  1:1381191/148‑1 (MQ=255)
                                                                                                               aCCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTACAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTgccgcc                                    >  2:3279950/1‑148 (MQ=255)
                                                                                                                   tCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAggg                                <  1:620859/148‑1 (MQ=255)
                                                                                                                   tCGAGAATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAggg                                >  2:1604552/1‑148 (MQ=255)
                                                                                                                      agaATCAGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCa                             >  2:957019/1‑148 (MQ=255)
                                                                                                                          tCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACAt                         <  1:2345259/148‑1 (MQ=255)
                                                                                                                          tCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACAt                         <  2:3405226/148‑1 (MQ=255)
                                                                                                                           cAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATa                        >  1:2615194/1‑148 (MQ=255)
                                                                                                                           cAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATa                        >  2:361129/1‑148 (MQ=255)
                                                                                                                            aGCGAGGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATaa                       <  2:1535731/148‑1 (MQ=255)
                                                                                                                            aGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATaa                       <  1:178894/148‑1 (MQ=255)
                                                                                                                            aGCGACGCAACAAAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATaa                       <  2:2777735/148‑1 (MQ=255)
                                                                                                                                   caacaaAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAg                >  2:261511/1‑148 (MQ=255)
                                                                                                                                   caacaaAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAg                <  1:3278606/148‑1 (MQ=255)
                                                                                                                                     acaaAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCa              >  1:1142838/1‑148 (MQ=255)
                                                                                                                                      caaAGAAGGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCAt             <  2:1653844/148‑1 (MQ=255)
                                                                                                                                       aaaGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATc            >  2:500112/1‑148 (MQ=255)
                                                                                                                                          gaagGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGcc         >  1:3565018/1‑148 (MQ=255)
                                                                                                                                          gaagGGAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGcc         >  2:2136058/1‑148 (MQ=255)
                                                                                                                                             gggAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCcaca      >  1:1162159/1‑148 (MQ=255)
                                                                                                                                              ggAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACAt     >  2:914767/1‑148 (MQ=255)
                                                                                                                                               gAATTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATg    <  1:3183807/148‑1 (MQ=255)
                                                                                                                                               gAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATg    <  1:1525853/148‑1 (MQ=255)
                                                                                                                                               gAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATg    <  1:162845/148‑1 (MQ=255)
                                                                                                                                                aaTTATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATGc   <  1:957019/148‑1 (MQ=255)
                                                                                                                                                aaTGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATGc   <  2:1708038/148‑1 (MQ=255)
                                                                                                                                                 aTGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATGCt  >  1:2193742/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAACGTCCCCCAGATCCCATAACCAACAATTGTCGCCACAAACGCCAGATACATCAGAGACAAGATGGTGGTCATATCGATAGTAACCAGACTGTGAATCATGGTTGCGGAACCATCGAGAATCAGCGAGGCAACAAAGAAGGGAATGATTGGGATTAAAGCGCTCCAGATTACCAGCGACATCACCGCCGGACGCGTTGAGTGCGACATGATCTTTTTATTGAAGATGTTGCCACACGCCCAACTAAATGCTGCCGCCAGGGTCAACATAAAGCCGAGCATCGCCACATGCT  >  CP009273/1614690‑1614982

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: