Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,581,646 G→A 9.8% S275S (TCG→TCA PP_2262 → sugar ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,581,6460GA 9.8% 217.3 / 11.0 81S275S (TCG→TCAPP_2262sugar ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (32/41);  new base A (5/3);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ACCGAACCGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCAACGCCGAGATCCT  >  NC_002947/2581497‑2581795
                                                                                                                                                     |                                                                                                                                                       
acCGAACCGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCg                                                                                                                                                         >  2:13852/1‑150 (MQ=255)
 ccGAACCGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACACTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTACTGATCATCCTGCTGGTGTCGt                                                                                                                                                        <  2:95740/150‑1 (MQ=255)
   gAACCGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTgg                                                                                                                                                      <  1:1186094/150‑1 (MQ=255)
   gAACCGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTgg                                                                                                                                                      <  2:1023406/150‑1 (MQ=255)
       cGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGTGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCATGGCTGt                                                                                                                                                  >  2:385884/1‑150 (MQ=255)
         ttCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTc                                                                                                                                                >  1:1058122/1‑150 (MQ=255)
          tCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCt                                                                                                                                               <  2:1963872/150‑1 (MQ=255)
           cGTGCTGACCGGCCGCGGGCCATGCAATTCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTa                                                                                                                                              <  1:111720/150‑1 (MQ=255)
                tGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAg                                                                                                                                       <  1:1894045/150‑1 (MQ=255)
                    cggcggcggGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCAt                                                                                                                                   <  2:2258402/150‑1 (MQ=255)
                    cggcggcggGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCAt                                                                                                                                   <  1:1229188/150‑1 (MQ=255)
                      gcggcggGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGa                                                                                                                                 <  1:13852/150‑1 (MQ=255)
                      gcggcggGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGa                                                                                                                                 <  2:1576223/150‑1 (MQ=255)
                       cggcggGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTAGTGGCTGTTCTACACAGCCATGAc                                                                                                                                <  1:420955/150‑1 (MQ=255)
                       atctcggcgttgctcttgaacgacatgttcagcaaccagtagatcggcactagcaggaagaagaagtacagcagcagggccattgatttgcgtttgctcatggcctagtccttgtcggcatagGTCATGGCTGTTCTACACAGCCATGAc                                                                                                                                >  2:1253846‑M1/124‑150 (MQ=255)
                               cAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCGAGAGTGTCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCg                                                                                                                        >  1:1329374/1‑150 (MQ=255)
                               cAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCAt                                                                                                                                   >  2:2358985/1‑139 (MQ=255)
                               cAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCAt                                                                                                                                   <  1:2358985/139‑1 (MQ=255)
                                 ggCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGTGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGAc                                                                                                                      >  2:2122777/1‑150 (MQ=255)
                                    aaTGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAg                                                                                                                   >  1:423856/1‑150 (MQ=255)
                                       gCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGAc                                                                                                                <  2:2202718/150‑1 (MQ=255)
                                          aCCACGTTCCTCAGCAAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCTATGCCGACAAGGACTAg                                                                                                             <  2:2041496/150‑1 (MQ=255)
                                            cgcGTTCCTCAGCCAGTTCGTCGCGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTACTGATCTTCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGc                                                                                                           <  2:802392/148‑1 (MQ=255)
                                            cACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGt                                                                                                                                                        <  1:1957868/107‑1 (MQ=255)
                                            cACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGt                                                                                                                                                        >  2:1957868/1‑107 (MQ=255)
                                            cACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCt                                                                                                                                               >  1:1574813/1‑116 (MQ=255)
                                            cACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCt                                                                                                                                               <  2:1574813/116‑1 (MQ=255)
                                            cACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGc                                                                                                           >  1:1832463/1‑150 (MQ=255)
                                               gTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCAt                                                                                                        <  2:1783070/150‑1 (MQ=255)
                                               gTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCAt                                                                                                        >  1:1293893/1‑150 (MQ=255)
                                               gTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCAt                                                                                                        <  1:140116/150‑1 (MQ=255)
                                                   tcagcaaccagtagatcggcaccagcaggaagaagaagtacagcagcagggccattgatttgcgtgtgctcatggcctagtccttgtcggcatggGTCGTGGCTGTTCTACACAGCCATGa                                                                                                                                 <  1:404433‑M1/26‑1 (MQ=255)
                                                   tcagcaaccagtagatcggcaccagcaggaagaagaagtacagcagcagggccattgatttgcgtgtgctcatggcctagtccttgtcggcatggGTCGTGGCTGTTCTACACAGCCATGa                                                                                                                                 >  2:404433‑M1/96‑121 (MQ=255)
                                                        ccAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGtt                                                                                                                                                 >  2:1521459/1‑102 (MQ=255)
                                                        ccAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGtt                                                                                                                                                 <  1:1521459/102‑1 (MQ=255)
                                                        ccAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTaca                                                                                                                                           >  2:2448343/1‑107 (MQ=255)
                                                        ccAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTaca                                                                                                                                           <  1:2448343/107‑1 (MQ=255)
                                                         cAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGc                                                                                              >  2:290268/1‑150 (MQ=255)
                                                          aGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGAAAAGAGGTGAGGTGCCGGGGCTGTTCTACACAGCCATGACCCATGACGCCATGTACTTGGCCATGAGGACAAGCa                                                                                             >  1:644544/1‑150 (MQ=255)
                                                           gACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCaa                                                                                            >  2:1013028/1‑150 (MQ=255)
                                                           cagtagatcggcaccagcaggaagaagaagtacagcagcagggccattgatttgcgtgtgctcatggcctagtccttgtcggcatggGTCGTGGCTGTTCTACACAGCCATGAc                                                                                                                                >  1:1627463‑M1/88‑114 (MQ=255)
                                                           cagtagatcggcaccagcaggaagaagaagtacagcagcagggccattgatttgcgtgtgctcatggcctagtccttgtcggcatggGTCGTGGCTGTTCTACACAGCCATGAc                                                                                                                                <  2:1627463‑M1/27‑1 (MQ=255)
                                                             cccTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAAt                                                                                          <  2:2218404/150‑1 (MQ=255)
                                                                tCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCaa                                                                                       >  1:440746/1‑150 (MQ=255)
                                                                  aCGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATg                                                                                     <  2:2115306/150‑1 (MQ=255)
                                                                     cgTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGcc                                                                                  <  1:2360687/150‑1 (MQ=255)
                                                                              ggaagaagaagtacagcagcagggccattgatttgcgtgtgctcatgacctagtccttgtcggcatggGTCATGGCTGTTCTACACAGCCATGa                                                                                                                                 <  1:1586171‑M1/26‑1 (MQ=255)
                                                                              ggaagaagaagtacagcagcagggccattgatttgcgtgtgctcatgacctagtccttgtcggcatggGTCATGGCTGTTCTACACAGCCATGa                                                                                                                                 >  2:1586171‑M1/69‑94 (MQ=255)
                                                                                 gggCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTAc                                                                      <  2:1832463/150‑1 (MQ=255)
                                                                                  ggCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTgg                                                                                                                                                      <  1:1187791/71‑1 (MQ=255)
                                                                                  ggCAGTTCGACCTGGGGCCGGCGCCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTgg                                                                                                                                                      >  2:1187791/1‑71 (MQ=255)
                                                                                   gCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTActt                                                                    >  1:829533/1‑150 (MQ=255)
                                                                                       ttCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTActtcttc                                                                >  1:1511876/1‑150 (MQ=255)
                                                                                       ttCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTActtcttc                                                                >  2:2263932/1‑150 (MQ=255)
                                                                                       ttCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTActtcttc                                                                <  2:423856/150‑1 (MQ=255)
                                                                                                gggCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCctgctg                                                       >  1:847399/1‑150 (MQ=255)
                                                                                                            gCGTTTTCACTGCTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCCTGCCGACAAGTACTCGGCTATGCGCACACGCCAATCACTGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTAc                                           <  2:1329374/150‑1 (MQ=255)
                                                                                                              gTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGCCAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTg                                         <  1:290268/150‑1 (MQ=255)
                                                                                                              gTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTg                                         <  1:2442722/150‑1 (MQ=255)
                                                                                                               ttttCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTgg                                        >  2:1416019/1‑150 (MQ=255)
                                                                                                                 ttCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGtt                                      >  2:371908/1‑150 (MQ=255)
                                                                                                                 ttCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACAAGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGtt                                      <  2:644544/150‑1 (MQ=255)
                                                                                                                 ttCACTGGTGTACTTCCTGATCATCCTGCTCGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGtt                                      <  1:1938273/150‑1 (MQ=255)
                                                                                                                      tGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGa                                 >  2:2283155/1‑150 (MQ=255)
                                                                                                                       ggTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGaa                                >  2:2131044/1‑150 (MQ=255)
                                                                                                                          gtACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACttcttcttc                                                             >  2:2032631/1‑118 (MQ=255)
                                                                                                                          gtACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACttcttcttc                                                             <  1:2032631/118‑1 (MQ=255)
                                                                                                                              ttCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCg                         >  1:1726080/1‑150 (MQ=255)
                                                                                                                                tagtccttgtcggcatggGTCATGGCTGTGTAGAACAGCCACGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGtt                       <  2:2031668‑M1/132‑1 (MQ=255)
                                                                                                                                   tccttgtcggcatggGTCATGGCTGTGTAGAACAGCCACGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCaa                    >  1:2031668‑M1/16‑150 (MQ=255)
                                                                                                                                   gATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACCCGCAAAGCAATGGCCCTGCTGCTGTACTTCTTCATCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCaa                    >  1:1114632/1‑150 (MQ=255)
                                                                                                                                    atcatcCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAATGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAag                   >  2:2179248/1‑150 (MQ=255)
                                                                                                                                     tcatcCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAaga                  <  1:371908/150‑1 (MQ=255)
                                                                                                                                     tcatcCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCAATGAGCACACGCCGGTCAATGGCCCTGCTTCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAaga                  <  1:2122777/150‑1 (MQ=255)
                                                                                                                                      ttgtcggcatggGTCATGGCTGTGTAGAACAGCCACGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTgg                                        >  1:999337‑M1/13‑127 (MQ=255)
                                                                                                                                      ttgtcggcatggGTCATGGCTGTGTAGAACAGCCACGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTgg                                        <  2:999337‑M1/115‑1 (MQ=255)
                                                                                                                                        tcCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCa               >  1:527971/1‑150 (MQ=255)
                                                                                                                                              tgGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGa                                                                                                      <  2:1869908/57‑1 (MQ=255)
                                                                                                                                              tgGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGa                                                                                                      >  1:1869908/1‑57 (MQ=255)
                                                                                                                                               ggTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCAACGNCga        <  2:1511876/150‑1 (MQ=255)
                                                                                                                                               ggTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCAACGCCga        <  2:440746/150‑1 (MQ=255)
                                                                                                                                                 tgtCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCAACGCCgaga      <  2:527971/150‑1 (MQ=255)
                                                                                                                                                     gTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCAACGCCGAGATCCt  >  2:808247/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ACCGAACCGTTCGTGCTGACCGGCGGCGGGCCAGGCAATGCGACCACGTTCCTCAGCCAGACCCTCACGCGTATGGCTGTAGGGCAGTTCGACCTGGGGCCGGCGGCAGCGTTTTCACTGGTGTACTTCCTGATCATCCTGCTGGTGTCGTGGCTGTTCTACACAGCCATGACCCATGCCGACAAGGACTAGGCCATGAGCACACGCAAATCAATGGCCCTGCTGCTGTACTTCTTCTTCCTGCTGGTGCCGATCTACTGGTTGCTGAACATGTCGTTCAAGAGCAACGCCGAGATCCT  >  NC_002947/2581497‑2581795

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: