Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,357,653 A→G 7.6% L250L (CTA→CTG yfaD → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,6530AG 7.6% 270.2 / 7.9 91L250L (CTA→CTGyfaDtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base G (2/5);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAA  >  NC_000913/2357510‑2357793
                                                                                                                                               |                                                                                                                                            
ttCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACgg                                                                                                                                           >  5:193233/1‑147 (MQ=255)
ttCACTAATGACAGCCAGCTACAAACACTGTTAAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACgg                                                                                                                                           >  6:370380/1‑147 (MQ=255)
     tAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACAATTGCTGAACGGCTACTTCAGGAAg                                                                                                                                    >  4:152804/1‑149 (MQ=255)
       aTGACAGCCAGCTACTAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAg                                                                                                                                    <  7:336648/147‑1 (MQ=255)
       aTGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                  >  8:344734/1‑149 (MQ=255)
        tGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAggg                                                                                                                                  >  1:310052/1‑148 (MQ=255)
          aCAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCAt                                                                                                                               <  3:102711/149‑1 (MQ=255)
                   tACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggc                                                                                                                      >  3:265878/1‑149 (MQ=255)
                     cAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggc                                                                                                                      <  6:193233/147‑1 (MQ=255)
                      aaaCACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATtggctgg                                                                                                                   <  4:22200/149‑1 (MQ=255)
                              tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                             <  3:231070/147‑1 (MQ=255)
                              tttAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAgg                                                                                                           >  4:274240/1‑149 (MQ=255)
                                tAATTATCTGCTGCAATGCGGCGATACCTCGCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACGACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGTCATCAAATTGGCTGGCAGGAAGGt                                                                                                          <  4:176954/148‑1 (MQ=255)
                                 aaTTATCTGCTGCAATGCAGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAt                                                                                                        >  8:117529/1‑149 (MQ=255)
                                     aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAg                                                                                                            <  6:162499/141‑1 (MQ=255)
                                     aTCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAg                                                                                                            >  5:162499/1‑141 (MQ=255)
                                      tCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCAt                                                                                                    <  7:344734/148‑1 (MQ=255)
                                        tgctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                 >  7:293602/1‑149 (MQ=255)
                                         gctgcAATGCGGCGATACCTCCCGTTTCTCCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                 >  1:170049/1‑148 (MQ=255)
                                         gctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGaa                                                                                                 >  1:301929/1‑148 (MQ=255)
                                          ctgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                >  8:116852/1‑148 (MQ=255)
                                           tgcAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACaa                                                                                              >  6:285211/1‑149 (MQ=255)
                                             cAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAg                                                                                             <  4:265878/148‑1 (MQ=255)
                                             cAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATCACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGc                                                                                            >  6:303061/1‑149 (MQ=255)
                                                         aCCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAACGTATGCATGAACAAGCCATTAAAATTGc                                                                                >  6:58091/1‑149 (MQ=255)
                                                          ccTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCt                                                                               <  2:310052/149‑1 (MQ=255)
                                                           cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                    >  6:212937/1‑143 (MQ=255)
                                                           cTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTaaaa                                                                                    <  5:212937/143‑1 (MQ=255)
                                                              ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTTGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                             <  5:370380/147‑1 (MQ=255)
                                                              ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                             <  7:117529/147‑1 (MQ=255)
                                                              ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCttt                                                                             <  7:293313/147‑1 (MQ=255)
                                                              ccGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGATGAACGGCTACGTAAGGAAGAGCAACAAAATGAAAGGCAGTAAGGAATGCAATAACACTCAAATAAAATTGCTTTg                                                                            >  5:135340/1‑148 (MQ=255)
                                                                gTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgcgc                                                                         <  1:325374/149‑1 (MQ=255)
                                                                 tttCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGc                                                                                                      <  2:13101/119‑1 (MQ=255)
                                                                 tttCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGc                                                                                                      >  1:13101/1‑119 (MQ=255)
                                                                    cACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGGTTGAACATGCCATTAAACTTGCTTTGAGCATGc                                                                     >  7:279348/1‑149 (MQ=25)
                                                                     aCCCGTTTTATTTATGAGATTGCCGAACGTTCACCACTATAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaaa                                                                                                                           <  2:101371/94‑1 (MQ=18)
                                                                     aCCCGTTTTATTGATGAGTTTGCCTAACGTTCATCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACAt                                                                                                                               <  6:135951/90‑1 (MQ=18)
                                                                     aCCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATtggctgg                                                                                                                   <  8:108682/102‑1 (MQ=255)
                                                                     aCCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGTCAGGAAGGGCATCAAATTGGCTTGCAGGAATGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCt                                                                    <  5:303061/149‑1 (MQ=255)
                                                                     aCCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGACATCAAATtggctgg                                                                                                                   >  7:108682/1‑102 (MQ=255)
                                                                      cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCAt                                                                       >  1:99689/1‑145 (MQ=255)
                                                                      cccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCAt                                                                       <  2:99689/145‑1 (MQ=255)
                                                                       ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGTGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATCCTg                                                                   >  6:223207/1‑148 (MQ=255)
                                                                       ccGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                                   >  8:337236/1‑148 (MQ=255)
                                                                          ttttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGcc                                                                                           <  8:378685/121‑1 (MQ=255)
                                                                          ttttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGcc                                                                                           >  7:378685/1‑121 (MQ=255)
                                                                          ttttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                           >  6:383549/1‑137 (MQ=255)
                                                                          ttttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTgc                                                                           <  5:383549/137‑1 (MQ=255)
                                                                          ttttATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGAATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAg                                                                >  1:109684/1‑148 (MQ=255)
                                                                              aTTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAgg                                                            <  6:279308/148‑1 (MQ=255)
                                                                               ttGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGc                                                          >  8:310516/1‑149 (MQ=255)
                                                                                  aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCttt                                                       <  5:285211/149‑1 (MQ=255)
                                                                                  aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCttt                                                       >  6:298204/1‑149 (MQ=255)
                                                                                  aggagATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCtt                                                        <  2:109684/148‑1 (MQ=255)
                                                                                          gCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTg                                                <  5:58091/148‑1 (MQ=255)
                                                                                           ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCa                                                                                                                             <  4:46500/70‑1 (MQ=255)
                                                                                           ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCa                                                                                                                             >  3:46500/1‑70 (MQ=255)
                                                                                           ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGt                                                 >  5:256560/1‑146 (MQ=255)
                                                                                           ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTga                                               <  2:32113/148‑1 (MQ=255)
                                                                                           ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTg                                                <  8:206901/147‑1 (MQ=255)
                                                                                           ccGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTg                                                >  7:206901/1‑147 (MQ=255)
                                                                                            cGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTgaga                                             <  3:274240/149‑1 (MQ=255)
                                                                                             gAACGTTCACCACTACAACAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATTAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAAc                                                   <  1:98238/142‑1 (MQ=255)
                                                                                             gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaaa                                                                                                                           >  3:203492/1‑70 (MQ=255)
                                                                                             gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCaaa                                                                                                                           <  4:203492/70‑1 (MQ=255)
                                                                                             gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGTCAGGACGGGCATAAAATTGGCTGGCATGAAGGTATGCAGGCACAAGCCAGTAACATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGagagcg                                            >  2:98238/1‑147 (MQ=25)
                                                                                             gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                <  7:387952/97‑1 (MQ=255)
                                                                                             gAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAAc                                                                                                >  8:387952/1‑97 (MQ=255)
                                                                                              aaCGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTgaccag                                           <  2:170049/149‑5 (MQ=255)
                                                                                                  ttCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGc                                       <  7:310516/149‑1 (MQ=255)
                                                                                                     accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                    <  7:337236/149‑1 (MQ=255)
                                                                                                     accacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTgg                                    >  7:58114/1‑149 (MQ=255)
                                                                                                      ccacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGTAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGTCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                   >  2:334342/1‑149 (MQ=255)
                                                                                                       cacTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGc                                   >  7:236470/1‑148 (MQ=255)
                                                                                                            cAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGTAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGTCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACaa                              <  1:334342/148‑1 (MQ=255)
                                                                                                             aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAcc                            >  4:323658/1‑149 (MQ=255)
                                                                                                             aaaaGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAAcc                            <  6:256560/149‑1 (MQ=255)
                                                                                                                     agaTTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCAc                    <  8:58114/149‑1 (MQ=255)
                                                                                                                        ttAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGGTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGAGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTgc                 >  5:185169/1‑148 (MQ=255)
                                                                                                                        ttAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGa                 >  4:312821/1‑149 (MQ=255)
                                                                                                                           aTGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGc              <  3:323658/149‑1 (MQ=255)
                                                                                                                            tGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCt             <  7:116852/149‑1 (MQ=255)
                                                                                                                              aCTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGGATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACAGATGCTGa           <  1:77911/149‑1 (MQ=255)
                                                                                                                                tATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGata         >  3:293115/1‑149 (MQ=255)
                                                                                                                                     cTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTcc     <  8:293602/148‑1 (MQ=255)
                                                                                                                                      tGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCg    >  4:238487/1‑148 (MQ=255)
                                                                                                                                       gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCa                     <  6:226290/130‑1 (MQ=255)
                                                                                                                                       gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCa                     >  5:226290/1‑130 (MQ=255)
                                                                                                                                       gAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa  >  7:230691/1‑149 (MQ=255)
                                                                                                                                        aaCGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGaa  <  8:236470/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                            
TTCACTAATGACAGCCAGCTACAAACACTGTTTAATTATCTGCTGCAATGCGGCGATACCTCCCGTTTCACCCGTTTTATTGAGGAGATTGCCGAACGTTCACCACTACAAAAGGAGAGATTAATGACTATTGCTGAACGGCTACGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTATGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAA  >  NC_000913/2357510‑2357793

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: