Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,933,204 C→T 5.3% D257N (GAT→AAT)  hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,933,2040CT 5.3% 955.6 / 25.3 263D257N (GAT→AAT) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base C (129/120);  new base T (6/8);  total (135/128)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 265 total aligned reads displayed.

CGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGG  >  CP009273/3933055‑3933352
                                                                                                                                                     |                                                                                                                                                    
cGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATc                                                                                                                                                      <  2:3347907/150‑1 (MQ=255)
 ggTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCt                                                                                                                                                     >  1:4654/1‑150 (MQ=255)
 ggTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCt                                                                                                                                                     >  1:1372621/1‑150 (MQ=255)
 ggTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCt                                                                                                                                                     >  1:23732/1‑150 (MQ=255)
    gCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAAt                                                                                                                                                  <  2:1492994/150‑1 (MQ=255)
        tcatACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGaaa                                                                                                                                              <  2:3410542/150‑1 (MQ=255)
            aCAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTaa                                                                                                                                          >  2:3213469/1‑150 (MQ=255)
            aCAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTaa                                                                                                                                          >  2:2266717/1‑150 (MQ=255)
             cagccatcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAAt                                                                                                                                         <  2:1556651/150‑1 (MQ=255)
             cagccagcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAAt                                                                                                                                         <  2:1805132/150‑1 (MQ=255)
              agccagcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATg                                                                                                                                        <  1:2586733/150‑1 (MQ=255)
                 cagcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGt                                                                                                                                     <  1:2004372/150‑1 (MQ=255)
                  agcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGtt                                                                                                                                    <  2:1889273/150‑1 (MQ=255)
                   gcAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGttt                                                                                                                                   <  1:573429/150‑1 (MQ=255)
                    cAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTg                                                                                                                                  <  1:4017262/150‑1 (MQ=255)
                     aaTAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTgg                                                                                                                                 <  1:3424243/150‑1 (MQ=255)
                       tAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGgcg                                                                                                                               <  1:1440456/150‑1 (MQ=255)
                        aaaCGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGgcgt                                                                                                                              <  1:1267515/150‑1 (MQ=255)
                         aaCGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGgcgtg                                                                                                                             >  1:947675/1‑150 (MQ=255)
                          aCGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGgcgtgc                                                                                                                            <  1:2031815/150‑1 (MQ=255)
                            gCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGgcgtgcgt                                                                                                                          <  2:4039878/150‑1 (MQ=255)
                               tGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATTTAATGAAATTAATGGGTTTGGCGTGcgtcgt                                                                                                                       >  1:877753/1‑150 (MQ=255)
                               tGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGcgtcgt                                                                                                                       >  2:3676208/1‑150 (MQ=255)
                               tGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGcgtcgt                                                                                                                       >  1:2494600/1‑150 (MQ=255)
                                gATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTg                                                                                                                      >  1:791884/1‑150 (MQ=255)
                                   aaCCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCAt                                                                                                                   <  2:2412100/150‑1 (MQ=255)
                                      cAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGga                                                                                                                <  1:1032088/150‑1 (MQ=255)
                                       aaaTCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGgag                                                                                                               >  1:3554014/1‑150 (MQ=255)
                                        aaTCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGgaga                                                                                                              >  1:4364477/1‑150 (MQ=255)
                                                cTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATagag                                                                                                      <  1:2124912/150‑1 (MQ=255)
                                                 tGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAgaga                                                                                                     >  1:964648/1‑150 (MQ=255)
                                                 tGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAgaga                                                                                                     <  1:2965736/150‑1 (MQ=255)
                                                     aCATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCTTGCGTCGTGCATGGAGAATATAGAGAAGCa                                                                                                 <  1:4286836/150‑1 (MQ=255)
                                                     aCATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCTTGCGTCGTGCATGGAGAATATAGAGAAGCa                                                                                                 <  1:4286954/150‑1 (MQ=255)
                                                     aCATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCa                                                                                                 >  2:700079/1‑150 (MQ=255)
                                                          aaTGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAAc                                                                                            >  1:1483485/1‑150 (MQ=255)
                                                              gCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGa                                                                                        <  2:3877065/150‑1 (MQ=255)
                                                                 tAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGc                                                                                     <  1:812725/150‑1 (MQ=255)
                                                                  aaGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCt                                                                                    >  1:1118094/1‑150 (MQ=255)
                                                                  aaGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCt                                                                                    <  2:2940180/150‑1 (MQ=255)
                                                                  aaGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATAGTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCt                                                                                    <  2:1632290/150‑1 (MQ=255)
                                                                   aGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTg                                                                                   >  1:4486040/1‑150 (MQ=255)
                                                                    ggAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTgg                                                                                  <  2:351345/150‑1 (MQ=255)
                                                                     gAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGTGAAGCAGTAACCCGATGCTGGt                                                                                 >  2:2716245/1‑150 (MQ=255)
                                                                     gAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGt                                                                                 >  2:4023203/1‑150 (MQ=255)
                                                                      aaCGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCCTGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTg                                                                                >  2:2328713/1‑150 (MQ=255)
                                                                      aaCGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTg                                                                                >  2:20978/1‑150 (MQ=255)
                                                                         ggTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACa                                                                             >  2:3173945/1‑150 (MQ=255)
                                                                              acaGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTAc                                                                        >  1:4011173/1‑150 (MQ=255)
                                                                                aGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCg                                                                      >  1:4112610/1‑150 (MQ=255)
                                                                                   cGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCa                                                                   >  1:829397/1‑150 (MQ=255)
                                                                                     gTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAg                                                                 <  1:1146732/150‑1 (MQ=255)
                                                                                      tCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAgg                                                                <  1:3117130/150‑1 (MQ=255)
                                                                                       ccccAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGc                                                               >  2:2541177/1‑150 (MQ=255)
                                                                                           aGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATTTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATc                                                           <  2:4337206/150‑1 (MQ=255)
                                                                                            gACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATcc                                                          >  1:1687973/1‑150 (MQ=255)
                                                                                             aCGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCa                                                         <  1:3246130/150‑1 (MQ=255)
                                                                                              cGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAg                                                        >  1:1597847/1‑150 (MQ=255)
                                                                                              cGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGACAATCCAg                                                        <  1:626858/150‑1 (MQ=255)
                                                                                                 cttGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTATAAATAAATCTAATGACATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAAGTACCGTCAAGGCAATCCAGCTg                                                     >  2:28163/2‑150 (MQ=255)
                                                                                                  ttGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTgg                                                    >  2:787416/1‑150 (MQ=255)
                                                                                                   tGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGc                                                   >  2:2661308/1‑150 (MQ=255)
                                                                                                   tGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGc                                                   >  1:3698147/1‑150 (MQ=255)
                                                                                                    gCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCg                                                  <  2:3804551/150‑1 (MQ=255)
                                                                                                         atTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAAt                                             <  2:4401815/150‑1 (MQ=255)
                                                                                                               tACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTc                                       >  1:4443255/1‑150 (MQ=255)
                                                                                                                aCGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTcc                                      <  2:2049826/150‑1 (MQ=255)
                                                                                                                  gATACCGATCGAGAAAGTGCGGGTTTTAAATAAATTTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTcccc                                    >  2:2500995/1‑150 (MQ=255)
                                                                                                                  gATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTcccc                                    >  2:3006608/1‑150 (MQ=255)
                                                                                                                   aTACCGATCGAGAAAGTGCGGTTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCg                                   <  1:1409911/150‑1 (MQ=255)
                                                                                                                    tACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGa                                  >  1:741767/1‑150 (MQ=255)
                                                                                                                     aCCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGaa                                 <  2:1707470/150‑1 (MQ=255)
                                                                                                                        gATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAg                              >  2:3583118/1‑150 (MQ=255)
                                                                                                                         aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGc                             <  2:4132944/150‑1 (MQ=255)
                                                                                                                         aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGc                             >  2:216980/1‑150 (MQ=255)
                                                                                                                          tCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGct                            >  2:3639546/1‑150 (MQ=255)
                                                                                                                           cGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGctc                           <  2:672289/150‑1 (MQ=255)
                                                                                                                            gagaAAGTGCGGGTTTTAAATAAATTTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCATCTGGCGACAATCTTTTCCCCGAATAGctct                          >  2:334421/1‑150 (MQ=255)
                                                                                                                             agaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTa                         <  2:1423145/150‑1 (MQ=255)
                                                                                                                             agaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTa                         <  1:224944/150‑1 (MQ=255)
                                                                                                                              gaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAt                        >  2:273689/1‑150 (MQ=255)
                                                                                                                               aaaGTGCGGGTTTTAAATAAATTTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                       >  2:116718/1‑150 (MQ=255)
                                                                                                                               aaaGTGCGGGTTTTAAATAAATTTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTATTCCCCGAATAGCTCTAtt                       >  2:1312843/1‑150 (MQ=255)
                                                                                                                               aaaGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                       <  2:1767726/150‑1 (MQ=255)
                                                                                                                               aaaGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                       >  2:3856086/1‑150 (MQ=255)
                                                                                                                                   tGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGc                   >  1:1484518/1‑150 (MQ=255)
                                                                                                                                   tGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGc                   <  2:589287/150‑1 (MQ=255)
                                                                                                                                     cGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCtt                 >  2:144263/1‑150 (MQ=255)
                                                                                                                                     cGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCtt                 >  2:4611751/1‑150 (MQ=255)
                                                                                                                                       ggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTga               <  2:924138/150‑1 (MQ=255)
                                                                                                                                           ttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaag           <  1:974052/150‑1 (MQ=255)
                                                                                                                                            taaataaaTCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaaga          >  2:1410590/1‑150 (MQ=255)
                                                                                                                                             aaataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaagag         >  2:3256482/1‑150 (MQ=255)
                                                                                                                                              aataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaagaga        <  1:2072418/150‑1 (MQ=255)
                                                                                                                                               ataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaagagaa       >  2:4187977/1‑150 (MQ=255)
                                                                                                                                                taaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAAc      <  1:2735006/150‑1 (MQ=255)
                                                                                                                                                 aaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACa     >  1:100741/1‑150 (MQ=255)
                                                                                                                                                 aaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACa     >  2:3102377/1‑150 (MQ=255)
                                                                                                                                                 aaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGAAGTAACCCGATGATGGTGACAATTACCGTCAAGGCAAACCAGCTGGCGACAATCTTTTCCACGAATAGCTCTATACCGCTTGAGAAGAGAACa     >  2:400642/1‑150 (MQ=255)
                                                                                                                                                    tCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCTATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAgg  >  1:1858547/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGTGCCATCATACAGCCAGCAATAAACGCCTGATAACCAAATCCTACCTGTAACATCAATGGCATAAGGAACGGTACACAGCCGGTCCCCAGACGGGTTGCAATATTGCCTACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGAACAAGG  >  CP009273/3933055‑3933352

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: