Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 167,651 (G)5→6 coding (168/2244 nt) fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_000091167,6461.G98.6% 540.5 / 1.2 140coding (163/2244 nt)fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (80/58);  total (81/59)
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.35e-01

Only 100 of 150 total aligned reads displayed.

GCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑T‑GGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTA  >  AC_000091/167557‑167730
                                                                                           |                                                                                     
gCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑t                                                                                        <  2:2831878/90‑1 (MQ=255)
gCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑t                                                                                        <  1:545557/90‑1 (MQ=255)
 cGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑T‑g                                                                                      >  1:1072401/1‑90 (MQ=255)
 cGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑T‑g                                                                                      <  2:2144578/90‑1 (MQ=255)
    cATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑T‑gggg                                                                                   >  1:2083214/1‑90 (MQ=255)
    cATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGAGCCGCAAGAAAGCGCA‑T‑gggg                                                                                   >  1:1251274/1‑90 (MQ=255)
            tttATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑Cctgct                                                                           >  2:2955505/1‑90 (MQ=255)
             ttATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑Cctgctg                                                                          >  1:3306498/1‑90 (MQ=255)
             ttATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑Cctgctg                                                                          <  1:504480/90‑1 (MQ=255)
               aTGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCa                                                                        >  2:1497815/1‑90 (MQ=255)
                  cacaGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACt                                                                     <  2:1793046/90‑1 (MQ=255)
                    caGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTAt                                                                   >  1:1606863/1‑90 (MQ=255)
                    caGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTAt                                                                   <  2:1610806/90‑1 (MQ=255)
                     aGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTAtt                                                                  >  1:1970185/1‑90 (MQ=255)
                      ggCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTg                                                                 <  1:76344/90‑1 (MQ=255)
                      ggCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTg                                                                 >  2:1601369/1‑90 (MQ=255)
                        cAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcg                                                               >  1:851906/1‑90 (MQ=255)
                        cAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcg                                                               <  2:900544/90‑1 (MQ=255)
                         aGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcgg                                                              >  1:1577230/1‑90 (MQ=255)
                         aGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcgg                                                              <  2:827887/90‑1 (MQ=255)
                         aGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcgg                                                              >  1:1837922/1‑90 (MQ=255)
                         aGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcgg                                                              >  1:3385424/1‑90 (MQ=255)
                          gCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcggc                                                             >  2:1701107/1‑90 (MQ=255)
                           cGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcggcg                                                            >  1:1376713/1‑90 (MQ=255)
                           cGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTgcggcg                                                            >  1:2874921/1‑90 (MQ=255)
                             gTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGcgcg                                                          >  1:662587/1‑90 (MQ=255)
                              ttGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGa                                                         <  2:319655/90‑1 (MQ=255)
                              ttGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGa                                                         >  1:2007810/1‑90 (MQ=255)
                               tGAACCGAAAGAAGACACTATCCCCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGAc                                                        <  1:3301513/90‑1 (MQ=255)
                               tGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGAc                                                        >  2:2526974/1‑90 (MQ=255)
                               tGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCAGCAACTATTGCGGCGCGAc                                                        >  1:1110321/1‑90 (MQ=255)
                                gAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACa                                                       >  2:2865565/1‑90 (MQ=255)
                                gAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACa                                                       <  2:664061/90‑1 (MQ=255)
                                 aaCCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAg                                                      <  1:3261057/90‑1 (MQ=255)
                                 aaCCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAg                                                      >  1:2412915/1‑90 (MQ=255)
                                  aCCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGt                                                     >  1:2853688/1‑90 (MQ=255)
                                   ccGAAAGAAGACACTATCACCGTTACCGCTGCACCTGTGTCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTc                                                    <  2:1363665/90‑1 (MQ=255)
                                   ccGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTc                                                    >  1:3354571/1‑90 (MQ=255)
                                      aaaGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGc                                                 >  1:2140640/1‑90 (MQ=255)
                                       aagaagACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCt                                                <  2:3130178/90‑1 (MQ=255)
                                        agaggaCACTATCTCCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTa                                               >  2:3359234/1‑90 (MQ=255)
                                        agaagaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTa                                               >  1:2505217/1‑90 (MQ=255)
                                         gaagaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTAc                                              >  2:1480638/1‑90 (MQ=255)
                                           agaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTCTTGCGGCGCGACAGTCTGCTACCg                                            <  1:986278/90‑1 (MQ=255)
                                           agaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCg                                            <  2:1890525/90‑1 (MQ=255)
                                            gaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCgg                                           >  1:2232222/1‑90 (MQ=255)
                                            gaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCgg                                           >  1:753726/1‑90 (MQ=255)
                                                 tATCCCCGTTACCGCTGCCCCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTa                                      <  1:513383/90‑1 (MQ=255)
                                                 tATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGACTGCTACCGGCACTa                                      >  2:1247865/1‑90 (MQ=255)
                                                  aTCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTaa                                     >  1:2553195/1‑90 (MQ=255)
                                                   tCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTaaa                                    <  1:91644/90‑1 (MQ=255)
                                                   tCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTaaa                                    >  2:2635499/1‑90 (MQ=255)
                                                   tCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTaaa                                    >  1:279267/1‑90 (MQ=255)
                                                    cACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTaaaa                                   >  2:1325718/1‑90 (MQ=255)
                                                      ccGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAAcc                                 <  1:3113601/90‑1 (MQ=255)
                                                      ccGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAAcc                                 <  2:2159317/90‑1 (MQ=255)
                                                      ccGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAAcc                                 >  1:3271589/1‑90 (MQ=255)
                                                      ccGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAAcc                                 <  1:235298/90‑1 (MQ=255)
                                                       cGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCg                                <  2:794652/90‑1 (MQ=255)
                                                       cGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCg                                >  1:2483388/1‑90 (MQ=255)
                                                        gTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGa                               >  1:3252025/1‑90 (MQ=255)
                                                         ttACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGAt                              >  1:1643431/1‑90 (MQ=255)
                                                            ccGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACg                           <  1:2789992/90‑1 (MQ=255)
                                                            ccGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACg                           <  2:1410782/90‑1 (MQ=255)
                                                             cGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGc                          >  2:2425671/1‑90 (MQ=255)
                                                             cGCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGc                          >  1:485355/1‑90 (MQ=255)
                                                              gCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGcc                         <  2:3380969/90‑1 (MQ=255)
                                                              gCTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGcc                         <  2:1989970/90‑1 (MQ=255)
                                                               cTGCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCg                        >  1:3073927/1‑90 (MQ=255)
                                                                 gCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGAt                      >  2:2052005/1‑90 (MQ=255)
                                                                 gCACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATAAGCCGAt                      <  1:1105232/90‑1 (MQ=255)
                                                                  cACCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGAtt                     <  1:1179733/90‑1 (MQ=255)
                                                                   aCCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTc                    >  2:2972083/1‑90 (MQ=255)
                                                                   aCCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTc                    >  2:670930/1‑90 (MQ=255)
                                                                   aCCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTc                    <  1:1470193/90‑1 (MQ=255)
                                                                   aCCTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTc                    >  1:434650/1‑90 (MQ=255)
                                                                    ccTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCa                   >  2:809862/1‑90 (MQ=255)
                                                                     cTGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaa                  <  2:2739683/90‑1 (MQ=255)
                                                                      tGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaaa                 >  1:2538529/1‑90 (MQ=255)
                                                                      tGCGCCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaaa                 >  2:127644/1‑90 (MQ=255)
                                                                       gcgcCGCAAGAGAGCGCAGG‑GGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaaaa                <  2:1288563/90‑1 (MQ=255)
                                                                       gcgcCGCAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaaaa                <  1:884370/90‑1 (MQ=255)
                                                                        ggccgcAAGAAAGCGCA‑T‑GGGGGCCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaaaaa               <  1:844788/89‑1 (MQ=255)
                                                                           cgcAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTg            <  1:3134352/90‑1 (MQ=255)
                                                                           cgcAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTg            <  2:391777/90‑1 (MQ=255)
                                                                           cgcAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTg            >  1:2663510/1‑90 (MQ=255)
                                                                             cAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGcc          >  1:2220871/1‑90 (MQ=255)
                                                                             cAAGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGcc          <  1:2677067/90‑1 (MQ=255)
                                                                              aaGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCca         <  2:2773996/90‑1 (MQ=255)
                                                                              aaGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCca         <  2:1015287/90‑1 (MQ=255)
                                                                               aGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCcac        <  2:2637015/90‑1 (MQ=255)
                                                                               aGAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCcac        >  2:1994156/1‑90 (MQ=255)
                                                                                gAAAGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCcaca       <  2:1160269/90‑1 (MQ=255)
                                                                                 aaaGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAg      <  2:1111566/90‑1 (MQ=255)
                                                                                  aaGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGt     >  1:2835074/1‑90 (MQ=255)
                                                                                   aGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTc    <  2:201671/90‑1 (MQ=255)
                                                                                   aGCGCA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTc    <  1:2907215/90‑1 (MQ=255)
                                                                                    gcgcA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTTCCACAGTCt   <  2:355620/90‑1 (MQ=255)
                                                                                    gcgcA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCt   <  2:1140807/90‑1 (MQ=255)
                                                                                     cgcA‑TGGGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTa  <  1:103386/90‑1 (MQ=255)
                                                                                           |                                                                                     
GCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA‑T‑GGGGG‑CCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTA  >  AC_000091/167557‑167730

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