Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,350,508 T→A 9.7% K11I (AAA→ATA) ‡ PP_1176 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,350,5080TA 9.7% 233.8 / 12.6 81K11I (AAA→ATA) ‡PP_1176hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base A (4/4);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GGAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTA  >  NC_002947/1350359‑1350651
                                                                                                                                                     |                                                                                                                                               
ggAGATCCGTTTCTGACGATAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAtt                                                                                                                                                 <  2:233813/150‑1 (MQ=255)
   gagCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCACAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgg                                                                                                                                              <  2:826559/150‑1 (MQ=255)
   gagCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgg                                                                                                                                              <  1:1035166/150‑1 (MQ=255)
   gagCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgg                                                                                                                                              >  2:538972/1‑150 (MQ=255)
      ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTg                                                                                                                                               <  1:607920/146‑1 (MQ=255)
      ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTg                                                                                                                                               >  2:607920/1‑146 (MQ=255)
      ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAg                                                                                                                                           >  2:1262768/1‑150 (MQ=255)
      ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAg                                                                                                                                           >  1:549716/1‑150 (MQ=255)
      ccGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAg                                                                                                                                           >  1:1511338/1‑150 (MQ=255)
        gTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGcc                                                                                                                                         >  2:159525/1‑150 (MQ=255)
          ttCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATACAACGGATTTTAAATGGCAGCCtt                                                                                                                                       <  2:1019132/150‑1 (MQ=255)
            cTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgt                                                                                                                                     <  1:1905477/150‑1 (MQ=255)
                cGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGAt                                                                                                                                 >  1:1502830/1‑150 (MQ=255)
                 gagaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgg                                                                                                                                              <  1:155814/136‑1 (MQ=255)
                 gagaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTgg                                                                                                                                              >  2:155814/1‑136 (MQ=255)
                    aaCGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttg                                                                                                                                      <  2:693563/141‑1 (MQ=255)
                    aaCGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttg                                                                                                                                      >  1:693563/1‑141 (MQ=255)
                    aaCGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTc                                                                                                                             <  1:538972/150‑1 (MQ=255)
                       gggATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGc                                                                                                                          >  2:463577/1‑150 (MQ=255)
                         gATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGccc                                                                                                                        <  1:1903410/150‑1 (MQ=255)
                            tCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGt                                                                                                                     <  1:716988/150‑1 (MQ=255)
                            tCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGt                                                                                                                     <  1:713587/150‑1 (MQ=255)
                               gCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaa                                                                                                                  >  2:277756/1‑150 (MQ=255)
                                cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                 >  2:132031/1‑150 (MQ=255)
                                cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGACAGTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                 <  1:1401218/150‑1 (MQ=255)
                                    gtgAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTcaatcaat                                                                                                             <  2:308324/150‑1 (MQ=255)
                                    gtgAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTcaatcaat                                                                                                             <  1:132031/150‑1 (MQ=255)
                                      gAGCGTCCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATcc                                                                                                           <  2:380389/150‑1 (MQ=255)
                                        gtgctCTTAAGCCACGAAAAACGTTAACGCGATGTCCAAAAACAGGCAACCGCGGCAGCCGTAGGAAACGTCCTGTTCTGCTCCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCa                                                                                                         <  2:576310/146‑1 (MQ=255)
                                          gTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGc                                                                                                       <  1:1262768/150‑1 (MQ=255)
                                           ttCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCt                                                                                                      >  1:1431315/1‑150 (MQ=255)
                                                  gcgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGAAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGcc                                                                                               >  2:31963/1‑150 (MQ=255)
                                                   cggcGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGccc                                                                                              <  1:325900/147‑1 (MQ=255)
                                                      cGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAAAGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTAAATCAATCCCAGCTCAACGCCCCAc                                                                                           >  1:793987/1‑150 (MQ=255)
                                                        gAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCAccc                                                                                         <  1:122941/150‑1 (MQ=255)
                                                         aaaaCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTTAGACAATCCCAGCTCAACGCCCCACCCg                                                                                        >  1:941218/1‑150 (MQ=255)
                                                              gTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaatcaa                                                                                                              <  2:1809814/123‑1 (MQ=255)
                                                              gTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaatcaa                                                                                                              >  1:1809814/1‑123 (MQ=255)
                                                                tATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATa                                                                                 <  2:549716/150‑1 (MQ=255)
                                                                tATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATa                                                                                 <  2:1511338/150‑1 (MQ=255)
                                                                      gATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtttttatt                                                                                                                                 <  1:472077/96‑4 (MQ=38)
                                                                      gATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtttttatt                                                                                                                                 >  2:472077/1‑93 (MQ=38)
                                                                          tCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                       <  1:159525/150‑1 (MQ=255)
                                                                            cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                     >  2:1605643/1‑150 (MQ=255)
                                                                             caaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                    >  2:1370660/2‑150 (MQ=255)
                                                                               aaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGGCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCGATCAATCCCAGCTCAACGCCACACCCGTCTGATACTCGATCACAcgctt                                                                  >  2:981507/1‑148 (MQ=255)
                                                                                          gTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTc                                                                              >  1:189951/1‑127 (MQ=255)
                                                                                          gTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTc                                                                              <  2:189951/127‑1 (MQ=255)
                                                                                                  ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                         >  2:767569/1‑140 (MQ=255)
                                                                                                  ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                         <  1:767569/140‑1 (MQ=255)
                                                                                                   cGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCct                                              >  1:420556/1‑150 (MQ=255)
                                                                                                      aGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGc                                                                                                                               >  1:1696139/1‑66 (MQ=255)
                                                                                                      aGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGc                                                                                                                               <  2:1696139/66‑1 (MQ=255)
                                                                                                       ggAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCcttctt                                          >  1:1733571/1‑150 (MQ=255)
                                                                                                             gTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTc                                    <  2:420556/150‑1 (MQ=255)
                                                                                                                       gcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCg                          <  2:1431315/150‑1 (MQ=255)
                                                                                                                        cttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGc                         >  1:1541791/1‑150 (MQ=255)
                                                                                                                        cttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGc                         <  1:315291/150‑1 (MQ=255)
                                                                                                                         ttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                <  1:176148/126‑1 (MQ=255)
                                                                                                                         ttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTtctc                                                >  2:176148/1‑126 (MQ=255)
                                                                                                                             gcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATc                    <  1:588208/150‑1 (MQ=255)
                                                                                                                              ccAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATcc                   >  1:567769/1‑150 (MQ=255)
                                                                                                                              ccAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATcc                   >  1:567552/1‑150 (MQ=255)
                                                                                                                               cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCAt                                 <  1:31532/135‑1 (MQ=255)
                                                                                                                               cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCAt                                 >  2:31532/1‑135 (MQ=255)
                                                                                                                                 ggCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                <  2:567552/150‑1 (MQ=255)
                                                                                                                                 ggCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                <  2:567769/150‑1 (MQ=255)
                                                                                                                                   cAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATggg              >  2:362982/1‑150 (MQ=255)
                                                                                                                                    aaTAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGa             >  1:1403506/1‑150 (MQ=255)
                                                                                                                                       aaaacGGTTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACg          <  2:941218/150‑1 (MQ=255)
                                                                                                                                        acaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGa                                                                            >  1:544792/3‑83 (MQ=255)
                                                                                                                                        acaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGa                                                                            <  2:544792/81‑1 (MQ=39)
                                                                                                                                        acaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACgg         >  1:1833302/3‑150 (MQ=255)
                                                                                                                                         aaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGtc                                                                 >  2:1714555/1‑93 (MQ=255)
                                                                                                                                         aaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGtc                                                                 <  1:1714555/93‑1 (MQ=255)
                                                                                                                                          aaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGAt                                                                           >  2:173682/1‑82 (MQ=255)
                                                                                                                                          aaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGAt                                                                           <  1:173682/82‑1 (MQ=255)
                                                                                                                                           aGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGtt      >  1:718887/1‑150 (MQ=255)
                                                                                                                                              cTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCCACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                         >  1:1078810/1‑96 (MQ=255)
                                                                                                                                              cTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                         <  2:1078810/96‑1 (MQ=255)
                                                                                                                                               tGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTa  <  2:1502830/150‑1 (MQ=255)
                                                                                                                                               tGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGGTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTa  >  2:1034536/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GGAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTA  >  NC_002947/1350359‑1350651

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: