Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,321,863 A→T 12.1% intergenic (+42/‑70) proP → / → pmrR proline/glycine betaine transporter/putative membrane‑bound BasS regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,321,8630AT12.1% 543.6 / 55.0 189intergenic (+42/‑70)proP/pmrRproline/glycine betaine transporter/putative membrane‑bound BasS regulator
Reads supporting (aligned to +/- strand):  ref base A (83/83);  new base T (14/9);  total (97/92)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 193 total aligned reads displayed.

TACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT  >  CP009273/4321717‑4322010
                                                                                                                                                  |                                                                                                                                                   
tACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCaa                                                                                                                                                    <  1:585554/148‑1 (MQ=255)
tACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCaa                                                                                                                                                    >  1:1162237/1‑148 (MQ=255)
  cGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAATCGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAAtt                                                                                                                                                  >  2:602040/1‑148 (MQ=255)
  cGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAAtt                                                                                                                                                  >  2:3239951/1‑148 (MQ=255)
     taataTCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGaa                                                                                                                                               >  1:3166534/1‑148 (MQ=255)
     taataTCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGaa                                                                                                                                               >  1:2831211/1‑148 (MQ=255)
     taataTCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGaa                                                                                                                                               >  1:2160090/1‑148 (MQ=255)
       atatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAct                                                                                                                                             <  1:2939771/148‑1 (MQ=255)
       atatCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAct                                                                                                                                             >  1:2746461/1‑148 (MQ=255)
           cGAGCAGAAACGCGCTGCTAGTGACCACGAGATTGCCGATTTGCAGGCGAAACGGACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACCGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACGCTTa                                                                                                                                         <  2:971035/148‑1 (MQ=255)
           cGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGAACTCTTa                                                                                                                                         <  1:3271247/148‑1 (MQ=255)
           cGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTa                                                                                                                                         <  2:3397742/148‑1 (MQ=255)
              gCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                                                                                                      >  1:844614/1‑148 (MQ=255)
              gCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                                                                                                      >  2:743638/1‑148 (MQ=255)
              gCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                                                                                                      <  2:123651/148‑1 (MQ=255)
              gCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                                                                                                      >  2:2526001/1‑148 (MQ=255)
              gCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAgg                                                                                                                                      >  2:2018869/1‑148 (MQ=255)
               cAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGt                                                                                                                                     <  2:3452992/148‑1 (MQ=255)
                 gAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTc                                                                                                                                   >  2:3262598/1‑148 (MQ=255)
                  aaaaTCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCa                                                                                                                                  >  1:3493056/1‑148 (MQ=255)
                      tCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTa                                                                                                                              <  2:3346209/148‑1 (MQ=255)
                      tCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTa                                                                                                                              <  2:2210585/148‑1 (MQ=255)
                      tCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTa                                                                                                                              <  2:419598/148‑1 (MQ=255)
                      tCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTa                                                                                                                              >  2:711252/1‑148 (MQ=255)
                       cGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTaa                                                                                                                             >  1:1796933/1‑148 (MQ=255)
                       cGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTaa                                                                                                                             >  1:1797244/1‑148 (MQ=255)
                        gatgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAt                                                                                                                            <  2:1110422/148‑1 (MQ=255)
                         atgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTAAtc                                                                                                                           >  2:3095000/1‑148 (MQ=255)
                          tgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtct                                                                                                                          <  1:3133518/148‑1 (MQ=255)
                          tgatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtct                                                                                                                          <  1:1364750/148‑1 (MQ=255)
                           gatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtctc                                                                                                                         >  2:2418647/1‑148 (MQ=255)
                           gatATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAAtctc                                                                                                                         <  2:2344791/148‑1 (MQ=255)
                            atatTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATctct                                                                                                                        >  1:3036913/1‑148 (MQ=255)
                              atTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGa                                                                                                                      >  2:1065521/1‑148 (MQ=255)
                                   ccACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACgcgc                                                                                                                 >  2:2851932/1‑148 (MQ=255)
                                   ccACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACgcgc                                                                                                                 >  2:942003/1‑148 (MQ=255)
                                        agaTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAct                                                                                                            >  2:211727/1‑148 (MQ=255)
                                          aTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATActct                                                                                                          >  2:1707616/1‑148 (MQ=255)
                                          aTTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATActct                                                                                                          <  2:1153986/148‑1 (MQ=255)
                                            tGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcc                                                                                                        >  2:1200886/1‑148 (MQ=255)
                                             gCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcct                                                                                                       >  1:2570579/1‑148 (MQ=255)
                                              ccGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcctc                                                                                                      <  2:79990/148‑1 (MQ=255)
                                                gATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCa                                                                                                    <  1:2803857/148‑1 (MQ=255)
                                                  tttGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAgg                                                                                                  >  1:1025872/1‑148 (MQ=255)
                                                  tttGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAgg                                                                                                  <  1:685856/148‑1 (MQ=255)
                                                   ttGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGt                                                                                                 <  2:3421872/148‑1 (MQ=255)
                                                     gCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTa                                                                                               <  2:2953579/148‑1 (MQ=255)
                                                         gCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACgg                                                                                           >  1:1232696/1‑148 (MQ=255)
                                                           gAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggag                                                                                         <  2:2861257/148‑1 (MQ=255)
                                                             aaCGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagga                                                                                       >  1:1581990/1‑148 (MQ=255)
                                                                 tACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGt                                                                                   >  2:1900448/1‑148 (MQ=255)
                                                                  aCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTg                                                                                  >  1:2334830/1‑148 (MQ=255)
                                                                   cccGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGc                                                                                 <  1:1509220/148‑1 (MQ=255)
                                                                   cccGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGc                                                                                 >  1:3290313/1‑148 (MQ=255)
                                                                    ccGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCa                                                                                <  1:2679428/148‑1 (MQ=255)
                                                                    ccGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCa                                                                                <  1:1079828/148‑1 (MQ=255)
                                                                     cGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCaa                                                                               >  2:3535903/1‑148 (MQ=255)
                                                                          ggTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaaa                                                                          >  2:323116/1‑148 (MQ=255)
                                                                          ggTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaaa                                                                          >  2:2641566/1‑148 (MQ=255)
                                                                            tGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaaaaa                                                                        >  2:2616856/1‑148 (MQ=255)
                                                                              cagcaACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAAcc                                                                      >  1:3009002/1‑148 (MQ=255)
                                                                                gcaACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGGTAACGGAGGAGAGTGCAATGAAAAACCgt                                                                    >  2:3478914/1‑148 (MQ=255)
                                                                                 caACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtg                                                                   >  2:809601/1‑148 (MQ=255)
                                                                                  aaCATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtgt                                                                  <  1:1650276/148‑1 (MQ=255)
                                                                                  aaCATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtgt                                                                  <  1:1425162/148‑1 (MQ=255)
                                                                                    cATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGttt                                                                >  1:876022/1‑148 (MQ=255)
                                                                                     aTCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTa                                                               <  1:3414667/148‑1 (MQ=255)
                                                                                         gcgAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGaa                                                           <  2:1186322/148‑1 (MQ=255)
                                                                                          cgAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGaaa                                                          <  2:1263173/148‑1 (MQ=255)
                                                                                            aaTTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGt                                                        <  1:3267108/148‑1 (MQ=255)
                                                                                               tgatgaATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTa                                                     >  1:1019164/1‑148 (MQ=255)
                                                                                                gatgaATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTaa                                                    <  2:2085246/148‑1 (MQ=255)
                                                                                                   gaATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTa                                                 >  1:1554485/1‑148 (MQ=255)
                                                                                                   gaATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTa                                                 >  1:1858644/1‑148 (MQ=255)
                                                                                                   gaATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCCCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTa                                                 <  2:1358949/148‑1 (MQ=255)
                                                                                                     aTAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTAcc                                               >  1:2164640/1‑148 (MQ=255)
                                                                                                        aGCTGAGACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCgtg                                            <  2:299944/148‑1 (MQ=255)
                                                                                                          cTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGtt                                          <  2:2905983/148‑1 (MQ=255)
                                                                                                           tGAAACGGAGGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTc                                         <  2:65024/148‑1 (MQ=255)
                                                                                                             aaaCGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAg                                       >  1:3202893/1‑148 (MQ=255)
                                                                                                                  gATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGc                                  >  1:355632/1‑148 (MQ=255)
                                                                                                                             gtgACGCTGTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTcagc                       <  2:1338180/148‑1 (MQ=255)
                                                                                                                             gtgACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTcagc                       <  2:3022670/148‑1 (MQ=255)
                                                                                                                             gtgACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTcagc                       <  1:565293/148‑1 (MQ=255)
                                                                                                                             gtgACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTcagc                       >  1:982288/1‑148 (MQ=255)
                                                                                                                             gtgACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACGACCGTGTTCAGCGTGCTGGTGGTcagc                       <  2:2450026/148‑1 (MQ=255)
                                                                                                                              tgACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTcagca                      >  1:1399963/1‑148 (MQ=255)
                                                                                                                                  gCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGctt                  >  2:1851287/1‑148 (MQ=255)
                                                                                                                                    tGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGctttc                <  1:2370547/148‑1 (MQ=255)
                                                                                                                                     gTCTTATCAGGCCTATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGctttct               >  2:3162172/1‑148 (MQ=255)
                                                                                                                                      tCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGctttctt              >  1:412743/1‑148 (MQ=255)
                                                                                                                                           tCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTatat         <  1:1384363/148‑1 (MQ=255)
                                                                                                                                            cAGGCCTATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATc        <  1:1889552/148‑1 (MQ=255)
                                                                                                                                             aGGCCTATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCt       <  2:28085/148‑1 (MQ=255)
                                                                                                                                              ggCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTg      >  1:2929796/1‑148 (MQ=255)
                                                                                                                                               gCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTgg     >  2:581594/1‑148 (MQ=255)
                                                                                                                                                ccAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGt    <  2:2831211/148‑1 (MQ=255)
                                                                                                                                                ccAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGt    <  1:1079272/148‑1 (MQ=255)
                                                                                                                                                ccAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGt    <  2:1440750/148‑1 (MQ=255)
                                                                                                                                                  aaTTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGttt  >  1:1991961/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTTTATGAAAGTTTAACTACCGTGTTCAGCGTGCTGGTGGTCAGCAGCTTTCTTTATATCTGGTTT  >  CP009273/4321717‑4322010

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: