BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,106,363	A→G	K282K (AAA→AAG)	ymdC →	putative synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,106,363	0	A	G	86.3%	638.8 / 87.7	263	K282K (AAA→AAG)	ymdC	putative synthase
Reads supporting (aligned to +/- strand): ref base A (20/16); new base G (111/116); total (131/132)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 269 total aligned reads displayed.

TAAATCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTACCGACAC  >  NC_000913/1106272‑1106454

taAATCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAAc                                                                                            <  2:686677/93‑1 (MQ=255)

 aaaTCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCg                                                                                           <  1:3177898/93‑1 (MQ=255)

 aaaTCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCg                                                                                           >  1:6730537/1‑93 (MQ=255)

 aaaTCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACg                                                                                           >  1:6828898/1‑93 (MQ=255)

  aaTCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCgg                                                                                          >  2:5764005/1‑93 (MQ=255)

  aaTCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCgg                                                                                          >  1:1442483/1‑93 (MQ=255)

    tcatcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCa                                                                                        <  2:2111966/93‑1 (MQ=255)

     catcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCAt                                                                                       <  1:1251678/93‑1 (MQ=255)

     catcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCAt                                                                                       >  2:4289197/1‑93 (MQ=255)

      atcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCAtt                                                                                      >  2:1623505/1‑93 (MQ=255)

      atcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCAtt                                                                                      <  2:6792651/93‑1 (MQ=255)

      atcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCAtt                                                                                      >  1:4697231/1‑93 (MQ=255)

       tcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTc                                                                                     >  1:175188/1‑93 (MQ=255)

       tcTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTc                                                                                     >  2:2852239/1‑93 (MQ=255)

         tGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCAc                                                                                   >  1:101533/1‑93 (MQ=255)

           gTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTg                                                                                 <  1:5201306/93‑1 (MQ=255)

             tGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGtt                                                                               <  2:6121165/93‑1 (MQ=255)

               aTGGAACATTGCCGCTTATCTGGGCGAAGAGACGTTTAGTAAGTGATGATCCGGCGAAAGGGGAGGGGAAGTCAAAGCGGCATTCACTGTTAc                                                                             <  2:748449/93‑1 (MQ=255)

               aTGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTAc                                                                             <  2:1742230/93‑1 (MQ=255)

                  gAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGc                                                                          <  2:7014888/93‑1 (MQ=255)

                   aaCATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCgga                                                                         <  1:4920902/93‑3 (MQ=255)

                   aaCATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCa                                                                         <  1:826136/93‑1 (MQ=255)

                    aCATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAg                                                                        <  1:735682/93‑1 (MQ=255)

                     cATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAgc                                                                       >  1:383354/1‑93 (MQ=255)

                       ttGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAgcgc                                                                     <  1:1415656/93‑1 (MQ=255)

                       ttGCCGCTTATCTGGGCGAAGACACGTTTATTAAGGGATGATCCGGCGAAAGGGGTGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAgcgc                                                                     >  1:1128358/1‑93 (MQ=255)

                         gCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCt                                                                   <  2:4795578/93‑1 (MQ=255)

                            gCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGtt                                                                <  1:1764264/93‑1 (MQ=255)

                            gCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGtt                                                                >  1:4215990/1‑93 (MQ=255)

                            gCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGtt                                                                >  2:4387398/1‑93 (MQ=255)

                             cTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTc                                                               >  1:3037367/1‑93 (MQ=255)

                             cTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTc                                                               <  1:5922348/93‑1 (MQ=255)

                               tATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGa                                                             >  1:913090/1‑93 (MQ=255)

                                 tCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGata                                                           <  1:4028704/93‑1 (MQ=255)

                                 tCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGata                                                           <  1:6304422/93‑1 (MQ=255)

                                   tGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATc                                                         >  1:4439607/1‑93 (MQ=255)

                                    gggCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCa                                                        <  1:3508351/93‑1 (MQ=255)

                                      gCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATg                                                      <  2:2997668/93‑1 (MQ=255)

                                      gCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATg                                                      <  2:1226629/93‑1 (MQ=255)

                                       cGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATgg                                                     <  1:6651762/93‑1 (MQ=255)

                                       cGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATgg                                                     >  2:7375744/1‑93 (MQ=255)

                                        gAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATggg                                                    >  1:3306467/1‑93 (MQ=255)

                                        gAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATggg                                                    >  2:1945909/1‑93 (MQ=255)

                                        gAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATggg                                                    >  1:1987053/1‑93 (MQ=255)

                                         aaGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGc                                                   >  2:3382670/1‑93 (MQ=255)

                                         aaGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGc                                                   >  1:2127536/1‑93 (MQ=255)

                                          aGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCt                                                  >  1:747787/1‑93 (MQ=255)

                                          aGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCt                                                  >  2:7349040/1‑93 (MQ=255)

                                           gACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTc                                                 >  1:2291704/1‑93 (MQ=255)

                                           gACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTc                                                 >  2:4132999/1‑93 (MQ=255)

                                           gACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTc                                                 >  1:1839915/1‑93 (MQ=255)

                                            acacGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCa                                                >  1:2670732/1‑93 (MQ=255)

                                            acacGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCa                                                >  2:3993354/1‑93 (MQ=255)

                                             cacGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCAc                                               >  2:4497520/1‑93 (MQ=255)

                                              acGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCAcc                                              >  1:3658709/1‑93 (MQ=255)

                                              acGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCAcc                                              >  2:277269/1‑93 (MQ=255)

                                              acGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCAcc                                              <  2:2280873/93‑1 (MQ=255)

                                                gTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCa                                            >  1:5078039/1‑93 (MQ=255)

                                                 tttATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAg                                           >  2:3515988/1‑93 (MQ=255)

                                                 tttATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAg                                           >  1:2975629/1‑93 (MQ=255)

                                                 tttATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAg                                           <  1:5315557/93‑1 (MQ=255)

                                                  ttattaAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGt                                          <  2:1536722/93‑1 (MQ=255)

                                                  ttattaAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGt                                          <  1:3720279/93‑1 (MQ=255)

                                                  ttattaAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGt                                          >  2:3810288/1‑93 (MQ=255)

                                                   tattaAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTg                                         >  1:4881710/1‑93 (MQ=255)

                                                     ttaAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGaa                                       <  2:3128580/93‑1 (MQ=255)

                                                      taAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAAc                                      <  1:5698428/93‑1 (MQ=255)

                                                         gTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCa                                   >  2:6262597/1‑93 (MQ=255)

                                                          tgatgaTCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCAt                                  <  2:4900425/93‑1 (MQ=255)

                                                           gatgatCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATc                                 <  2:2589236/93‑1 (MQ=255)

                                                           gatgatCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATc                                 <  2:1350827/93‑1 (MQ=255)

                                                            atgatCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCg                                >  1:5517741/1‑93 (MQ=255)

                                                            atgatCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCACTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCg                                >  2:3082354/1‑93 (MQ=255)

                                                               atCCGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGata                             <  1:5383676/93‑1 (MQ=255)

                                                                  cGGCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGAtatta                          <  2:2382667/93‑1 (MQ=255)

                                                                  cGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGAtatta                          >  2:5675391/1‑93 (MQ=255)

                                                                   ggCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATAGAtattat                         <  2:4298497/93‑1 (MQ=255)

                                                                    gCGAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTAtc                        <  1:4114644/93‑1 (MQ=255)

                                                                      gAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTAtctc                      >  1:6070239/1‑93 (MQ=255)

                                                                      gAAAGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTAtctc                      >  1:3579482/1‑93 (MQ=255)

                                                                       aaaGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATctct                     >  2:3356476/1‑93 (MQ=255)

                                                                         aGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTc                   >  2:3469469/1‑93 (MQ=255)

                                                                         aGGGGAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTc                   >  1:5868404/1‑93 (MQ=255)

                                                                          ggggAGGGCAAGGCAAAGCGGCATTTAACGGTAACGGCGCCCCTGTTTGATATTAAGGGGTCAACCAAGGGACGCATCGATATTTTTTCtttc                  >  1:6180042/1‑91 (MQ=255)

                                                                          ggggAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTcc                  >  2:3188299/1‑93 (MQ=255)

                                                                             gAGGGCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTAt               <  2:4677633/93‑1 (MQ=255)

                                                                               gggCAAGGCAAAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTAttt             >  1:7156348/1‑93 (MQ=255)

                                                                                ggcaaggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTAtttt            <  1:5178973/93‑1 (MQ=255)

                                                                                 gcaaggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTg           >  2:867573/1‑93 (MQ=255)

                                                                                  caaggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGt          >  2:6316263/1‑93 (MQ=255)

                                                                                   aaggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTa         <  2:5038136/93‑1 (MQ=255)

                                                                                   aaggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTa         >  1:6550895/1‑93 (MQ=255)

                                                                                    aggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTAc        >  2:6634090/1‑93 (MQ=255)

                                                                                    aggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTAc        >  1:3379268/1‑93 (MQ=255)

                                                                                    aggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTAc        >  2:1861305/1‑93 (MQ=255)

                                                                                    aggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTAc        >  1:1680642/1‑93 (MQ=255)

                                                                                     ggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTAcc       >  2:5457105/1‑93 (MQ=255)

                                                                                     ggcaaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTAcc       <  1:7318254/93‑1 (MQ=255)

                                                                                       caaAGCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTACCGa     <  1:3928822/93‑1 (MQ=255)

                                                                                         aaaCGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTACCGaca   >  2:5955495/1‑93 (MQ=255)

TAAATCATCTGGTTGATGGAACATTGCCGCTTATCTGGGCGAAGACACGTTTATTAAGTGATGATCCGGCGAAAGGGGAGGGCAAGGCAAAACGGCATTCACTGTTACCGCAGCGCCTGTTCGATATCATGGGCTCACCCAGTGAACGCATCGATATTATCTCTTCCTATTTTGTACCGACAC  >  NC_000913/1106272‑1106454

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