BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	953,802	A→G	100%	V222A (GTA→GCA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	953,802	0	A	G	99.1%	406.9 / ‑4.6	115	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (50/64); total (51/64)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 115 total aligned reads displayed.

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTG  >  NC_000913/953705‑953897

gTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCa                                                                                               <  1:2601711/100‑1 (MQ=255)

  tgaccgTGATGAATTTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                             >  1:1985515/6‑100 (MQ=255)

    tCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           <  2:18960/100‑1 (MQ=255)

     ccgTGATGAATTTCATCACGGTCAGGTGAGAAAAATTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           >  1:1025679/3‑100 (MQ=255)

     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          >  1:234430/1‑100 (MQ=255)

      aGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                         >  1:23956/1‑100 (MQ=255)

      aGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAccc                                                                                         >  2:767996/1‑100 (MQ=255)

          atgaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCg                                                                                     <  1:159059/100‑1 (MQ=255)

           tgaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                    <  1:1051477/100‑1 (MQ=255)

           tgaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                    >  2:1909129/1‑100 (MQ=255)

           tgaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                    >  2:1226724/1‑100 (MQ=255)

           tgaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                    >  1:2915966/1‑100 (MQ=255)

            gaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                   >  1:231943/1‑100 (MQ=255)

            gaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                   >  2:1559790/1‑100 (MQ=255)

            gaAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                   >  2:535571/1‑100 (MQ=255)

            gaAATTCATCACGGTCAGGTGAGAAAAA‑GTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCggg                                                                                   >  1:2439882/1‑100 (MQ=255)

             aaaTTCATCACGGTCAGGTGAGAAAAA‑TTCTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGTATGCGAAGTCGCGGATTGCAATACCCATCGGGa                                                                                  <  1:1839986/100‑1 (MQ=255)

              aaTTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                                 <  1:163546/100‑1 (MQ=255)

              aaTTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGAGTGCAATACCCATCGGGAt                                                                                 >  1:3050332/1‑100 (MQ=255)

                ttCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                               <  2:392670/100‑1 (MQ=255)

                ttCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                               <  1:87155/100‑1 (MQ=255)

                 tcatcaCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCAt                                                                              <  2:1315420/100‑1 (MQ=255)

                  catcaCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                             >  1:985366/1‑100 (MQ=255)

                   atcaCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGCACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                            >  1:2789011/1‑100 (MQ=255)

                   atcaCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                            <  1:767991/100‑1 (MQ=255)

                    tcaCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGATCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGCTCATaaa                                                                           <  1:1606626/100‑1 (MQ=255)

                     caCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                          <  1:2439405/100‑1 (MQ=255)

                     caCGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                          <  1:1734297/100‑1 (MQ=255)

                       cGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACAt                                                                        >  1:951979/1‑100 (MQ=255)

                         gTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGt                                                                      <  2:10095/100‑1 (MQ=255)

                           cAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCCTCGGGATCATAAACATGttt                                                                    >  2:140252/1‑100 (MQ=255)

                             ggTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                  <  2:541630/100‑1 (MQ=255)

                             ggTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                  >  2:2949721/1‑100 (MQ=255)

                             ggTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                  >  1:1279953/1‑100 (MQ=255)

                              gTGAGAAAAA‑TTTTCCGTTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                                 >  2:643652/1‑100 (MQ=255)

                              gTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                                 <  1:38453/100‑1 (MQ=255)

                              gTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCg                                                                 <  2:1794496/100‑1 (MQ=255)

                               tGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGttt                                                                    >  1:2459025/1‑96 (MQ=255)

                               tGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCAATCGGGATCATAAACATGTTTGCGa                                                                >  1:541938/1‑100 (MQ=255)

                                gagaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               <  1:377291/100‑1 (MQ=255)

                                gagaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               >  1:484062/1‑100 (MQ=255)

                                gagaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               >  2:879067/1‑100 (MQ=255)

                                gagaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               >  2:693405/1‑100 (MQ=255)

                                gagaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               <  2:621097/100‑1 (MQ=255)

                                gagaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               <  2:480192/100‑1 (MQ=255)

                                  gaAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             >  1:1103413/1‑100 (MQ=255)

                                   aaaaa‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACt                                                            >  1:322445/1‑100 (MQ=255)

                                    aaaa‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           >  1:1402340/1‑100 (MQ=255)

                                    aaaa‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           <  1:982035/100‑1 (MQ=255)

                                            tCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTc                                                      <  1:2051389/98‑1 (MQ=255)

                                              cGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCaaaa                                                  >  2:1681564/1‑100 (MQ=255)

                                               ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                                 >  2:169568/1‑100 (MQ=255)

                                               ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                                 >  2:1064688/1‑100 (MQ=255)

                                                gTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAcc                                                >  2:1810816/1‑100 (MQ=255)

                                                 tGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg                                               >  1:1157925/1‑100 (MQ=255)

                                                  gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                              <  1:535567/100‑1 (MQ=255)

                                                  gCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGc                                              <  1:1559782/100‑1 (MQ=255)

                                                   cAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                             >  2:2810949/1‑100 (MQ=255)

                                                   cAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCt                                             <  1:2384563/100‑1 (MQ=255)

                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                            <  2:1014937/100‑1 (MQ=255)

                                                    aGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                            >  1:1149468/1‑100 (MQ=255)

                                                     gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           <  2:2915984/100‑1 (MQ=255)

                                                     gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           <  2:356079/100‑1 (MQ=255)

                                                     gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTATGCTCAAAACCGCTgg                                           <  2:1234222/100‑1 (MQ=255)

                                                      aaCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  2:2883697/100‑1 (MQ=255)

                                                      aaCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  1:140250/100‑1 (MQ=255)

                                                      aaCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  2:2352193/100‑1 (MQ=255)

                                                      aaCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  2:1261349/100‑1 (MQ=255)

                                                       aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCTGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGc                                          <  1:918300/99‑1 (MQ=255)

                                                       aCCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGca                                         <  2:2338671/100‑1 (MQ=255)

                                                          gACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaaca                                      <  2:231945/100‑1 (MQ=255)

                                                            cTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa                                    >  1:504051/1‑100 (MQ=255)

                                                            cTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaa                                    >  2:200652/1‑100 (MQ=255)

                                                             tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                   <  1:972533/100‑1 (MQ=255)

                                                             tGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCaacaaac                                   <  2:1413378/100‑1 (MQ=255)

                                                                 gTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                               <  2:1680288/100‑1 (MQ=255)

                                                                  tCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgc                              <  1:1034288/100‑1 (MQ=255)

                                                                  tCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAGACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgc                              <  2:2146787/100‑1 (MQ=255)

                                                                   ccAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                             >  2:587338/1‑100 (MQ=255)

                                                                    cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGGTCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                            <  1:2594497/100‑1 (MQ=255)

                                                                     aaaaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGAc                           <  2:562673/100‑1 (MQ=255)

                                                                       aaTTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCg                         <  2:2293727/100‑1 (MQ=255)

                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                       >  1:19222/1‑100 (MQ=255)

                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                       <  2:2846243/100‑1 (MQ=255)

                                                                         ttCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgc                       <  2:1882170/100‑1 (MQ=255)

                                                                          tCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                      <  2:794577/100‑1 (MQ=255)

                                                                          tCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                      <  1:1226718/100‑1 (MQ=255)

                                                                           ccGGGGATCGTAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAACCCGCTGGCAACAAACATCGCGACCGgca                      <  1:2596311/99‑1 (MQ=255)

                                                                            cGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgcagc                    >  2:1896221/1‑100 (MQ=255)

                                                                               ggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCAcc                 <  2:91280/100‑1 (MQ=255)

                                                                                gATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCa                <  2:1873138/100‑1 (MQ=255)

                                                                                  tGCTAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCAt               <  1:1026906/99‑1 (MQ=255)

                                                                                  tGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATa              >  1:379115/1‑100 (MQ=255)

                                                                                  tGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATa              <  1:357512/100‑1 (MQ=255)

                                                                                         tCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACg       <  2:1622203/100‑1 (MQ=255)

                                                                                          cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc      >  1:2801865/1‑100 (MQ=255)

                                                                                           gcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCt     >  2:738365/1‑100 (MQ=255)

                                                                                             ggATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCttt   <  2:2326403/100‑1 (MQ=255)

                                                                                              gATTGCAATACCCATCGAGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTg  <  1:624389/100‑1 (MQ=255)

                                                                                              gATTGCAATACCCAGCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTg  <  2:2539470/100‑1 (MQ=255)

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAA‑TTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTG  >  NC_000913/953705‑953897

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