BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,011,688	G→A	P146P (CCG→CCA)	fliF →	FliF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,011,688	0	G	A	96.4%	371.2 / ‑0.1	112	P146P (CCG→CCA)	fliF	FliF
Reads supporting (aligned to +/- strand): ref base G (2/2); new base A (47/61); total (49/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

Only 100 of 114 total aligned reads displayed.

AGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAA  >  NC_000913/2011600‑2011777

aGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaa                                                                                       <  1:3363915/93‑1 (MQ=255)

 gTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaaa                                                                                      >  2:4425374/1‑93 (MQ=255)

 gTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaaa                                                                                      >  2:1167936/1‑93 (MQ=255)

  tttGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTaaaa                                                                                     <  1:2814986/93‑1 (MQ=255)

      gTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAgggg                                                                                 <  1:1193276/93‑1 (MQ=255)

      gTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAgggg                                                                                 <  1:1218937/93‑1 (MQ=255)

       tATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgc                                                                                >  1:290231/1‑93 (MQ=255)

       tATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgc                                                                                <  2:3386469/93‑1 (MQ=255)

        aTCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTAACATCGAAACTATCGGCCCAGTAAAAggggtg                                                                               >  2:2513079/1‑91 (MQ=255)

         tCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgcgc                                                                              >  1:1923051/1‑93 (MQ=255)

           aGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgcgcgc                                                                            >  2:1365053/1‑93 (MQ=255)

           aGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGgcgcgc                                                                            >  1:1075623/1‑93 (MQ=255)

            gCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGcgcgcg                                                                           >  2:3652406/1‑93 (MQ=255)

             ccAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGt                                                                          <  2:3164153/93‑1 (MQ=255)

              cAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACGATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTa                                                                         <  2:1540279/93‑1 (MQ=255)

               aGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTAc                                                                        <  2:1688080/93‑1 (MQ=255)

                gTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACa                                                                       <  1:3218479/93‑1 (MQ=255)

                gTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACa                                                                       >  1:2015604/1‑93 (MQ=255)

                   cAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCt                                                                    >  2:1785478/1‑93 (MQ=255)

                   cAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCt                                                                    <  1:4598482/93‑1 (MQ=255)

                    aGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTg                                                                   <  2:3267637/93‑1 (MQ=255)

                     gCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTgg                                                                  <  2:1211991/93‑1 (MQ=255)

                      cGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGc                                                                 <  2:4183920/93‑1 (MQ=255)

                       gAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCa                                                                >  2:2271163/1‑93 (MQ=255)

                       gAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCa                                                                >  1:4015593/1‑93 (MQ=255)

                        aaCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCaa                                                               >  2:10589/1‑93 (MQ=255)

                         aCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                              <  1:3101659/93‑1 (MQ=255)

                         aCAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAt                                                              <  2:3888521/93‑1 (MQ=255)

                           aGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGc                                                            <  1:4421187/93‑1 (MQ=255)

                            ggTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGcc                                                           >  2:3541036/1‑93 (MQ=255)

                            ggTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGcc                                                           >  1:2475466/1‑93 (MQ=255)

                             gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTGCCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                          <  1:294124/93‑1 (MQ=255)

                             gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                          <  1:3719440/93‑1 (MQ=255)

                             gTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCg                                                          >  2:550492/1‑93 (MQ=255)

                                aaTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGaaa                                                       >  1:3422303/1‑93 (MQ=255)

                                aaTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGaaa                                                       >  1:1605872/1‑93 (MQ=255)

                                 aTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAAc                                                      <  1:882794/93‑1 (MQ=255)

                                 aTTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAAc                                                      <  1:3547984/93‑1 (MQ=255)

                                  ttATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAAcc                                                     >  2:4398406/1‑93 (MQ=255)

                                   tATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCg                                                    <  2:61492/93‑1 (MQ=255)

                                    aTCAGGGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGt                                                   <  1:1681569/93‑1 (MQ=255)

                                    aTCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTATCATCGAAACTATCGGGCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGGAACCGt                                                   >  1:1361221/1‑93 (MQ=255)

                                    aTCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGt                                                   >  1:165418/1‑93 (MQ=255)

                                    aTCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGt                                                   >  2:741738/1‑93 (MQ=255)

                                     tCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTc                                                  <  2:530125/93‑1 (MQ=255)

                                      cAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCt                                                 <  1:2894775/93‑1 (MQ=255)

                                       aGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCtt                                                >  1:3244050/1‑93 (MQ=255)

                                          gggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGTGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAt                                             <  2:3281153/93‑1 (MQ=255)

                                          gggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAt                                             <  2:938887/93‑1 (MQ=255)

                                           ggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAtt                                            >  1:259724/1‑93 (MQ=255)

                                           ggCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTAtt                                            <  1:477521/93‑1 (MQ=255)

                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       <  2:4476889/93‑1 (MQ=255)

                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       <  2:2305975/93‑1 (MQ=255)

                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       >  2:3577790/1‑93 (MQ=255)

                                                tGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTcc                                       <  1:180638/93‑1 (MQ=255)

                                                 ggAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCg                                      >  1:746646/1‑93 (MQ=255)

                                                 ggAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCg                                      >  2:4221644/1‑93 (MQ=255)

                                                    aGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGa                                   <  2:2617685/93‑1 (MQ=255)

                                                     ggCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGaa                                  <  2:2785281/93‑1 (MQ=255)

                                                     ggCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGaa                                  <  2:1426310/93‑1 (MQ=255)

                                                     ggCGAGCTTTCTCGTACCATCGAAACTATCGGACCAGTTAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGaa                                  >  1:391926/1‑93 (MQ=255)

                                                      gCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTTTGTGAAc                                 >  1:3726094/1‑93 (MQ=255)

                                                      gCGAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAAc                                 <  2:183956/93‑1 (MQ=255)

                                                        gAGCTTTCTCGTACCATCGAACCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCTTCTTTATTCGTCCGTGAACaa                               <  2:697823/93‑1 (MQ=255)

                                                        gAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaa                               <  1:3933348/93‑1 (MQ=255)

                                                        gAGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaa                               <  1:2951265/93‑1 (MQ=255)

                                                         aGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaa                              >  1:64195/1‑93 (MQ=255)

                                                         aGCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaa                              >  2:3619181/1‑93 (MQ=255)

                                                          gCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaa                             <  2:2894445/93‑1 (MQ=255)

                                                          gCTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaa                             >  1:344586/1‑93 (MQ=255)

                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            >  1:4412650/1‑93 (MQ=255)

                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            <  1:1156457/93‑1 (MQ=255)

                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            <  1:2155962/93‑1 (MQ=255)

                                                           cTTTCTCGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACaaaaa                            >  1:2688709/1‑93 (MQ=255)

                                                              tcncGTACCATCGAAACTATCNNCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAccc                         <  1:2302048/93‑1 (MQ=255)

                                                               ctcGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAAcccc                        <  1:1369124/93‑1 (MQ=255)

                                                                 cGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCt                      <  2:3424677/93‑1 (MQ=255)

                                                                 cGTACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCt                      <  1:3170664/93‑1 (MQ=255)

                                                                   tACCATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTc                    >  1:1106437/1‑93 (MQ=255)

                                                                     ccATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCtc                  >  1:1828666/1‑92 (MQ=255)

                                                                      cATCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGc                 <  2:4466205/93‑1 (MQ=255)

                                                                       aTCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCa                >  1:818780/1‑93 (MQ=255)

                                                                        tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt               <  2:4310177/93‑1 (MQ=255)

                                                                        tCGAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCAt               >  1:831938/1‑93 (MQ=255)

                                                                         cGAATCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATc              <  2:347246/93‑1 (MQ=255)

                                                                          gAAACTATCGGCCCGGTAAAAGGGGCGGGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCg             <  1:2244686/93‑1 (MQ=255)

                                                                          gAAACTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCg             >  1:4300811/1‑93 (MQ=255)

                                                                           aaaCTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCgg            >  1:4207220/1‑93 (MQ=255)

                                                                            aaCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGt           >  2:3230393/1‑93 (MQ=255)

                                                                            aaCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAAGGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAGCCCCCTTCTGCATCGGt           <  1:2202901/93‑1 (MQ=255)

                                                                             aCTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTg          <  1:744145/93‑1 (MQ=255)

                                                                              cTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGa         >  1:4043376/1‑93 (MQ=255)

                                                                              cTATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGa         <  1:2455307/93‑1 (MQ=255)

                                                                               tATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGAc        >  2:2661414/1‑93 (MQ=255)

                                                                               tATCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGAc        <  2:977233/93‑1 (MQ=255)

                                                                                aTCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACg       <  1:135750/93‑1 (MQ=255)

                                                                                aTCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACACCCTTCTGCATCGGTGACg       >  2:2871216/1‑93 (MQ=255)

                                                                                 tCGGCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACgg      >  2:2447307/1‑93 (MQ=255)

                                                                                    gCCCAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTaa   >  1:1592289/1‑93 (MQ=255)

                                                                                     cccAGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAACCCCCTTCTGCATCGGTGACGGTaaa  <  2:1955411/93‑1 (MQ=255)

AGTTTGGTATCAGCCAGTTCAGCGAACAGGTGAATTATCAGCGGGCGCTGGAAGGCGAGCTTTCTCGTACCATCGAAACTATCGGCCCGGTAAAAGGGGCGCGCGTACATCTGGCAATGCCGAAACCGTCTTTATTCGTCCGTGAACAAAAATCCCCTTCTGCATCGGTGACGGTAAA  >  NC_000913/2011600‑2011777

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