BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,795,233	T→C	L220P (CTG→CCG)	rfaL →	O‑antigen ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,795,233	0	T	C	96.8%	386.5 / 1.6	126	L220P (CTG→CCG)	rfaL	O‑antigen ligase
Reads supporting (aligned to +/- strand): ref base T (3/1); new base C (81/41); total (84/42)
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.49e-01

Only 100 of 128 total aligned reads displayed.

CTAAGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCTGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTGCTATA  >  NC_000913/3795143‑3795324

cTAAGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGt                                                                                           >  2:4355147/1‑93 (MQ=255)

 tAAGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGtt                                                                                          >  1:2612571/1‑93 (MQ=255)

     aaaaaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTACCt                                                                                      <  2:2995594/93‑1 (MQ=255)

      aaaaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCta                                                                                     <  2:1507241/93‑1 (MQ=255)

      aaaaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCta                                                                                     <  2:1922694/93‑1 (MQ=255)

         atcatcCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCATCCCTACTCCCGTTCCCTATaa                                                                                  <  2:3263773/93‑1 (MQ=255)

            atcCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAAttt                                                                               <  2:4663690/93‑1 (MQ=255)

             tcCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTtg                                                                              <  1:119973/93‑1 (MQ=255)

                ttttttATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTgtgt                                                                           <  1:3563645/93‑1 (MQ=255)

                    ttatttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtgct                                                                       <  1:2070760/93‑1 (MQ=255)

                       tttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCa                                                                    >  2:2420043/1‑93 (MQ=255)

                       tttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCa                                                                    >  1:925710/1‑93 (MQ=255)

                       tttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCa                                                                    <  1:3344831/93‑1 (MQ=255)

                       tttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCa                                                                    >  2:3325776/1‑93 (MQ=255)

                        ttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCAt                                                                   >  2:1644737/1‑93 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 >  1:112471/1‑93 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 >  2:1015255/1‑93 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 <  1:3188564/93‑1 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 <  1:3188565/93‑1 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 <  2:3636263/93‑1 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 >  2:4044600/1‑93 (MQ=255)

                          ttattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAAACCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTa                                                                 >  2:1015254/1‑93 (MQ=255)

                           tattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTaa                                                                <  1:1626586/93‑1 (MQ=255)

                           tattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTaa                                                                <  1:1626588/93‑1 (MQ=255)

                           tattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTaa                                                                >  1:3513135/1‑93 (MQ=255)

                           tattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTaa                                                                >  2:4566740/1‑93 (MQ=255)

                            attaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAAt                                                               >  2:1613056/1‑93 (MQ=255)

                            attaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAAt                                                               >  2:4631815/1‑93 (MQ=255)

                             ttaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATa                                                              >  1:3034436/1‑93 (MQ=255)

                             ttaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATa                                                              >  1:1529823/1‑93 (MQ=255)

                              taAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAg                                                             >  1:3929034/1‑93 (MQ=255)

                               aaaTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAgt                                                            >  1:2388899/1‑92 (MQ=255)

                               aaaTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGc                                                            >  2:1793784/1‑93 (MQ=255)

                                aaTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCt                                                           >  1:3362287/1‑93 (MQ=255)

                                 aTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCtt                                                          >  1:1538254/1‑93 (MQ=255)

                                  tAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCtta                                                         <  2:3146350/93‑1 (MQ=255)

                                  tAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCtta                                                         >  1:1420932/1‑93 (MQ=255)

                                  tAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCtta                                                         >  1:172760/1‑93 (MQ=255)

                                  tAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCtta                                                         <  2:3146349/93‑1 (MQ=255)

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATGGCttat                                                        >  1:1654954/1‑93 (MQ=255)

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttat                                                        >  1:3766512/1‑93 (MQ=255)

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttat                                                        >  2:3685577/1‑93 (MQ=255)

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttat                                                        >  1:3941284/1‑93 (MQ=255)

                                    gTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatt                                                       >  1:3357306/1‑93 (MQ=255)

                                    gTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatt                                                       <  1:3135502/93‑1 (MQ=255)

                                    gTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatt                                                       >  1:630766/1‑93 (MQ=255)

                                    gTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAAACCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatt                                                       <  1:3033324/93‑1 (MQ=255)

                                     ttGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatta                                                      <  2:2534793/93‑1 (MQ=255)

                                     ttGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatta                                                      <  1:3338104/93‑1 (MQ=255)

                                     ttGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatta                                                      <  2:427825/93‑1 (MQ=255)

                                      tGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTtattat                                                     <  2:3668160/93‑1 (MQ=255)

                                           ggTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTAtaataa                                                >  2:2700802/1‑93 (MQ=255)

                                           ggTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTAtaataa                                                >  1:2972035/1‑93 (MQ=255)

                                            gTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATaaa                                               >  1:633921/1‑93 (MQ=255)

                                             tACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAAt                                              >  2:1279161/1‑93 (MQ=255)

                                                tttATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAc                                           >  2:1661975/1‑93 (MQ=255)

                                                 ttATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAcc                                          <  1:545887/93‑1 (MQ=255)

                                                 ttATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAcc                                          <  1:3532379/93‑1 (MQ=255)

                                                 ttATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAcc                                          <  2:3816957/93‑1 (MQ=255)

                                                  tATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAccc                                         >  2:3022146/1‑93 (MQ=255)

                                                   aTGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCa                                        >  2:2112471/1‑93 (MQ=255)

                                                   aTGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCa                                        >  1:4743637/1‑93 (MQ=255)

                                                    tGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCaa                                       >  1:2003869/1‑93 (MQ=255)

                                                    tGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCaa                                       >  1:2150066/1‑93 (MQ=255)

                                                    tGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCaa                                       >  2:3475987/1‑93 (MQ=255)

                                                     gTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAg                                      <  2:660425/93‑1 (MQ=255)

                                                     gTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAg                                      <  2:660424/93‑1 (MQ=255)

                                                      ttCTGGCGCTAACACAAACCAGAGCAACCCTACTCCTGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGa                                     >  2:4015363/1‑93 (MQ=255)

                                                      ttCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGa                                     <  2:3991805/93‑1 (MQ=255)

                                                      ttCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGa                                     >  2:3217484/1‑93 (MQ=255)

                                                        cTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGaaa                                   >  1:2181685/1‑93 (MQ=255)

                                                        cTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGaaa                                   >  1:2928213/1‑93 (MQ=255)

                                                         tGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAAt                                  >  1:3706663/1‑93 (MQ=255)

                                                         tGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAAt                                  >  1:2607929/1‑93 (MQ=255)

                                                         tGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAAt                                  <  2:3501910/93‑1 (MQ=255)

                                                          ggCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAAtt                                 >  2:1569836/1‑93 (MQ=255)

                                                            cgcTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCa                               >  1:2599567/1‑93 (MQ=255)

                                                             gcTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCAc                              >  1:2356531/1‑93 (MQ=255)

                                                             gcTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCAc                              >  2:3445085/1‑93 (MQ=255)

                                                              cTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACt                             >  1:3733624/1‑93 (MQ=255)

                                                               tAACACAAACCAGAGCAACCCTACTCCTGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACtt                            <  1:2409904/93‑1 (MQ=255)

                                                               tAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACtt                            >  1:2448092/1‑93 (MQ=255)

                                                                 acacAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTcc                          >  2:1905087/1‑93 (MQ=255)

                                                                  cacaAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGTTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCct                         <  1:2172637/93‑1 (MQ=255)

                                                                  cacaAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCct                         >  1:141257/1‑93 (MQ=255)

                                                                   acaAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCctc                        <  2:2377236/93‑1 (MQ=255)

                                                                   acaAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCctc                        <  1:1831188/93‑1 (MQ=255)

                                                                    caAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCctct                       <  2:238094/93‑1 (MQ=255)

                                                                     aaaCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTa                      <  1:1470623/93‑1 (MQ=255)

                                                                     aaaCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTa                      >  1:919465/1‑93 (MQ=255)

                                                                      aaCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTAt                     >  2:535903/1‑93 (MQ=255)

                                                                          agagCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGtt                 >  1:141109/1‑93 (MQ=255)

                                                                          agagAAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGtt                 >  2:624462/1‑93 (MQ=255)

                                                                           gagCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTc                >  1:4419455/1‑93 (MQ=255)

                                                                              cagcccTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac             <  1:3466125/90‑1 (MQ=255)

                                                                              cAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATATCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac             >  1:3993082/1‑93 (MQ=255)

                                                                              cAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac             >  1:3993081/1‑93 (MQ=255)

                                                                              cAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac             >  1:4363535/1‑93 (MQ=255)

                                                                                   cTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCCCTATTGTTCTACTAAtt        >  1:4357083/1‑93 (MQ=255)

                                                                                      cTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTGCt     >  2:4307196/1‑93 (MQ=255)

CTAAGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCTGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTGCTATA  >  NC_000913/3795143‑3795324

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