BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	953,802	A→G	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.3%	503.7 / NA	136	V222A (GTA→GCA)	focA	formate channel
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (70/65); total (70/66)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTG  >  NC_000913/953705‑953897

  tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                             <  2:917603/100‑1 (MQ=255)

  tATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAAt                                                                                             <  2:629295/100‑1 (MQ=255)

   aTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATa                                                                                            <  1:1128946/100‑1 (MQ=255)

   aTCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATa                                                                                            >  2:2912705/1‑100 (MQ=255)

     taGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          <  1:837657/99‑1 (MQ=255)

     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          <  1:116276/100‑1 (MQ=255)

     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          >  1:1562468/1‑100 (MQ=255)

     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcc                                                                                          <  2:672843/100‑1 (MQ=255)

     cAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAc                                                                                           >  1:1570885/1‑99 (MQ=255)

      aGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATAcca                                                                                         >  2:837654/1‑99 (MQ=255)

         gatgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATc                                                                                      >  2:74093/1‑100 (MQ=255)

          atgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCg                                                                                     >  1:2705167/1‑100 (MQ=255)

           tgaAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCgg                                                                                    <  1:74094/100‑1 (MQ=255)

              aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                                 >  1:1911839/1‑100 (MQ=255)

              aaTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGAt                                                                                 >  2:2591629/1‑100 (MQ=255)

               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                                <  2:2551398/100‑1 (MQ=255)

               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                                <  2:1280017/100‑1 (MQ=255)

               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                                >  1:1426176/1‑100 (MQ=255)

               aTTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCAGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATc                                                                                <  2:1835630/100‑1 (MQ=255)

                ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                               <  1:126896/100‑1 (MQ=255)

                ttCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCa                                                                               <  2:2331172/100‑1 (MQ=255)

                  catcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                             <  1:686620/100‑1 (MQ=255)

                  catcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATa                                                                             >  1:1050340/1‑100 (MQ=255)

                   atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                            <  2:318530/100‑1 (MQ=255)

                   atcaCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATaa                                                                            >  1:377409/1‑100 (MQ=255)

                     caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                          >  2:516971/1‑100 (MQ=255)

                     caCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAAc                                                                          <  2:67610/100‑1 (MQ=255)

                      aCGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACa                                                                         <  2:2313252/100‑1 (MQ=255)

                       cGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACAt                                                                        >  1:987505/1‑100 (MQ=255)

                          tCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCagatcggaaga                                                                     >  1:875802/1‑90 (MQ=255)

                            aGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTg                                                                   <  1:2555693/100‑1 (MQ=255)

                             ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                  >  2:470100/1‑100 (MQ=255)

                             ggTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGc                                                                  <  2:1050335/100‑1 (MQ=255)

                                gagaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAt                                                               >  2:1836172/1‑100 (MQ=255)

                                 agaAAAATTTTCCGGTGCGGAACCGACTGCGGTCCAAAATTCTGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGAtn                                                              <  2:504554/100‑2 (MQ=255)

                                 agaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATa                                                              >  2:2758894/1‑100 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             >  2:2641448/1‑100 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             <  1:929159/100‑1 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             >  1:2380034/1‑100 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             >  1:726120/1‑100 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             >  1:1751260/1‑100 (MQ=255)

                                  gaAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATAc                                                             >  1:1781872/1‑100 (MQ=255)

                                    aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           <  2:2597382/100‑1 (MQ=255)

                                    aaaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtg                                                           <  1:2313819/100‑1 (MQ=255)

                                     aaaTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGCCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACtgt                                                          >  1:2112598/1‑100 (MQ=255)

                                       aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                        <  1:532632/100‑1 (MQ=255)

                                       aTTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGc                                                        <  2:1490283/100‑1 (MQ=255)

                                       aTTTTCCGGTGCAGAACCGACTGCGGTCC‑AAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCt                                                       >  1:1542905/1‑100 (MQ=255)

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

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

                                              ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                                 <  1:265381/100‑1 (MQ=255)

                                              ggTGCAGAACCGACTGCGGTCCAAAATTCCGGGAATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAAc                                                 >  2:1942766/1‑100 (MQ=255)

                                               gTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCtgccggtc                                                <  2:2774167/100‑8 (MQ=255)

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

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

                                                tGCAGAACCGACTGCGGTCCAAAATGCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCg                                               <  2:475846/100‑1 (MQ=255)

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

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

                                                   aGAACCGACTGCGGTCCAAAAGTCCGGGGATGCGAAGTCGCGTATTGCAATACCCATCGTGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTg                                            <  2:554901/100‑1 (MQ=255)

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

                                                    gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           >  2:2435956/1‑100 (MQ=255)

                                                    gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           >  2:2301756/1‑100 (MQ=255)

                                                    gAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           >  2:1016013/1‑100 (MQ=255)

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

                                                    gAACCGACTACGGTCCAAAATTTCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTgg                                           <  1:550297/100‑1 (MQ=255)

                                                        cGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac                                       <  2:353331/100‑1 (MQ=255)

                                                         gACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGCTTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaac                                       >  2:1461686/1‑99 (MQ=255)

                                                          aCTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGcaacaa                                     >  1:188759/1‑100 (MQ=255)

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

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

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

                                                                gTCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcg                               >  2:2055734/1‑100 (MQ=255)

                                                                 tCCAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgc                              >  1:2533433/1‑100 (MQ=255)

                                                                  gaaaagTTCCGGGGATGGGAAGTAGAGGATTGCAATAACAATAGGGATCATAAACATGTTTGTGATAATGTGTTAAAAAACGATGGCAACAAAAATagcg                             >  2:2839393/3‑100 (MQ=255)

                                                                  ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                             >  2:2697663/1‑100 (MQ=255)

                                                                  ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                             >  2:180725/1‑100 (MQ=255)

                                                                  ccAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATcgcg                             >  2:2697529/1‑100 (MQ=255)

                                                                   cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                            >  1:651833/1‑100 (MQ=255)

                                                                   cAAAATTCCGGGGATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGa                            >  2:2856894/1‑100 (MQ=255)

                                                                       gatctggcgATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATATTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCgg                        <  2:1603673/92‑1 (MQ=255)

                                                                         tCCGGGGATGCGAAGTCGCGGATTGCAATACCTATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGgca                      <  1:2709722/100‑1 (MQ=255)

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

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

                                                                              ggATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCCACCGGCAGCAcc                 <  1:1957864/100‑1 (MQ=255)

                                                                               gATGCGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCa                >  1:878494/1‑100 (MQ=255)

                                                                                   cGAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAAt            <  1:909003/100‑1 (MQ=255)

                                                                                    gAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg           <  2:1544306/100‑1 (MQ=255)

                                                                                    gAAGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg           <  2:1542896/100‑1 (MQ=255)

                                                                                     aaGTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATg           <  2:2454340/99‑1 (MQ=255)

                                                                                       gTCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAAc        <  1:1266203/100‑1 (MQ=255)

                                                                                        tCGCGGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAGTGAACg       <  2:1327180/100‑1 (MQ=255)

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

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

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

                                                                                         cgcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGc      <  1:2697545/100‑1 (MQ=255)

                                                                                          gcgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGGTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCt     >  1:293085/1‑100 (MQ=255)

                                                                                           cgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCtt    <  2:39309/100‑1 (MQ=255)

                                                                                           cgGATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCtt    <  1:2711787/100‑1 (MQ=255)

                                                                                            ggATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCttt   >  2:1241164/1‑100 (MQ=255)

                                                                                             gATTGCAATACCCATCGGGATCATAAACATGTTTGCGATACTGCGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTg  <  2:1911828/100‑1 (MQ=255)

GTTATCAGTGATGAAATTCATCACGGTCAGGTGAGAAAAATTTTCCGGTGCAGAACCGACTGCGGTCCAAAATTCCGGGGATGCGAAGTCGCGGATTACAATACCCATCGGGATCATAAACATGTTTGCGATACTGTGCTCAAAACCGCTGGCAACAAACATCGCGACCGGCAGCACCATAATGAACGCTTTG  >  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: `‑`