BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,440,313	A→G	T119A (ACC→GCC)	ytfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,440,313	0	A	G	99.3%	439.2 / NA	145	T119A (ACC→GCC)	ytfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (56/87); total (57/88)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

AAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCC  >  NC_000913/4440223‑4440402

aaagGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGcc                                                                                         <  1:1125598/93‑1 (MQ=255)

aaagGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGcc                                                                                         >  2:4233806/1‑93 (MQ=255)

aaagGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGcc                                                                                         <  1:1320932/93‑1 (MQ=255)

aaagGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGcc                                                                                         <  2:2383967/93‑1 (MQ=255)

 aagGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCg                                                                                        <  1:4344833/93‑1 (MQ=255)

  agGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGa                                                                                       <  1:654413/93‑1 (MQ=255)

   gggCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGAt                                                                                      >  2:1449681/1‑93 (MQ=255)

   gggCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGAt                                                                                      <  1:1333834/93‑1 (MQ=255)

      cggcAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATaaa                                                                                   >  2:3859289/1‑93 (MQ=255)

      cggcAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATaaa                                                                                   <  2:3909299/93‑1 (MQ=255)

       ggcAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAg                                                                                  >  2:3099058/1‑93 (MQ=255)

        gcAGGTATTGATCGCGAAAGTCACGCCAGGCGTGCCGGTGTTAATGGGCGGCACCGATGTGGTATTGCGCGTCGGCGCGCGGGCCGATAAAGa                                                                                 <  2:695520/93‑1 (MQ=255)

         cAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGAc                                                                                <  2:4356993/93‑1 (MQ=255)

         cAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGAc                                                                                <  1:3288972/93‑1 (MQ=255)

         cAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGAc                                                                                >  1:2675393/1‑93 (MQ=255)

          aGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACt                                                                               >  1:3130284/1‑93 (MQ=255)

           ggTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTa                                                                              <  1:1297664/93‑1 (MQ=255)

            gTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAt                                                                             >  1:2934542/1‑93 (MQ=255)

             tATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCCCCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAtt                                                                            <  2:944932/93‑1 (MQ=255)

             tATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAtt                                                                            <  1:790389/93‑1 (MQ=255)

             tATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAtt                                                                            <  1:528246/93‑1 (MQ=255)

              aTTGATCTCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAttt                                                                           >  1:1357432/1‑93 (MQ=255)

              aTTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAttt                                                                           >  2:3326478/1‑93 (MQ=255)

              aTTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTAttt                                                                           <  2:689640/93‑1 (MQ=255)

               ttGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTg                                                                          >  1:2800460/1‑93 (MQ=255)

               ttGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCAGATAAAGACTATTTg                                                                          >  1:3248998/1‑93 (MQ=255)

               ttGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTGAATTGGCGGCACGGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTg                                                                          <  1:1438107/93‑1 (MQ=255)

                  aTCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGaaa                                                                       <  1:3449041/93‑1 (MQ=255)

                  aTCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGaaa                                                                       <  2:2831908/93‑1 (MQ=255)

                   tCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAAt                                                                      >  2:4541210/1‑93 (MQ=255)

                   tCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAAt                                                                      <  1:2828478/93‑1 (MQ=255)

                    cGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCAGAGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAAtt                                                                     <  1:2397249/93‑1 (MQ=255)

                      ccAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGc                                                                   <  1:4629702/93‑1 (MQ=255)

                       cAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCt                                                                  >  1:4312640/1‑93 (MQ=255)

                        aaaGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTc                                                                 >  1:1717239/1‑93 (MQ=255)

                           gTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGAt                                                              >  1:3792391/1‑93 (MQ=255)

                           gTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGAt                                                              >  1:4226402/1‑93 (MQ=255)

                           gTCACGCCAGCCGTGCCGTTGTTAATGGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGAt                                                              <  2:622082/93‑1 (MQ=255)

                             cACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATAc                                                            <  2:468699/93‑1 (MQ=255)

                              aCGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACt                                                           >  2:1231663/1‑93 (MQ=255)

                               cGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTc                                                          <  1:1330977/93‑1 (MQ=255)

                                 ccAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGc                                                        >  1:2323403/1‑93 (MQ=255)

                                   aGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGccc                                                      <  1:3912216/93‑1 (MQ=255)

                                   aGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGccc                                                      <  1:260511/93‑1 (MQ=255)

                                   aGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGGATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGgcc                                                      >  1:1639374/1‑93 (MQ=255)

                                    ggCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCg                                                     >  2:2852364/1‑93 (MQ=255)

                                    ggCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCg                                                     >  2:356922/1‑93 (MQ=255)

                                    ggCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCg                                                     >  2:3945920/1‑93 (MQ=255)

                                      ggTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGGGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGc                                                   <  1:3264280/92‑1 (MQ=255)

                                      ggTGCCGCTGTTAATGGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGc                                                   <  2:2568198/92‑1 (MQ=255)

                                        tGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTa                                                 <  1:1252259/93‑1 (MQ=255)

                                          ccGGTGTTAATTGGCGGCACGGTTTTGGTATTGCGCGGCGGCGCCCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTAtt                                               <  2:198514/93‑1 (MQ=255)

                                          ccGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTAtt                                               <  1:2148095/93‑1 (MQ=255)

                                          ccGGTGTTAAGTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGTCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTAtt                                               <  1:2721005/93‑1 (MQ=255)

                                           cGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTg                                              >  2:367523/1‑93 (MQ=255)

                                            ggTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTgg                                             <  2:2373236/93‑1 (MQ=255)

                                            ggTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTgg                                             >  1:2514803/1‑93 (MQ=255)

                                             gtgtTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGc                                            >  1:647400/1‑93 (MQ=255)

                                             gtgtTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGc                                            >  2:2689631/1‑93 (MQ=255)

                                               gtTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCAc                                          >  1:566868/1‑93 (MQ=255)

                                               gtTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCAc                                          >  1:179630/1‑93 (MQ=255)

                                                 tatttGGCGACACCGTTGTGGTATTGCGCGGCGGCGCGCGGGCCTATAAATACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACgg                                        <  2:3880075/90‑1 (MQ=255)

                                                  aaTTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGt                                       <  2:3675683/93‑1 (MQ=255)

                                                  aaTTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGt                                       <  1:1191983/93‑1 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     <  2:4356832/93‑1 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     >  2:2955766/1‑93 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     >  2:4064448/1‑93 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     >  1:1118269/1‑93 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     <  1:2893864/93‑1 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     >  1:4294440/1‑93 (MQ=255)

                                                    ttGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTAc                                     <  1:3812942/93‑1 (MQ=255)

                                                     tGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACt                                    <  1:3846844/93‑1 (MQ=255)

                                                        cggcACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGaa                                 <  2:4633520/93‑1 (MQ=255)

                                                        cgccACGGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAATTTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGaa                                 <  2:4505932/90‑1 (MQ=255)

                                                         ggcACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAAc                                <  2:2027409/93‑1 (MQ=255)

                                                         ggcACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAAc                                <  1:3007152/93‑1 (MQ=255)

                                                         ggcACCGATGTAGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAAc                                >  1:2656057/1‑93 (MQ=255)

                                                          gcACCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAAcc                               >  1:3534627/1‑93 (MQ=255)

                                                            aCCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAg                             >  2:723967/1‑93 (MQ=255)

                                                            aCCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAg                             <  1:3647467/93‑1 (MQ=255)

                                                            aCCGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAg                             <  1:3781564/93‑1 (MQ=255)

                                                             ccGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAgg                            <  2:4451766/93‑1 (MQ=255)

                                                             ccGATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAgg                            <  1:1278052/93‑1 (MQ=255)

                                                               gATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGc                          <  2:4245151/93‑1 (MQ=255)

                                                               gATGTGGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGc                          <  2:1431081/93‑1 (MQ=255)

                                                                 tgtgGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGa                        >  1:82888/1‑93 (MQ=255)

                                                                 tgtgGTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGa                        >  2:1038889/1‑93 (MQ=255)

                                                                   gggTATTGCGCGGCGGCGCGTGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTTTTGGCACGGTCCTGAACCAGGGCGAtt                      <  2:4665685/92‑1 (MQ=255)

                                                                    ggTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTa                     <  2:1283125/93‑1 (MQ=255)

                                                                    ggTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTa                     <  2:4614107/93‑1 (MQ=255)

                                                                     gTATTGCGCGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTAt                    <  1:4010636/93‑1 (MQ=255)

                                                                          gcgcGGCGGCGCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGaaaa               <  1:1561759/93‑1 (MQ=255)

                                                                                 ggcgCGCGGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCaaa        <  1:1849124/93‑1 (MQ=255)

                                                                                  gcgcgcGGGCCGATAAAGACTATTTGAAATTGGTCGATACTCGCCCGGGTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCaaaa       <  1:1963389/93‑1 (MQ=255)

                                                                                   cgcgcgGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCaaaaa      <  2:2263919/93‑1 (MQ=255)

                                                                                   cgcgcgGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCaaaaa      <  2:834834/93‑1 (MQ=255)

                                                                                      gcgGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTc   >  1:4266082/1‑93 (MQ=255)

                                                                                      ccgGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGCCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTc   <  2:3302753/92‑1 (MQ=255)

                                                                                       cgGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTcc  <  2:2739244/93‑1 (MQ=255)

                                                                                       cgGGCCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTcc  <  1:1569940/93‑1 (MQ=255)

AAAGGGCGGCAGGTATTGATCGCCAAAGTCACGCCAGGCGTGCCGGTGTTAATTGGCGGCACCGATGTGGTATTGCGCGGCGGCGCGCGGACCGATAAAGACTATTTGAAATTGCTCGATACTCGCCCGGCTATTGGCACGGTACTGAACCAGGGCGATTATGAAAATTTCAAAAAGTCC  >  NC_000913/4440223‑4440402

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 18 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`