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.7% 968.4 / ‑5.6 308intergenic (+256/‑525)smpB/intAsmall protein B/prophage CP4‑57 integrase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (138/169);  total (139/169)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 311 total aligned reads displayed.

AGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753566‑2753744
                                                                                         |                                                                                         
aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  2:2952443/93‑1 (MQ=255)
 ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       <  1:6871973/93‑1 (MQ=255)
 ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       >  2:6681448/1‑93 (MQ=255)
   gtGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGc                                                                                     <  1:3122451/92‑1 (MQ=255)
     ggTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCg                                                                                   >  2:4581837/1‑93 (MQ=255)
     ggTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCg                                                                                   >  2:5580334/1‑93 (MQ=255)
      gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  <  1:4459998/93‑1 (MQ=255)
      gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  <  1:2352532/93‑1 (MQ=255)
      gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  <  2:6186735/93‑1 (MQ=255)
       tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGAGGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 >  1:2723374/1‑93 (MQ=255)
       tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 <  1:1474116/93‑1 (MQ=255)
       tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 <  2:69396/93‑1 (MQ=255)
         aTGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACTCAAGGCGCATGCCGAggg                                                                               <  1:1219631/93‑1 (MQ=255)
          tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              <  1:4752711/93‑1 (MQ=255)
          tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              >  1:3345833/1‑93 (MQ=255)
           ggCGCTCATAAATCTGGTATACTTACCTTTACGCATTGGGGCTGACTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGc                                                                             <  1:4022613/93‑1 (MQ=255)
           ggCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGc                                                                             <  2:153675/93‑1 (MQ=255)
            gcgcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCg                                                                            <  2:2367616/93‑1 (MQ=255)
              gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          <  1:3250657/93‑1 (MQ=255)
              gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          <  1:2440471/93‑1 (MQ=255)
                tCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTg                                                                        >  1:4359320/1‑93 (MQ=255)
                   tcaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGcc                                                                     <  1:6148833/91‑1 (MQ=255)
                    aaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCt                                                                    >  2:2438180/1‑93 (MQ=255)
                     aaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTc                                                                   >  1:3906687/1‑93 (MQ=255)
                       tCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGt                                                                 >  1:7391615/1‑93 (MQ=255)
                        cTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTa                                                                <  2:6523521/93‑1 (MQ=255)
                         tGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaa                                                               >  2:757408/1‑93 (MQ=255)
                          ggTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaa                                                              >  1:889118/1‑93 (MQ=255)
                           gTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaa                                                             <  2:4156308/93‑1 (MQ=255)
                           gTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaa                                                             <  1:4850826/93‑1 (MQ=255)
                            tataCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaaa                                                            <  1:2263026/93‑1 (MQ=255)
                             ataCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAg                                                           <  2:5937168/93‑1 (MQ=255)
                              taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          <  1:6532958/93‑1 (MQ=255)
                              taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  1:4287464/1‑93 (MQ=255)
                              taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  2:3295008/1‑93 (MQ=255)
                               aCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGcc                                                         <  1:2809764/93‑1 (MQ=255)
                               aCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGcc                                                         <  1:2556552/93‑1 (MQ=255)
                                cTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCg                                                        <  2:1739907/93‑1 (MQ=255)
                                 ttACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGc                                                       >  2:373158/1‑93 (MQ=255)
                                 ttACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGc                                                       <  1:2936811/93‑1 (MQ=255)
                                    ccTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaa                                                    >  2:3624479/1‑93 (MQ=255)
                                    ccTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaa                                                    <  1:5790043/93‑1 (MQ=255)
                                     cTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaa                                                   <  1:1643266/93‑1 (MQ=255)
                                      tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  >  2:1637524/1‑93 (MQ=255)
                                      tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  >  2:6267543/1‑93 (MQ=255)
                                      tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  <  2:7137652/93‑1 (MQ=255)
                                      tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  <  1:6939745/93‑1 (MQ=255)
                                         acacATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               <  1:7898040/93‑1 (MQ=255)
                                           acaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAgg                                             <  2:4841475/93‑2 (MQ=255)
                                           acaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGt                                             >  2:111638/1‑93 (MQ=255)
                                            caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            <  1:7126115/93‑1 (MQ=255)
                                            caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            >  2:3856938/1‑93 (MQ=255)
                                             aTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCg                                           >  2:438944/1‑93 (MQ=255)
                                              ttGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGc                                          <  2:7862906/93‑1 (MQ=255)
                                              ttGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGc                                          >  1:4934276/1‑93 (MQ=255)
                                               tGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCa                                         >  2:1444080/1‑93 (MQ=255)
                                                 gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       <  1:5381238/93‑1 (MQ=255)
                                                  ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      >  2:5023521/1‑93 (MQ=255)
                                                  gccTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      <  2:5440158/91‑1 (MQ=255)
                                                    cTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacga                                    <  2:442959/93‑1 (MQ=255)
                                                     tGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgac                                   <  2:52555/93‑1 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTTCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGTCGCAAAAAATAGTCGCAAatgacg                                  >  1:1832733/1‑93 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgacg                                  >  1:5824165/1‑93 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgacg                                  <  1:130769/93‑1 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgaca                                  >  1:4813494/1‑92 (MQ=255)
                                                         tCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaa                               >  1:5573264/1‑93 (MQ=255)
                                                          cTGGATTCGACGGGATTTGCGAAACCCAGGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              <  2:7050034/93‑1 (MQ=255)
                                                          cTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              >  2:4274062/1‑93 (MQ=255)
                                                          cTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGACCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              <  1:3595757/93‑1 (MQ=255)
                                                            ggATTCGACGGGATTTGCGAAACCCAAGGCCCATGCCGAGGGGCGGGTGGCCTCGTAAAAAGCCGCCAAAAATAGGAGCAAACGACGAAAACt                            >  2:1819160/1‑93 (MQ=255)
                                                             gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  2:7263071/93‑1 (MQ=255)
                                                             gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:2126651/93‑1 (MQ=255)
                                                             gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:5023990/93‑1 (MQ=255)
                                                              aTTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTAc                          <  2:4715529/93‑1 (MQ=255)
                                                                tCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGc                        <  2:2577345/93‑1 (MQ=255)
                                                                  gACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCtt                      >  2:6802172/1‑93 (MQ=255)
                                                                   aCGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCttt                     >  2:4478728/1‑93 (MQ=255)
                                                                     gggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTag                   <  1:3205402/93‑1 (MQ=255)
                                                                     gggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACAACGAAAACTACGCTTTag                   >  1:2382691/1‑93 (MQ=255)
                                                                      ggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagc                  >  1:3174900/1‑93 (MQ=255)
                                                                      ggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagc                  <  1:5713026/93‑1 (MQ=255)
                                                                        aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                <  2:6442967/93‑1 (MQ=255)
                                                                        aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                >  2:3959201/1‑93 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               <  1:4805492/93‑1 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               <  2:1042681/93‑1 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               >  2:1198980/1‑93 (MQ=255)
                                                                             cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           <  1:5980881/93‑1 (MQ=255)
                                                                             cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           <  2:5659313/93‑1 (MQ=255)
                                                                             cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           <  1:4607165/93‑1 (MQ=255)
                                                                             cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           <  1:6521705/93‑1 (MQ=255)
                                                                             cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           <  2:1955384/93‑1 (MQ=255)
                                                                              gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  1:102727/93‑1 (MQ=255)
                                                                              gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  2:3135986/93‑1 (MQ=255)
                                                                               aaaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaata         >  1:6613430/1‑93 (MQ=255)
                                                                                 aCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAc       >  1:3532139/1‑93 (MQ=255)
                                                                                 aCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAc       <  2:5113102/93‑1 (MQ=255)
                                                                                    cAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTg    >  2:4222249/1‑93 (MQ=255)
                                                                                     aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGc   <  1:5134559/93‑1 (MQ=255)
                                                                                     aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGc   >  2:7529300/1‑93 (MQ=255)
                                                                                     aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACAAAAACTACGCTTTAGCAGCTTAATAACCTGc   <  1:5281787/93‑1 (MQ=255)
                                                                                         |                                                                                         
AGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753566‑2753744

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