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	99.4%	579.1 / ‑4.6	179	L220P (CTG→CCG)	rfaL	O‑antigen ligase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (97/81); total (98/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 182 total aligned reads displayed.

CTAAGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCTGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTGCTA  >  NC_000913/3795143‑3795322

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

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

  aaGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCTGTTc                                                                                       >  1:1553631/1‑93 (MQ=255)

      aaaaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTGCTCCCGTTCCCta                                                                                   <  1:6438235/93‑1 (MQ=255)

      aaaaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCta                                                                                   >  2:3969406/1‑93 (MQ=255)

       aaaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCtat                                                                                  <  2:5599875/93‑1 (MQ=255)

        aaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCtata                                                                                 >  1:1404985/1‑93 (MQ=255)

        aaTCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCtata                                                                                 <  1:4504795/93‑1 (MQ=255)

          tcatcCTTTTTTATTTTTATTGAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAAt                                                                               >  1:5977525/1‑93 (MQ=255)

          tcatcCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAAt                                                                               <  1:5690039/93‑1 (MQ=255)

             tcCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTtg                                                                            >  2:934026/1‑93 (MQ=255)

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

             tcCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTtg                                                                            >  1:5331725/1‑93 (MQ=255)

                ttttttATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTgtgt                                                                         >  2:5636083/1‑93 (MQ=255)

                  ttttatttttatTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAATCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtg                                                                       <  1:4181260/93‑1 (MQ=255)

                  ttttatttttatTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtg                                                                       >  1:2033968/1‑93 (MQ=255)

                  ttttatttttatTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtg                                                                       >  2:4636494/1‑93 (MQ=255)

                     tatttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtgctg                                                                    <  1:2483451/93‑1 (MQ=255)

                      atttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtgctgc                                                                   >  1:3074410/1‑93 (MQ=255)

                      atttttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTtgctgc                                                                   >  1:340380/1‑93 (MQ=255)

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

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

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

                         tttattAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCAtt                                                                <  2:2423316/93‑1 (MQ=255)

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

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

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

                           tattaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTCGTGTTGCTGCATTaa                                                              >  2:6671014/1‑93 (MQ=255)

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

                             ttaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATa                                                            <  2:4446566/93‑1 (MQ=255)

                             ttaAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATa                                                            <  2:459675/93‑1 (MQ=255)

                              taAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAg                                                           <  2:2701360/93‑1 (MQ=255)

                               aaaTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGc                                                          <  1:3647539/93‑1 (MQ=255)

                                aaTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCt                                                         <  1:3507708/93‑1 (MQ=255)

                                 aTAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCtt                                                        <  2:1849822/93‑1 (MQ=255)

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

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

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

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttat                                                      <  2:480424/93‑1 (MQ=255)

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttat                                                      <  2:3871062/93‑1 (MQ=255)

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

                                   aGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttat                                                      <  2:4544041/93‑1 (MQ=255)

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

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

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

                                     ttGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatta                                                    >  2:6072032/1‑93 (MQ=255)

                                     ttGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCttatta                                                    >  2:6436946/1‑93 (MQ=255)

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

                                        cgcgGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTataa                                                 <  2:4257476/93‑1 (MQ=255)

                                         gcgGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTataat                                                >  1:802495/1‑93 (MQ=255)

                                          cgGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTataata                                               >  1:61250/1‑93 (MQ=255)

                                          cgGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTataata                                               >  2:1793498/1‑93 (MQ=255)

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

                                              aCTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATc                                           >  1:582651/1‑93 (MQ=255)

                                              aCTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATc                                           >  1:5464006/1‑93 (MQ=255)

                                               cTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCa                                          <  1:1282457/93‑1 (MQ=255)

                                               cTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCa                                          >  2:790291/1‑93 (MQ=255)

                                                tttATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAc                                         <  1:4711613/93‑1 (MQ=255)

                                                  tATGTTCTGGCGCTAACACAAACCAGAGCATCCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAccc                                       <  1:560958/93‑1 (MQ=255)

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

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

                                                  tATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCAccc                                       <  1:5294472/93‑1 (MQ=255)

                                                   aTGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCa                                      <  2:4991714/93‑1 (MQ=255)

                                                   aTGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCa                                      <  1:2135768/93‑1 (MQ=255)

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

                                                     gTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAg                                    >  2:1984566/1‑93 (MQ=255)

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

                                                      ttCTGGCGCTAACACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACACAAGa                                   <  1:4137307/93‑1 (MQ=255)

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

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

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

                                                                aaCACAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTc                         >  1:5326708/1‑93 (MQ=255)

                                                                 acacAAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTcc                        <  1:3775626/93‑1 (MQ=255)

                                                                  cacaAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCct                       <  2:2902600/93‑1 (MQ=255)

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

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

                                                                    caAACCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCctct                     >  2:323524/1‑93 (MQ=255)

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

                                                                      aaCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTAt                   <  2:4044048/93‑1 (MQ=255)

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

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

                                                                      aaCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTAt                   <  1:3913890/93‑1 (MQ=255)

                                                                       aCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTAtt                  >  1:2716272/1‑93 (MQ=255)

                                                                       aCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTAtt                  >  1:1563966/1‑93 (MQ=255)

                                                                       aCCAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTAtt                  >  2:1380230/1‑93 (MQ=255)

                                                                        ccAGAGCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTg                 <  1:2529859/93‑1 (MQ=255)

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

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

                                                                           gagCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTc              <  2:3669044/93‑1 (MQ=255)

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

                                                                            agCAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTct             >  1:3862048/1‑93 (MQ=255)

                                                                              cAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac           <  1:638448/93‑1 (MQ=255)

                                                                              cAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac           <  1:1867291/93‑1 (MQ=255)

                                                                              cAACCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctac           <  1:1238842/93‑1 (MQ=255)

                                                                                aCCCTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTctacta         >  2:4223356/1‑93 (MQ=255)

                                                                                  ccTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAAt       >  1:1169528/1‑93 (MQ=255)

                                                                                   cTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAAtt      <  2:3299600/93‑1 (MQ=255)

                                                                                   cTACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAAtt      <  1:5445877/93‑1 (MQ=255)

                                                                                    tACTCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTg     >  2:6546218/1‑93 (MQ=255)

                                                                                       tCCCGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTGCta  <  1:3928428/93‑1 (MQ=255)

CTAAGAAAAATCATCCTTTTTTATTTTTATTAAATAGTTGCGCGGTACTTTATGTTCTGGCGCTAACACAAACCAGAGCAACCCTACTCCTGTTCCCTATAATTTGTGTTGCTGCATTAATAGCTTATTATAATAAATCACCCAAGAAATTCACTTCCTCTATTGTTCTACTAATTGCTA  >  NC_000913/3795143‑3795322

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