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.8% 1294.3 / NA 405intergenic (+256/‑525)smpB/intAsmall protein B/prophage CP4‑57 integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (150/254);  minor base A (1/0);  total (151/254)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 418 total aligned reads displayed.

GAGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753565‑2753744
                                                                                          |                                                                                         
 aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  2:4967794/93‑1 (MQ=255)
 aGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  1:5605005/93‑1 (MQ=255)
 aGGCTGGTCATGGCACTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCa                                                                                        <  1:5524118/93‑1 (MQ=255)
  ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       <  1:7984516/93‑1 (MQ=255)
  ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCAt                                                                                       <  2:5583392/93‑1 (MQ=255)
  ggCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCAGAACCCAAGGCGCAt                                                                                       >  1:3711206/1‑93 (MQ=255)
    gtGGGCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGc                                                                                     <  2:8416133/92‑1 (MQ=255)
     tGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGcc                                                                                    <  1:9465959/93‑1 (MQ=255)
       gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  <  2:3323377/93‑1 (MQ=255)
       gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  <  2:7154474/93‑1 (MQ=255)
       gTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGa                                                                                  <  2:1035758/93‑1 (MQ=255)
        tCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAg                                                                                 <  2:292097/93‑1 (MQ=255)
         cATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgg                                                                                <  1:9286610/93‑1 (MQ=255)
          aTGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAggg                                                                               >  2:4239416/1‑93 (MQ=255)
          aTGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAggg                                                                               <  2:1599882/93‑1 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              <  1:2142571/93‑1 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              >  2:2376077/1‑93 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              >  2:6022740/1‑93 (MQ=255)
           tGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAgggg                                                                              >  1:1000840/1‑93 (MQ=255)
             gcgcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCg                                                                            <  1:8945693/93‑1 (MQ=255)
               gcTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGt                                                                          <  2:6635889/93‑1 (MQ=255)
                cTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGtt                                                                         >  1:1423188/1‑93 (MQ=255)
                 tCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTg                                                                        >  1:2065071/1‑93 (MQ=255)
                 tCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTg                                                                        >  2:4510006/1‑93 (MQ=255)
                  cATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGCTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTgg                                                                       >  1:5828402/1‑93 (MQ=255)
                  cATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTgg                                                                       >  2:6791204/1‑93 (MQ=255)
                  cATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTgg                                                                       >  1:1921987/1‑93 (MQ=255)
                   aTAAATCTGGTATACTTACCGTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGc                                                                      <  1:8496895/93‑1 (MQ=255)
                    tAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGcc                                                                     <  2:4752319/93‑1 (MQ=255)
                     aaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCt                                                                    <  2:4387418/93‑1 (MQ=255)
                     aaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCt                                                                    <  2:5095397/93‑1 (MQ=255)
                     aaaTCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCt                                                                    <  2:4207802/93‑1 (MQ=255)
                          tGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaa                                                               >  1:5354401/1‑93 (MQ=255)
                           ggTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaa                                                              >  1:4303354/1‑93 (MQ=255)
                           ggTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTATGCCTCGTaga                                                              >  2:3643629/1‑91 (MQ=255)
                            gTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTaaaa                                                             >  2:9201714/1‑93 (MQ=255)
                               taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  2:3824849/1‑93 (MQ=255)
                               taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  1:8825988/1‑93 (MQ=255)
                               taCTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGc                                                          >  1:4961374/1‑93 (MQ=255)
                                 cTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCg                                                        >  1:2400234/1‑93 (MQ=255)
                                 cTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCCTAAAAACCCg                                                        >  1:2553638/1‑93 (MQ=255)
                                   tACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGGCTCGTAAAAAGCCGCa                                                      >  2:8879495/1‑93 (MQ=255)
                                   tACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCa                                                      >  2:9805081/1‑93 (MQ=255)
                                    aCCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaa                                                     <  1:8058250/93‑1 (MQ=255)
                                    aCCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaa                                                     >  2:7892422/1‑93 (MQ=255)
                                     ccTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaa                                                    >  1:9078745/1‑93 (MQ=255)
                                      cTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaa                                                   >  1:9624372/1‑93 (MQ=255)
                                        ttACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCaaaaaa                                                 >  2:8813551/1‑93 (MQ=255)
                                         tACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAAt                                                <  2:6458005/93‑1 (MQ=255)
                                         tACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAAt                                                >  2:1261964/1‑93 (MQ=255)
                                          acacATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               <  2:7490579/93‑1 (MQ=255)
                                          acacATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATa                                               <  1:6980136/93‑1 (MQ=255)
                                             caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            >  2:9682120/1‑93 (MQ=255)
                                             caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            >  2:806002/1‑93 (MQ=255)
                                             caTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTc                                            <  2:8306973/93‑1 (MQ=255)
                                              aTTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCg                                           >  2:132720/1‑93 (MQ=255)
                                               atGGGGCCGATTCTGGATTCGACTGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGc                                          <  1:8577285/92‑1 (MQ=255)
                                                 ggggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaa                                        <  2:7265867/93‑1 (MQ=255)
                                                  gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       >  2:3448720/1‑93 (MQ=255)
                                                  gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       <  1:296444/93‑1 (MQ=255)
                                                  gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       <  2:401843/93‑1 (MQ=255)
                                                  gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       <  2:867144/93‑1 (MQ=255)
                                                  gggCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCaaa                                       >  2:2940942/1‑93 (MQ=255)
                                                    gCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacg                                     <  1:6330629/93‑1 (MQ=255)
                                                     gtGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacga                                    <  2:3068085/92‑1 (MQ=255)
                                                     gtGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacga                                    <  1:7273108/92‑1 (MQ=255)
                                                     cTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAacga                                    >  2:2335776/1‑93 (MQ=255)
                                                        aTTCTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAAcgacga                                 <  1:3630046/93‑1 (MQ=255)
                                                         ttCTGTACTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGTCGTTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaa                                <  1:344178/93‑1 (MQ=255)
                                                           cTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              <  2:7649869/93‑1 (MQ=255)
                                                           cTGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGaaaa                              <  2:9516468/93‑1 (MQ=255)
                                                            tGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAAc                             <  1:5466188/93‑1 (MQ=255)
                                                            tGGATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAAc                             <  1:5945608/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  2:1922260/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:2908304/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:3138548/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  1:1754697/93‑1 (MQ=255)
                                                              gATTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTa                           <  2:623761/93‑1 (MQ=255)
                                                               aTTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTAc                          <  2:8247902/93‑1 (MQ=255)
                                                               aTTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTAc                          <  1:8750401/93‑1 (MQ=255)
                                                               aTTCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTAc                          <  1:1169399/93‑1 (MQ=255)
                                                                 tCGACGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGc                        <  1:9166566/93‑1 (MQ=255)
                                                                    aCGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCttt                     <  1:4103211/93‑1 (MQ=255)
                                                                     cGGGATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTa                    >  1:7633683/1‑93 (MQ=255)
                                                                      gggATTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTag                   <  2:5397047/93‑1 (MQ=255)
                                                                        gATTTGCGAACCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagca                 <  2:6930473/93‑1 (MQ=255)
                                                                         aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                <  1:6843559/93‑1 (MQ=255)
                                                                         aTTTGCGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTagcag                <  1:7495953/93‑1 (MQ=255)
                                                                              cGAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaa           <  1:6235385/93‑1 (MQ=255)
                                                                               gAAACCCAAGGCGCATGCCGAGGGGCGGTTTGCCTTGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  1:7144646/93‑1 (MQ=255)
                                                                               gAAACCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaat          <  1:1226067/93‑1 (MQ=255)
                                                                                aaaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaata         >  2:4632063/1‑93 (MQ=255)
                                                                                 aaCCCAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTtaataa        >  1:5074520/1‑93 (MQ=255)
                                                                                   cccAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAAcc      <  1:4823368/93‑1 (MQ=255)
                                                                                     cAAGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTg    <  1:1663547/93‑1 (MQ=255)
                                                                                      aaGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGc   <  1:6437327/93‑1 (MQ=255)
                                                                                       aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  >  1:270246/1‑93 (MQ=255)
                                                                                       aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  <  2:2680934/93‑1 (MQ=255)
                                                                                       aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  >  2:1154469/1‑93 (MQ=255)
                                                                                       aGGCGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCt  <  1:418765/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAGGCTGGTCATGGCGCTCATAAATCTGGTATACTTACCTTTACACATTGGGGCTGATTCTGGATTCGACGGGATTTGCGAAACCCAAGGTGCATGCCGAGGGGCGGTTGGCCTCGTAAAAAGCCGCAAAAAATAGTCGCAAACGACGAAAACTACGCTTTAGCAGCTTAATAACCTGCT  >  NC_000913/2753565‑2753744

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