BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,041,632	T→A	V107E (GTG→GAG)	ygfZ →	folate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,041,632	0	T	A	97.5%	713.9 / 7.3	198	V107E (GTG→GAG)	ygfZ	folate‑binding protein
Reads supporting (aligned to +/- strand): ref base T (2/3); new base A (96/96); total (99/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 200 total aligned reads displayed.

CTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAA  >  NC_000913/3041535‑3041728

cTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAgaga                                                                                                <  1:3983574/100‑1 (MQ=255)

cTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAgaga                                                                                                <  1:1705243/100‑1 (MQ=255)

   ttCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCa                                                                                             >  1:3820145/1‑100 (MQ=255)

   ttCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCa                                                                                             >  1:815175/1‑100 (MQ=255)

          aCGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGcc                                                                                      <  1:1939290/100‑1 (MQ=255)

           cGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCa                                                                                     <  1:3247007/100‑1 (MQ=255)

            ggCGATGGCTTTGCATGGATTGAACGGCGTAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAg                                                                                    >  1:4930617/1‑100 (MQ=255)

             gCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAga                                                                                   <  1:1165526/100‑1 (MQ=255)

             gCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCCGTATTCTCTAAAGAGACCATCGCGCCAga                                                                                   <  1:4937530/100‑1 (MQ=255)

              cGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCa                                                                                             >  2:1978778/1‑89 (MQ=255)

              cGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCa                                                                                             <  2:1224464/89‑1 (MQ=255)

                  ggCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATcgc                                                                                         >  2:3323317/1‑89 (MQ=255)

                     tttGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGcc                                                                                      <  2:3591079/89‑1 (MQ=255)

                        gCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGGCAGACGACGAGCgtg                                                                        <  1:4726789/100‑1 (MQ=255)

                        gCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgtg                                                                        >  1:1441989/1‑100 (MQ=255)

                        gCATGGATTGAACGGCGCAGCGCGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgtg                                                                        <  1:530311/100‑1 (MQ=255)

                           tGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGtgc                                                                     >  1:698898/1‑100 (MQ=255)

                            ggATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAgacgac                                                                               <  2:411955/89‑1 (MQ=255)

                            ggATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGtgct                                                                    >  1:1583214/1‑100 (MQ=255)

                             gATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGacgacg                                                                              <  2:734514/89‑1 (MQ=255)

                                tGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGc                                                                           >  2:2245877/1‑89 (MQ=255)

                                tGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCtt                                                                <  1:980640/100‑1 (MQ=255)

                                 gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCg                                                                          >  2:1687368/1‑89 (MQ=255)

                                 gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCg                                                                          >  2:4888798/1‑89 (MQ=255)

                                 gAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTg                                                               >  1:4103465/1‑100 (MQ=255)

                                  aaCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTgg                                                              <  1:3194080/100‑1 (MQ=255)

                                   aCGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgt                                                             <  1:2255147/100‑1 (MQ=255)

                                    cGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgtgt                                                                       <  2:3806689/89‑1 (MQ=255)

                                    cGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCgtgt                                                                       >  2:884814/1‑89 (MQ=255)

                                     ggCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGtgtg                                                                      <  2:3266526/89‑1 (MQ=255)

                                     ggCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGtgtg                                                                      >  2:470274/1‑89 (MQ=255)

                                     ggCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgtgt                                                           >  1:2569325/1‑100 (MQ=255)

                                     ggCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCGCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTGGgtgt                                                           >  1:5026655/1‑100 (MQ=255)

                                        gcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGc                                                        <  1:4774911/100‑1 (MQ=255)

                                        gcAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACCAGCGTGtgctg                                                                   >  2:4946355/1‑89 (MQ=255)

                                           gcgtgcgtGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCtt                                                                >  2:1639025/1‑89 (MQ=255)

                                           gcgtgcgtGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCgg                                                     >  1:4169901/1‑100 (MQ=255)

                                              tgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGgt                                                             <  2:212172/89‑1 (MQ=255)

                                              tgcgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGttt                                                  <  1:1216728/100‑1 (MQ=255)

                                                cgtgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                <  1:3074977/100‑1 (MQ=255)

                                                  tgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTg                                                         >  2:4340094/1‑89 (MQ=255)

                                                  tgAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAg                                              <  1:2848214/100‑1 (MQ=255)

                                                   gAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGc                                                        >  2:183689/1‑89 (MQ=255)

                                                   gAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAgg                                             <  1:1382131/100‑1 (MQ=255)

                                                   gAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAgg                                             >  1:1273773/1‑100 (MQ=255)

                                                   gAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAgg                                             >  1:2774013/1‑100 (MQ=255)

                                                     aCCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcg                                           >  1:2527009/1‑100 (MQ=255)

                                                      ccGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCgg                                                     >  2:303893/1‑89 (MQ=255)

                                                      ccGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgc                                          >  1:3639209/1‑100 (MQ=255)

                                                       cGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGt                                                    >  2:237578/1‑89 (MQ=255)

                                                       cGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGt                                                    >  2:4714865/1‑89 (MQ=255)

                                                        gCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGtt                                                   <  2:4634742/89‑1 (MQ=255)

                                                           gCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTc                                                <  2:569477/89‑1 (MQ=255)

                                                           gCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgcc                                     >  1:2729194/1‑100 (MQ=255)

                                                           gCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgcc                                     <  1:1944742/100‑1 (MQ=255)

                                                            ctgactgaACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                               <  2:2431/89‑1 (MQ=255)

                                                            ctgactgaACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCa                                               <  2:4140360/89‑1 (MQ=255)

                                                            ctgactgaACTGAAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccg                                    >  1:2936275/1‑100 (MQ=255)

                                                               actgaactgaAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgc                                            <  2:3138325/89‑1 (MQ=255)

                                                               actgaactgaAAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgcgc                                 >  1:2607104/1‑100 (MQ=255)

                                                                ctgaactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                >  1:2114597/1‑100 (MQ=255)

                                                                ctgaactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                <  1:2741402/100‑1 (MQ=255)

                                                                 tgaactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTg                               >  1:4566316/1‑100 (MQ=255)

                                                                  gaactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTgg                              <  1:4433196/100‑1 (MQ=255)

                                                                  gaactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTgg                              >  1:3994379/1‑100 (MQ=255)

                                                                   aactgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGc                             <  1:632494/100‑1 (MQ=255)

                                                                     ctgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGgcgcgcgc                                      >  2:909788/1‑89 (MQ=255)

                                                                     ctgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                           <  1:4553859/100‑1 (MQ=255)

                                                                     ctgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                           >  1:100896/1‑100 (MQ=255)

                                                                     ctgaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCaa                           >  1:2654919/1‑100 (MQ=255)

                                                                       gaaAAAATATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAt                         >  1:4145471/1‑100 (MQ=255)

                                                                        aaaaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgc                                   >  2:4216077/1‑89 (MQ=255)

                                                                        aaaaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgc                                   >  2:2784305/1‑89 (MQ=255)

                                                                         aaaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGcgcgccgcg                                  >  2:3918016/1‑89 (MQ=255)

                                                                         aaaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtct                       <  1:1837920/100‑1 (MQ=255)

                                                                          aaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtctc                      >  1:3397843/1‑100 (MQ=255)

                                                                          aaaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAtctc                      <  1:1914791/100‑1 (MQ=255)

                                                                           aaaTATGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCt                                <  2:4089331/89‑1 (MQ=255)

                                                                              tatGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGc                             <  2:1479333/89‑1 (MQ=255)

                                                                               atGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAg                 <  1:37547/100‑1 (MQ=255)

                                                                                tGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAgn                <  1:2232973/100‑2 (MQ=255)

                                                                                tGCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGc                <  1:1136712/100‑1 (MQ=255)

                                                                                 gCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCg               <  1:3461403/100‑1 (MQ=255)

                                                                                 gCGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCg               >  1:2407302/1‑100 (MQ=255)

                                                                                  cGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAAt                         >  2:3380508/1‑89 (MQ=255)

                                                                                  cGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa              <  1:178943/100‑1 (MQ=255)

                                                                                  cGGTATTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa              >  1:4044671/1‑100 (MQ=255)

                                                                                      aTTCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTg          >  1:89362/1‑100 (MQ=255)

                                                                                       ttCTCTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCtt                    <  2:4473753/89‑1 (MQ=255)

                                                                                        tctcTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCttt                   <  2:3058797/89‑1 (MQ=255)

                                                                                        tctcTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCttt                   <  2:1407200/89‑1 (MQ=255)

                                                                                        tctcTAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGcc        <  1:4689623/100‑1 (MQ=255)

                                                                                         ctctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCt       >  1:4988111/1‑100 (MQ=255)

                                                                                          tctAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCtt      <  1:833882/100‑1 (MQ=255)

                                                                                            tAAAGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGTGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCg    <  1:5038186/100‑1 (MQ=255)

                                                                                             aaaGAGACCATCGCGCCAGACGCCGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa              <  2:1533762/89‑1 (MQ=255)

                                                                                             aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa              >  2:4166153/1‑89 (MQ=255)

                                                                                             aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa              >  2:784087/1‑89 (MQ=255)

                                                                                             aaaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGa              >  2:992571/1‑89 (MQ=255)

                                                                                              aaGAGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCgaa  <  1:3312618/100‑1 (MQ=255)

CTGTTCCGCGACGGCGATGGCTTTGCATGGATTGAACGGCGCAGCGTGCGTGAACCGCAGCTGACTGAACTGAAAAAATATGCGGTATTCTCTAAAGTGACCATCGCGCCAGACGACGAGCGTGTGCTGCTTGGTGTTGCCGGTTTTCAGGCGCGCGCCGCGCTGGCAAATCTCTTTAGCGAACTGCCTTCGAA  >  NC_000913/3041535‑3041728

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