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,6550TC98.5% 640.1 / ‑1.5 203intergenic (+256/‑525)smpB/intAsmall protein B/prophage CP4‑57 integrase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (61/139);  total (62/141)
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.37e-01

Only 100 of 206 total aligned reads displayed.

AGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753566‑2753744
                                                                                         |                                                                                         
aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  1:4176214/93‑1 (MQ=255)
aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  1:3449978/93‑1 (MQ=255)
 ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       <  2:1706036/93‑1 (MQ=255)
 ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       <  1:287491/93‑1 (MQ=255)
  gCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATg                                                                                      >  1:3736588/1‑93 (MQ=255)
  gCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATg                                                                                      <  1:2232321/93‑1 (MQ=255)
   cTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGc                                                                                     >  2:3054100/1‑93 (MQ=255)
   cTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGc                                                                                     <  2:1407726/93‑1 (MQ=255)
    tGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGcc                                                                                    <  2:1809011/93‑1 (MQ=255)
     ggTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCg                                                                                   <  1:2274827/93‑1 (MQ=255)
       tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 >  1:2902811/1‑93 (MQ=255)
       tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 >  2:3766147/1‑93 (MQ=255)
        cATGGCGCTCATAAATCTGGTATACTTACCTTTACACGTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgg                                                                                <  1:4336368/93‑1 (MQ=255)
        cATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgg                                                                                <  1:2153688/93‑1 (MQ=255)
         aTGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAggg                                                                               >  2:3405564/1‑93 (MQ=255)
         aTGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAggg                                                                               <  1:1140310/93‑1 (MQ=255)
          tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              <  2:3106885/93‑1 (MQ=255)
          tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              <  2:4382022/93‑1 (MQ=255)
           ggCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGc                                                                             <  2:4155663/93‑1 (MQ=255)
             cgcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCgg                                                                           <  1:2535958/93‑1 (MQ=255)
              gcTCATATATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATCCCGAGGGGCGGt                                                                          >  2:2978702/1‑93 (MQ=255)
              gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          <  1:3323612/93‑1 (MQ=255)
              gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          <  2:4373635/93‑1 (MQ=255)
              gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          >  2:1299929/1‑93 (MQ=255)
               cTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGtt                                                                         <  1:2594333/93‑1 (MQ=255)
                     aaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTc                                                                   >  2:1519843/1‑93 (MQ=255)
                      aTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCg                                                                  >  2:2735577/1‑93 (MQ=255)
                        cTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTa                                                                <  1:4197486/93‑1 (MQ=255)
                         tGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaa                                                               <  1:4627904/93‑1 (MQ=255)
                         tGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaa                                                               >  2:2434466/1‑93 (MQ=255)
                         tGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaa                                                               <  2:4277322/93‑1 (MQ=255)
                            tataCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaaa                                                            <  1:1688523/93‑1 (MQ=255)
                            gataCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaaa                                                            <  2:3282671/92‑1 (MQ=255)
                              taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          <  1:4249327/93‑1 (MQ=255)
                              taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  1:2955907/1‑93 (MQ=255)
                               aCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGcc                                                         <  2:2464432/93‑1 (MQ=255)
                                cTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCg                                                        >  2:819511/1‑93 (MQ=255)
                                 ttACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGc                                                       >  1:2858104/1‑93 (MQ=255)
                                 ttACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGc                                                       <  1:854193/93‑1 (MQ=255)
                                      tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  >  2:1353690/1‑93 (MQ=255)
                                      tttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaa                                                  >  1:3574467/1‑93 (MQ=255)
                                        tACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAAt                                                >  1:1800731/1‑93 (MQ=255)
                                        tACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAAt                                                <  1:2174160/93‑1 (MQ=255)
                                          cacaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAg                                              <  2:2738881/93‑1 (MQ=255)
                                           acaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGt                                             >  2:3711252/1‑93 (MQ=255)
                                           acaTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGt                                             <  2:1733165/93‑1 (MQ=255)
                                             aTTGGGGCTGATTCTGGATTCGACGGGATTTGCGGAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCg                                           >  1:1834590/1‑93 (MQ=255)
                                             aTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCg                                           <  1:1900239/93‑1 (MQ=255)
                                               tGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCa                                         <  2:2506075/93‑1 (MQ=255)
                                                ggggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaa                                        >  2:4672570/1‑93 (MQ=255)
                                                 gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       <  2:1943002/93‑1 (MQ=255)
                                                  ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      <  1:3610348/93‑1 (MQ=255)
                                                  ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      <  1:3045982/93‑1 (MQ=255)
                                                  ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      <  2:3681214/93‑1 (MQ=255)
                                                  ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAac                                      <  2:2766779/93‑1 (MQ=255)
                                                  ggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGAGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTTGCAAac                                      <  1:3843818/93‑1 (MQ=255)
                                                     tGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgac                                   >  1:996485/1‑93 (MQ=255)
                                                     tGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgac                                   <  1:2695230/93‑1 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgacg                                  >  2:3941088/1‑93 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgacg                                  >  2:878905/1‑93 (MQ=255)
                                                      gATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacgacg                                  >  1:4411536/1‑93 (MQ=255)
                                                       aTTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAAcgacga                                 >  1:3531335/1‑93 (MQ=255)
                                                           tGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGCAAAc                             >  1:3120881/1‑93 (MQ=255)
                                                             gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           >  1:156754/1‑93 (MQ=255)
                                                             gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:4372914/93‑1 (MQ=255)
                                                               ttCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACg                         <  2:1027318/93‑1 (MQ=255)
                                                                tCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGc                        <  2:1220581/93‑1 (MQ=255)
                                                                tCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGc                        <  2:4044154/93‑1 (MQ=255)
                                                                  gACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCtt                      <  2:59379/93‑1 (MQ=255)
                                                                   aCGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCttt                     <  1:1419332/93‑1 (MQ=255)
                                                                   aCGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCttt                     <  1:3797069/93‑1 (MQ=255)
                                                                   aCGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCttt                     <  2:2826961/93‑1 (MQ=255)
                                                                    cGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTa                    <  1:1992221/93‑1 (MQ=255)
                                                                    cGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTa                    <  2:2231598/93‑1 (MQ=255)
                                                                     gggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTag                   <  2:2914238/93‑1 (MQ=255)
                                                                     gggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTag                   <  2:2061627/93‑1 (MQ=255)
                                                                     gggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTag                   <  2:2195381/93‑1 (MQ=255)
                                                                       gATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagca                 <  1:2433034/93‑1 (MQ=255)
                                                                       gATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagca                 <  2:4648351/93‑1 (MQ=255)
                                                                       gATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagca                 <  2:1658740/93‑1 (MQ=255)
                                                                        aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                <  2:619000/93‑1 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               >  1:3420829/1‑93 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               <  1:2284030/93‑1 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               >  2:2922702/1‑93 (MQ=255)
                                                                         tttGCGAAACCCAAGGCGAATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcagc               <  1:1217176/93‑1 (MQ=255)
                                                                          ttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCt              <  1:1047351/93‑1 (MQ=255)
                                                                          ttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCt              <  2:44740/93‑1 (MQ=255)
                                                                          ttGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCt              <  1:714979/93‑1 (MQ=255)
                                                                            gCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTta            <  1:1324380/93‑1 (MQ=255)
                                                                              gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  1:4618725/93‑1 (MQ=255)
                                                                               aaaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaata         <  1:3375031/93‑1 (MQ=255)
                                                                               aaaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaata         <  2:3448447/93‑1 (MQ=255)
                                                                                 aCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAc       <  1:39479/93‑1 (MQ=255)
                                                                                 aCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAc       <  2:496333/93‑1 (MQ=255)
                                                                                 aCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAc       <  1:509562/93‑1 (MQ=255)
                                                                                 aCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAc       <  2:2698917/93‑1 (MQ=255)
                                                                                   ccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCt     <  1:3272798/93‑1 (MQ=255)
                                                                                   ccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCt     <  1:4125064/93‑1 (MQ=255)
                                                                                   ccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCt     >  1:4717328/1‑93 (MQ=255)
                                                                                      aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  <  1:4471308/93‑1 (MQ=255)
                                                                                         |                                                                                         
AGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753566‑2753744

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