BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,494	C→G	A261G (GCG→GGG)	yjdA →	putative vimentin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,494	0	C	G	86.4%	744.4 / 94.1	250	A261G (GCG→GGG)	yjdA	putative vimentin
Reads supporting (aligned to +/- strand): ref base C (14/20); new base G (102/113); total (116/134)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01

Only 100 of 253 total aligned reads displayed.

AGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGC  >  NC_000913/4325404‑4325583

aGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTcggggg                                                                                         >  2:4277884/1‑93 (MQ=255)

    gACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACt                                                                                     >  2:3748651/1‑93 (MQ=255)

    gACGTTAATGGATACCCCCGGGCCAAATGAAACCGGGCAACCGCATCTGCAAAAAATGCTTAACAAGCAGCTGGCACGCGCCTCGGGGGAACt                                                                                     >  1:911928/1‑93 (MQ=255)

      cGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTgg                                                                                   <  1:1445347/93‑1 (MQ=255)

       gTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGc                                                                                  <  2:4816234/93‑1 (MQ=255)

        ttACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCg                                                                                 >  2:2036896/1‑93 (MQ=255)

          aCTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGt                                                                               <  1:1213038/93‑1 (MQ=255)

           cTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTg                                                                              >  2:7171207/1‑93 (MQ=255)

             ggATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCt                                                                            >  1:6441163/1‑93 (MQ=255)

              gATCCCCCCGGGCAAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTg                                                                           <  2:6923959/93‑1 (MQ=255)

              gATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTg                                                                           <  2:5531493/93‑1 (MQ=255)

               aTACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTgg                                                                          >  1:2322594/1‑93 (MQ=255)

                tACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGa                                                                         >  2:2843089/1‑93 (MQ=255)

                 aCCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGAt                                                                        <  2:7306867/93‑1 (MQ=255)

                    cccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtat                                                                     >  1:1863809/1‑93 (MQ=255)

                    cccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATtat                                                                     >  1:4782341/1‑93 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    <  2:1232092/93‑1 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    >  1:76066/1‑93 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtata                                                                    >  1:1275029/1‑93 (MQ=255)

                     ccGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtann                                                                    >  2:1722775/1‑91 (MQ=255)

                      cNGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATtatan                                                                   <  2:6655003/93‑2 (MQ=255)

                      cGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATAc                                                                   >  2:6788496/1‑93 (MQ=255)

                       gggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACg                                                                  >  1:6761438/1‑93 (MQ=255)

                       gggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACg                                                                  >  1:1669871/1‑93 (MQ=255)

                       gggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACg                                                                  >  2:6535898/1‑93 (MQ=255)

                       gggCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACg                                                                  <  1:5104271/93‑1 (MQ=255)

                          ccAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCaa                                                               >  1:2166261/1‑93 (MQ=255)

                           cAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAAc                                                              >  2:4326889/1‑93 (MQ=255)

                            aaaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACt                                                             <  2:1821431/93‑1 (MQ=255)

                            aaaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACt                                                             >  2:5422645/1‑93 (MQ=255)

                            aaaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACt                                                             >  2:6289942/1‑93 (MQ=255)

                             aaTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTg                                                            >  1:2838766/1‑93 (MQ=255)

                              aTGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGa                                                           <  1:535101/93‑1 (MQ=255)

                               tGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGaa                                                          <  2:3918295/93‑1 (MQ=255)

                               tGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGaa                                                          <  2:3514022/93‑1 (MQ=255)

                                gAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGaaa                                                         <  2:5857803/93‑1 (MQ=255)

                                gAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGaaa                                                         <  2:4863289/93‑1 (MQ=255)

                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  2:867587/93‑1 (MQ=255)

                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  1:1038861/93‑1 (MQ=255)

                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  1:3999178/93‑1 (MQ=255)

                                      gggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtc                                                   <  1:487560/93‑1 (MQ=255)

                                       ggCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtct                                                  <  2:2208202/93‑1 (MQ=255)

                                        gCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGAtctc                                                 <  2:2947208/93‑1 (MQ=255)

                                         cAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTcc                                                >  1:2146047/1‑93 (MQ=255)

                                         cAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGAACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTcc                                                >  1:6809870/1‑93 (MQ=255)

                                          aaCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCg                                               >  1:3414102/1‑93 (MQ=255)

                                          aaCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCg                                               <  1:5841072/93‑1 (MQ=255)

                                           nncGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              <  2:1899523/91‑1 (MQ=255)

                                           aCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              >  1:7140406/1‑93 (MQ=255)

                                           aCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              >  2:629748/1‑93 (MQ=255)

                                           aCCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGa                                              <  2:4157459/93‑1 (MQ=255)

                                            ccGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGAt                                             <  1:7069878/93‑1 (MQ=255)

                                            ccGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTACACGCAACCGAAATCGATCTCCGAt                                             <  2:7359021/93‑1 (MQ=255)

                                             cGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCACCGATg                                            >  2:7154052/1‑93 (MQ=255)

                                              gCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGa                                           >  1:2215098/1‑93 (MQ=255)

                                               cATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGaa                                          >  2:3065279/1‑93 (MQ=255)

                                               cATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGaa                                          <  2:5050833/93‑1 (MQ=255)

                                               cATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGaa                                          <  1:7332433/93‑1 (MQ=255)

                                                aTCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAag                                         >  2:5829158/1‑93 (MQ=255)

                                                aTCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAag                                         >  2:2900257/1‑93 (MQ=255)

                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAaga                                        >  2:5099240/1‑93 (MQ=255)

                                                 tCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACACAACTGAAATCGATCTCCGATGAaga                                        >  2:1653485/1‑93 (MQ=255)

                                                  cTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAagag                                       >  1:2487238/1‑93 (MQ=255)

                                                   tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTGCTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAgg                                      <  2:1134474/93‑1 (MQ=255)

                                                   tGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAgg                                      >  2:2584670/1‑93 (MQ=255)

                                                    gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGt                                     <  1:7223698/93‑1 (MQ=255)

                                                    gCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGt                                     >  2:2095976/1‑93 (MQ=255)

                                                       aaaaaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCg                                  <  1:3552893/93‑1 (MQ=255)

                                                          aaTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGa                               <  2:7461793/93‑1 (MQ=255)

                                                          aaTGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGa                               >  2:5275534/1‑93 (MQ=255)

                                                           aTGCTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAg                              <  1:4422242/93‑1 (MQ=255)

                                                              cTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCg                           >  2:2693247/1‑93 (MQ=255)

                                                              cTTAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCg                           <  2:3215824/93‑1 (MQ=255)

                                                                tAACCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAt                         >  2:2373809/1‑93 (MQ=255)

                                                                 aaCCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAtt                        <  2:6727026/93‑1 (MQ=255)

                                                                  aCCAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAttt                       <  2:6016719/93‑1 (MQ=255)

                                                                   ccAGCAGCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGAtttt                      <  1:1105948/93‑1 (MQ=255)

                                                                    cagcagCTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTg                     <  1:5948649/93‑1 (MQ=255)

                                                                      gcagcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGc                   <  1:1594268/93‑1 (MQ=255)

                                                                      gcagcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGc                   >  2:3247381/1‑93 (MQ=255)

                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCgg                 <  1:7348368/93‑1 (MQ=255)

                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCgg                 <  2:4031709/93‑1 (MQ=255)

                                                                        agcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCgg                 <  1:3759804/93‑1 (MQ=255)

                                                                         gcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGt                <  1:6559975/93‑1 (MQ=255)

                                                                         gcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGt                >  2:3829770/1‑93 (MQ=255)

                                                                         gcTGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGt                <  1:6375018/93‑1 (MQ=255)

                                                                           tGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTgg              <  1:2997791/93‑1 (MQ=255)

                                                                           tGGCACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTgg              >  1:2004384/1‑93 (MQ=255)

                                                                           tGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTgg              <  1:4340874/93‑1 (MQ=255)

                                                                              cACGCGCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGc           <  1:1891088/93‑1 (MQ=255)

                                                                                cgcgCCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCaa         <  2:7528976/93‑1 (MQ=255)

                                                                                cgcgCCTCGGGGGTACTGACGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCaa         >  2:5339094/1‑93 (MQ=255)

                                                                                 gcgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAAt        <  2:3673470/93‑1 (MQ=255)

                                                                                 gcgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAAt        <  1:4375810/93‑1 (MQ=255)

                                                                                  cgcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATc       <  2:1577646/93‑1 (MQ=255)

                                                                                  cgcCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATc       <  1:3885513/93‑1 (MQ=255)

                                                                                   gcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCg      <  1:4624434/93‑1 (MQ=255)

                                                                                   gcCTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCg      >  2:7008231/1‑93 (MQ=255)

                                                                                    ccTCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCgg     <  1:4281291/93‑1 (MQ=255)

                                                                                      tCGGGGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTg   >  1:4043610/1‑93 (MQ=255)

AGTTGACGTTACTGGATACCCCCGGGCCAAATGAAGCCGGGCAACCGCATCTGCAAAAAATGCTTAACCAGCAGCTGGCACGCGCCTCGGCGGTACTGGCGGTGCTGGATTATACGCAACTGAAATCGATCTCCGATGAAGAGGTCCGTGAGGCGATTTTGGCGGTGGGGCAATCGGTGC  >  NC_000913/4325404‑4325583

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