BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	728,069	T→C	intergenic (‑25/‑288)	kdpF ← / → ybfA	KdpF/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	728,069	0	T	C	100.0%	407.2 / NA	130	intergenic (‑25/‑288)	kdpF/ybfA	KdpF/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (56/74); total (56/74)

Only 100 of 133 total aligned reads displayed.

GCATAAACCAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGTGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGCCAT  >  NC_000913/727980‑728158

gCATAAACCAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGAt                                                                                        <  1:2226876/93‑1 (MQ=255)

  aTAAACCAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATAc                                                                                      <  2:1743629/93‑1 (MQ=255)

      aCCAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGa                                                                                  >  2:1694946/1‑93 (MQ=255)

      aCCAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGa                                                                                  <  1:1849119/93‑1 (MQ=255)

       ccAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGaa                                                                                 <  1:3404948/93‑1 (MQ=255)

       ccAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGaa                                                                                 <  2:1877493/93‑1 (MQ=255)

        cAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAg                                                                                <  2:880152/93‑1 (MQ=255)

        cAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAg                                                                                >  1:864737/1‑93 (MQ=255)

        cAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAg                                                                                >  1:1680070/1‑93 (MQ=255)

          gATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAggg                                                                              <  1:3183590/93‑1 (MQ=255)

              aCCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGa                                                                          >  1:4416866/1‑93 (MQ=255)

               cccAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGaa                                                                         <  1:4176435/93‑1 (MQ=255)

                ccAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAAt                                                                        >  1:1461603/1‑93 (MQ=255)

                ccAGAAAAAAAACCACAAGAAATCCG‑CGGTTAACCCGCCCTGCACCCAAAG‑GCCCCCCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAAtt                                                                       <  2:252236/93‑1 (MQ=255)

                 cAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAAtt                                                                       <  1:2258536/93‑1 (MQ=255)

                 cAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAAtt                                                                       >  1:3029842/1‑93 (MQ=255)

                  aGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTc                                                                      <  1:994396/93‑1 (MQ=255)

                  aGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTc                                                                      >  1:795918/1‑93 (MQ=255)

                    taataaAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTg                                                                    <  1:4441460/93‑1 (MQ=255)

                    taataaAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTg                                                                    <  1:2758338/93‑1 (MQ=255)

                      ataaAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGc                                                                  <  1:1110461/93‑1 (MQ=255)

                       taaAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCt                                                                 >  1:2419171/1‑93 (MQ=255)

                       taaAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCt                                                                 >  2:1731430/1‑93 (MQ=255)

                       taaAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCt                                                                 >  2:800345/1‑93 (MQ=255)

                            acacCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGGATTCTGGCTGCaaa                                                            >  2:1636204/1‑93 (MQ=255)

                            acacCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCaaa                                                            <  1:3663931/93‑1 (MQ=255)

                             cacCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGAACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAg                                                           >  2:3825968/1‑93 (MQ=255)

                              acCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGa                                                          <  2:4059225/93‑1 (MQ=255)

                               ccAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCCAAGAt                                                         >  1:2279609/1‑93 (MQ=255)

                               ccAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGAt                                                         <  1:4087116/93‑1 (MQ=255)

                                cAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGAtt                                                        >  2:4002251/1‑93 (MQ=255)

                                cAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGAtt                                                        <  1:1062552/93‑1 (MQ=255)

                                  gCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTc                                                      <  1:2350092/93‑1 (MQ=255)

                                  gCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTc                                                      <  2:1251948/93‑1 (MQ=255)

                                  gCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTc                                                      <  1:4286821/93‑1 (MQ=255)

                                     aTACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCa                                                   <  2:2961018/93‑1 (MQ=255)

                                      tACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCaa                                                  >  2:1062532/1‑93 (MQ=255)

                                      tACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCaa                                                  <  2:2377596/93‑1 (MQ=255)

                                       aCGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCaaa                                                 >  1:3738655/1‑93 (MQ=255)

                                         gCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCaaaaa                                               <  1:3616469/93‑1 (MQ=255)

                                          ccGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAAt                                              >  1:2746972/1‑93 (MQ=255)

                                          ccGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAAt                                              >  2:4428370/1‑93 (MQ=255)

                                           cGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATc                                             <  1:702103/93‑1 (MQ=255)

                                           cGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATc                                             <  1:2707956/93‑1 (MQ=255)

                                            ggTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCt                                            <  1:1191033/93‑1 (MQ=255)

                                             gTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTg                                           >  2:2086927/1‑93 (MQ=255)

                                             gTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTg                                           >  2:891880/1‑93 (MQ=255)

                                              ttATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgc                                          <  1:361649/93‑1 (MQ=255)

                                              ttATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgc                                          <  2:1409925/93‑1 (MQ=255)

                                              ttATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTAGCTGCAAAGATTTCGCAAAAATCTgc                                          >  2:4316163/1‑93 (MQ=255)

                                                aTCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAACAATCTgcgg                                        >  1:3231336/1‑93 (MQ=255)

                                                aTCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgcgg                                        >  2:3291767/1‑93 (MQ=255)

                                                aTCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgcgg                                        >  1:957688/1‑93 (MQ=255)

                                                  cACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgcgggc                                      <  2:4080002/93‑1 (MQ=255)

                                                   aCG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgcgggcg                                     >  2:3446359/1‑93 (MQ=255)

                                                    cG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTgcgggcgg                                    <  2:1951057/93‑1 (MQ=255)

                                                     g‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGcgggcggg                                   <  2:553398/93‑1 (MQ=255)

                                                     g‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGcgggcggg                                   >  2:1793637/1‑93 (MQ=255)

                                                       ccTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCgggg                                  <  2:521670/93‑1 (MQ=255)

                                                        cTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGgt                                 >  1:2912087/1‑93 (MQ=255)

                                                          gCACTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGgtgt                               <  2:4302952/93‑1 (MQ=255)

                                                            aCTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTaa                             <  1:2970202/93‑1 (MQ=255)

                                                            aCTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTaa                             >  2:1488015/1‑93 (MQ=255)

                                                             cTCACAGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTaaa                            <  1:4568067/93‑1 (MQ=255)

                                                                 caGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTaaaaaaa                        <  2:661454/93‑1 (MQ=255)

                                                                 caGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTaaaaaaa                        >  2:818842/1‑93 (MQ=255)

                                                                 caGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTaaaaaaa                        <  2:4042908/93‑1 (MQ=255)

                                                                  aGTGCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAg                       >  1:488340/1‑93 (MQ=255)

                                                                     gCACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGtat                    <  2:4472335/93‑1 (MQ=255)

                                                                      cACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGtata                   <  2:3490699/93‑1 (MQ=255)

                                                                      cACCTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGtata                   <  1:4414740/93‑1 (MQ=255)

                                                                        ccTCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATaaa                 <  2:1945353/93‑1 (MQ=255)

                                                                          tCCAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATaaaaa               <  1:1203843/93‑1 (MQ=255)

                                                                           ccAGTGGCCTAAAAGCGATACCGGAAGGTTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAAt              <  2:1373426/93‑1 (MQ=255)

                                                                           ccAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAAt              <  2:2214144/93‑1 (MQ=255)

                                                                           ccAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAAt              <  1:945563/93‑1 (MQ=255)

                                                                            cAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATg             <  2:2319497/93‑1 (MQ=255)

                                                                            cAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATg             <  2:805531/93‑1 (MQ=255)

                                                                            cAGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATg             <  1:1383697/93‑1 (MQ=255)

                                                                             aGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATgg            <  2:935323/93‑1 (MQ=255)

                                                                             aGTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATgg            >  2:4149687/1‑93 (MQ=255)

                                                                              gTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGc           >  2:3699790/1‑93 (MQ=255)

                                                                              gTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGc           >  1:699577/1‑93 (MQ=255)

                                                                              gTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGc           >  2:4508833/1‑93 (MQ=255)

                                                                              gTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGc           <  2:2656733/93‑1 (MQ=255)

                                                                              gTGGCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCGGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGc           <  1:143475/93‑1 (MQ=255)

                                                                                ggCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCaa         <  2:2789597/93‑1 (MQ=255)

                                                                                ggCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCaa         <  1:860868/93‑1 (MQ=255)

                                                                                 gCCTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGAATAAAAATGGCaaa        >  1:3117476/1‑93 (MQ=255)

                                                                                  ccTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCaaaa       <  1:3793073/93‑1 (MQ=255)

                                                                                  ccTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCaaaa       <  2:2502341/93‑1 (MQ=255)

                                                                                   cTAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCaaaaa      >  1:3850071/1‑92 (MQ=255)

                                                                                    tAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGc     >  1:2338570/1‑93 (MQ=255)

                                                                                    tAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGc     <  1:3817818/93‑1 (MQ=255)

                                                                                    tAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGc     >  1:3742654/1‑93 (MQ=255)

                                                                                    tAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGc     >  1:4673712/1‑93 (MQ=255)

                                                                                    tAAAAGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGc     >  1:895082/1‑93 (MQ=255)

                                                                                      aaaGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGCCa   <  1:1521883/93‑1 (MQ=255)

                                                                                      aaaGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGCCa   <  2:109605/93‑1 (MQ=255)

                                                                                       aaGCGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGCCAt  >  2:1641782/1‑93 (MQ=255)

GCATAAACCAGATAACCCAGTAATAAAAACACCAGCAATACGCCGGTTATCACG‑CCTGCACTCACAGTGCACCTCCAGTGGCCTAAAAGTGATACCGGAAGGGTAGAATTCTGGCTGCAAAGATTTCGCAAAAATCTGCGGGCGGGGTGTAAAAAAAGTATAAAAATGGCAAAAGCCAT  >  NC_000913/727980‑728158

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