Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,753,655 T→C intergenic (+256/‑525) smpB → / → intA small protein B/prophage CP4‑57 integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,753,6550TC99.2% 829.7 / ‑5.7 262intergenic (+256/‑525)smpB/intAsmall protein B/prophage CP4‑57 integrase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (117/142);  total (120/142)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 267 total aligned reads displayed.

GAGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753565‑2753744
                                                                                          |                                                                                         
gAGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGgcgc                                                                                         >  2:1510440/1‑93 (MQ=255)
 aGGCTTGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  1:4722974/93‑1 (MQ=255)
 aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTCCCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  2:4014062/93‑1 (MQ=255)
 aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  1:5334242/93‑1 (MQ=255)
 aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        >  2:773653/1‑93 (MQ=255)
  ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       >  1:1036401/1‑93 (MQ=255)
  ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       <  2:468929/93‑1 (MQ=255)
      tgTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCg                                                                                   <  2:500793/92‑1 (MQ=255)
      ggTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCg                                                                                   <  2:1004265/93‑1 (MQ=255)
       gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  >  2:4261882/1‑93 (MQ=255)
        tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 >  1:1133822/1‑93 (MQ=255)
        tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 >  1:3464797/1‑93 (MQ=255)
        tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 >  1:6608986/1‑93 (MQ=255)
        tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 <  1:1320855/93‑1 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              >  2:4224616/1‑93 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              <  1:2292721/93‑1 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              <  2:4751731/93‑1 (MQ=255)
             gcgcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCg                                                                            <  2:5454559/93‑1 (MQ=255)
               gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          <  1:773379/93‑1 (MQ=255)
                cTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGtt                                                                         <  1:3358627/93‑1 (MQ=255)
                cTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGtt                                                                         >  1:4478380/1‑93 (MQ=255)
                  cATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTgg                                                                       >  1:5911186/1‑93 (MQ=255)
                     aaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCt                                                                    >  2:526323/1‑93 (MQ=255)
                     aaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCt                                                                    <  1:4183046/93‑1 (MQ=255)
                      aaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTc                                                                   >  1:5577284/1‑93 (MQ=255)
                      aaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTc                                                                   <  2:11526/93‑1 (MQ=255)
                       aTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCg                                                                  <  1:6396302/93‑1 (MQ=255)
                          tGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCACGGCGCATGCCGAGGGGCGGTTGGCCTCGTaa                                                               >  2:6682577/1‑93 (MQ=255)
                             tataCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCTAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaaa                                                            <  2:2279991/93‑1 (MQ=255)
                               taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  1:1466177/1‑93 (MQ=255)
                               taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          <  2:5861219/93‑1 (MQ=255)
                                aCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGcc                                                         >  1:3126845/1‑93 (MQ=255)
                                 cTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCg                                                        >  2:6146191/1‑93 (MQ=255)
                                  ttACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGc                                                       >  1:4609052/1‑93 (MQ=255)
                                   tACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCa                                                      >  2:1293177/1‑93 (MQ=255)
                                   tACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCa                                                      <  1:6127319/93‑1 (MQ=255)
                                    aCCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGCAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaa                                                     >  2:734349/1‑93 (MQ=255)
                                     ccTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaa                                                    >  2:3937845/1‑93 (MQ=255)
                                       tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  >  1:2677121/1‑93 (MQ=255)
                                       tttACACATGGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  <  1:2465887/93‑1 (MQ=255)
                                         tACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCTAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAATAAAt                                                >  2:50377/1‑93 (MQ=255)
                                         tACACATTGGGGCTGATTCTGGATGCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAAt                                                <  1:154627/93‑1 (MQ=255)
                                          tcacaTTGGGGCTGATGCTGGATTCGACGGGATTTGCGAATCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               <  1:2053268/92‑1 (MQ=255)
                                          acacATTGGGTCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               <  2:459035/93‑1 (MQ=255)
                                          acacATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               <  2:4433784/93‑1 (MQ=255)
                                          acacATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               >  1:3962809/1‑93 (MQ=255)
                                            acaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGt                                             >  2:292818/1‑93 (MQ=255)
                                            acaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGt                                             >  2:1700540/1‑93 (MQ=255)
                                             caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            <  2:4513132/93‑1 (MQ=255)
                                             caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            <  1:368667/93‑1 (MQ=255)
                                             caTTGGGGCTGATTCCGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGtg                                            >  1:2190548/1‑92 (MQ=255)
                                               ttGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGc                                          >  1:1698839/1‑93 (MQ=255)
                                                tGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCa                                         >  2:2604008/1‑93 (MQ=255)
                                                 ggggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaa                                        <  2:6012889/93‑1 (MQ=255)
                                                   ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      <  1:3117464/93‑1 (MQ=255)
                                                          tCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaa                               <  2:3546198/93‑1 (MQ=255)
                                                          tCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaa                               <  1:2702042/93‑1 (MQ=255)
                                                           cTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              <  2:1990489/93‑1 (MQ=255)
                                                           cTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              >  2:2861072/1‑93 (MQ=255)
                                                            tGGATTCGCCGGCATCTGCGAAACACAAGCCGCATGCAGAGTGTCGGTTGGCCTCTTAAAAAACCGCAAAAAATAGTCGCAAACGACGAAAAc                             <  2:4100166/93‑1 (MQ=255)
                                                            tGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAAc                             >  1:6510366/1‑93 (MQ=255)
                                                            tGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAAc                             >  1:5666243/1‑93 (MQ=255)
                                                             ggATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGCAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACt                            <  1:2501275/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:2353175/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           >  2:841748/1‑93 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  2:3318412/93‑1 (MQ=255)
                                                               aTTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTAc                          <  2:934395/93‑1 (MQ=255)
                                                                  cGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCt                       <  2:3049555/93‑1 (MQ=255)
                                                                   gACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCtt                      <  2:3747579/93‑1 (MQ=255)
                                                                   gACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCtt                      <  2:4394298/93‑1 (MQ=255)
                                                                   gACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCtt                      <  1:3269026/93‑1 (MQ=255)
                                                                     cGGTATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTa                    <  1:887810/93‑1 (MQ=255)
                                                                       ggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagc                  >  2:3244430/1‑93 (MQ=255)
                                                                       ggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagc                  >  1:4869876/1‑93 (MQ=255)
                                                                        gATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagca                 <  2:2133430/93‑1 (MQ=255)
                                                                        gATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagca                 <  1:3553543/93‑1 (MQ=255)
                                                                         aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                >  2:5660300/1‑93 (MQ=255)
                                                                         aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                <  2:3504626/93‑1 (MQ=255)
                                                                          tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               >  1:3893148/1‑93 (MQ=255)
                                                                           ttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCt              <  1:4352699/93‑1 (MQ=255)
                                                                           ttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCt              <  2:651798/93‑1 (MQ=255)
                                                                             gCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTta            >  1:3388756/1‑93 (MQ=255)
                                                                              cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           >  2:2050731/1‑93 (MQ=255)
                                                                               gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  1:5849388/93‑1 (MQ=255)
                                                                               gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  2:588452/93‑1 (MQ=255)
                                                                               gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          >  2:3144695/1‑93 (MQ=255)
                                                                               gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          >  2:2396538/1‑93 (MQ=255)
                                                                                aaaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAACAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaata         >  2:6538673/1‑93 (MQ=255)
                                                                                aaaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaata         >  2:4895315/1‑93 (MQ=255)
                                                                                 aaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaataa        >  1:6638211/1‑93 (MQ=255)
                                                                                 aaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaataa        <  1:5147216/93‑1 (MQ=255)
                                                                                 aaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaataa        >  2:628363/1‑93 (MQ=255)
                                                                                   cccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAcc      >  2:5542316/1‑93 (MQ=255)
                                                                                   cccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAcc      <  2:6408851/93‑1 (MQ=255)
                                                                                    ccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTGCGCTTTAGCAGCTTAATAACCt     <  1:6051759/93‑1 (MQ=255)
                                                                                      aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGc   <  1:623418/93‑1 (MQ=255)
                                                                                      aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGc   <  1:2867945/93‑1 (MQ=255)
                                                                                      aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGc   <  1:3591539/93‑1 (MQ=255)
                                                                                       aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  >  1:5043420/1‑93 (MQ=255)
                                                                                       aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  <  1:4062/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753565‑2753744

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