BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	719,820	A→G	L5L (TTA→TTG)	ybfK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	719,820	0	A	G	99.1%	1132.2 / ‑3.0	328	L5L (TTA→TTG)	ybfK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (173/152); total (173/155)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 330 total aligned reads displayed.

TGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGA  >  NC_000913/719731‑719909

tGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCa                                                                                        >  2:8439875/1‑93 (MQ=255)

tGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCa                                                                                        <  2:6440040/93‑1 (MQ=255)

   gaAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTc                                                                                     >  2:6238032/1‑93 (MQ=255)

   gaAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTc                                                                                     <  1:8662086/93‑1 (MQ=255)

     aaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAg                                                                                   >  1:3927684/1‑93 (MQ=255)

     aaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAg                                                                                   >  2:2801280/1‑93 (MQ=255)

     aaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAg                                                                                   >  2:415670/1‑93 (MQ=255)

      atcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAgt                                                                                  >  1:7296737/1‑92 (MQ=255)

      atcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAgg                                                                                  >  2:5631805/1‑93 (MQ=255)

       tcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGt                                                                                 >  2:3587605/1‑93 (MQ=255)

       acgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGt                                                                                 <  1:6708984/92‑1 (MQ=255)

         gatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTgg                                                                               >  1:4397361/1‑93 (MQ=255)

         gatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTgg                                                                               <  1:806095/93‑1 (MQ=255)

         gatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTgg                                                                               >  1:852391/1‑93 (MQ=255)

          atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGc                                                                              >  2:6076111/1‑93 (MQ=255)

          atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGc                                                                              >  2:3116934/1‑93 (MQ=255)

          atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGc                                                                              >  1:712613/1‑93 (MQ=255)

          atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGc                                                                              >  1:1410810/1‑93 (MQ=255)

          atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGc                                                                              >  2:1663048/1‑93 (MQ=255)

             gTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGa                                                                           >  2:7113891/1‑93 (MQ=255)

               gCCGAAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGata                                                                         <  2:6402455/93‑1 (MQ=255)

                ccGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGataa                                                                        <  2:6716539/93‑1 (MQ=255)

                    aaaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACa                                                                    >  1:1798265/1‑93 (MQ=255)

                    aaaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACa                                                                    <  2:1988506/93‑1 (MQ=255)

                     aaaaaGCAAAGTTACCCCGCTTATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACaa                                                                   <  1:4325640/93‑1 (MQ=255)

                     aaaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACaa                                                                   >  1:995164/1‑93 (MQ=255)

                     aaaaaGCAAAGCTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACaa                                                                   <  1:4295788/93‑1 (MQ=255)

                      aaaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAg                                                                  <  2:7899715/93‑1 (MQ=255)

                       aaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAgg                                                                 >  1:2595954/1‑93 (MQ=255)

                       aaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAgg                                                                 >  2:2563241/1‑93 (MQ=255)

                       aaaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAgg                                                                 <  1:7847871/93‑1 (MQ=255)

                          gCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGtt                                                              >  1:3840140/1‑93 (MQ=255)

                           cAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTa                                                             <  2:6841311/93‑1 (MQ=255)

                           cAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTa                                                             <  2:8603524/93‑1 (MQ=255)

                           cAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTa                                                             <  1:6067494/93‑1 (MQ=255)

                            aaaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTaa                                                            <  2:7302947/93‑1 (MQ=255)

                             aaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                                           <  2:4916403/93‑1 (MQ=255)

                             aaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                                           <  1:5983366/93‑1 (MQ=255)

                             aaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAAc                                                           <  1:3348530/93‑1 (MQ=255)

                              aGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACa                                                          <  2:8629287/93‑1 (MQ=255)

                              aGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACa                                                          >  1:6385688/1‑93 (MQ=255)

                                  aCCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAgaga                                                      >  1:779382/1‑93 (MQ=255)

                                   ccccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGaa                                                     >  1:6953321/1‑93 (MQ=255)

                                   ccccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGaa                                                     <  2:3221403/93‑1 (MQ=255)

                                   ccccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGaa                                                     >  1:4776545/1‑93 (MQ=255)

                                    cccGCTGATAAATAAATCAGGTAGATTGGTTTTTTCAGAATGCAGAAAACCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAAt                                                    >  1:5591631/1‑93 (MQ=255)

                                    cccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAAt                                                    >  1:7651331/1‑93 (MQ=255)

                                     ccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATg                                                   >  1:1014443/1‑93 (MQ=255)

                                     ccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATg                                                   >  1:3465399/1‑93 (MQ=255)

                                     ccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATg                                                   <  2:6013836/93‑1 (MQ=255)

                                     ccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATg                                                   <  2:8549622/93‑1 (MQ=255)

                                     ccGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATg                                                   >  2:3988330/1‑93 (MQ=255)

                                      cGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGa                                                  >  2:839680/1‑93 (MQ=255)

                                      cGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGa                                                  <  1:2068624/93‑1 (MQ=255)

                                       gCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGTTGGCTGATAATACAAGGGTTAACAGAGAATGAg                                                 <  2:2885258/93‑1 (MQ=255)

                                       gCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACATAGAATGAg                                                 <  2:7161925/93‑1 (MQ=255)

                                       gCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAg                                                 >  1:3703580/1‑93 (MQ=255)

                                         tGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCa                                               >  1:5312553/1‑93 (MQ=255)

                                          gATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCAc                                              >  2:5111320/1‑93 (MQ=255)

                                           ttaaataaaTCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                             >  2:5319146/2‑93 (MQ=255)

                                            taaataaaTCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTg                                            >  1:3050713/1‑93 (MQ=255)

                                             aaataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGc                                           >  2:43328/1‑93 (MQ=255)

                                             aaataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGc                                           >  2:4553801/1‑93 (MQ=255)

                                              aataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCg                                          >  2:1410762/1‑93 (MQ=255)

                                               ataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCgg                                         >  1:6446110/1‑93 (MQ=255)

                                                taaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGa                                        <  1:4572536/93‑1 (MQ=255)

                                                 aaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                       >  1:5002124/1‑93 (MQ=255)

                                                 aaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                       >  2:6981676/1‑93 (MQ=255)

                                                 aaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                       <  2:1696812/93‑1 (MQ=255)

                                                     cAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACa                                   <  2:539118/93‑1 (MQ=255)

                                                       ggTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTACCAGAGAATCAGCACTGCGGAATTCAgg                                 <  2:3508028/93‑1 (MQ=255)

                                                        gTAGATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGc                                >  1:4337769/1‑93 (MQ=255)

                                                           gATTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                             <  1:1656482/93‑1 (MQ=255)

                                                            aTTGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCatat                            >  2:237/1‑93 (MQ=255)

                                                             ttGCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATg                           >  2:6557347/1‑93 (MQ=255)

                                                               gCTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGcc                         >  2:747451/1‑93 (MQ=255)

                                                                cTTTTTTCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCt                        <  2:4394363/93‑1 (MQ=255)

                                                                 ttttttCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTg                       >  2:7347375/1‑93 (MQ=255)

                                                                  tttttCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGt                      <  2:1433394/93‑1 (MQ=255)

                                                                    tttCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTaa                    >  1:4101031/1‑93 (MQ=255)

                                                                    tttCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTaa                    <  1:1698124/93‑1 (MQ=255)

                                                                     ttCAGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAAt                   >  2:7750502/1‑93 (MQ=255)

                                                                        aGAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACa                <  2:3859188/93‑1 (MQ=255)

                                                                         gAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACag               >  1:2185352/1‑93 (MQ=255)

                                                                         gAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACag               >  1:1215563/1‑93 (MQ=255)

                                                                         gAATGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACag               <  2:2438507/93‑1 (MQ=255)

                                                                           aTGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagag             >  1:3086402/1‑93 (MQ=255)

                                                                           aTGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagag             >  1:28799/1‑93 (MQ=255)

                                                                            tGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga            <  1:4248098/93‑1 (MQ=255)

                                                                            tGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga            >  1:8298364/1‑93 (MQ=255)

                                                                            tGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga            <  2:8151822/93‑1 (MQ=255)

                                                                            tGCAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACagaga            <  2:2183838/93‑1 (MQ=255)

                                                                              cAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGt          <  1:684881/93‑1 (MQ=255)

                                                                              cAGAAACCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGt          <  1:247944/93‑1 (MQ=255)

                                                                                 aaaCCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAg       <  1:6675835/93‑1 (MQ=255)

                                                                                  aaCCCTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAgg      >  1:2227159/1‑93 (MQ=255)

                                                                                    cccTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGaa    >  2:5709190/1‑93 (MQ=255)

                                                                                    cccTTGCAAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGaa    >  1:6620036/1‑93 (MQ=255)

                                                                                     ccTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAg   >  2:1121659/1‑93 (MQ=255)

                                                                                      cTTGCCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGa  <  2:2735858/93‑1 (MQ=255)

TGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGA  >  NC_000913/719731‑719909

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