BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,568,665	A→G	L160P (CTC→CCC)	yjiR ←	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,568,665	0	A	G	87.6%	722.8 / 86.3	298	L160P (CTC→CCC)	yjiR	putative regulator
Reads supporting (aligned to +/- strand): ref base A (18/19); new base G (115/146); total (133/165)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

Only 100 of 299 total aligned reads displayed.

GTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACC  >  NC_000913/4568573‑4568754

  ttACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGc                                                                                         <  2:4352376/93‑1 (MQ=255)

  ttACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCCCTGTCGGGc                                                                                         >  1:4989935/1‑93 (MQ=255)

   tACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCa                                                                                        >  1:1883858/1‑93 (MQ=255)

    acccacTGCCATAAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCAt                                                                                       <  1:4942812/90‑1 (MQ=255)

    acacacTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCAt                                                                                       >  1:4653758/1‑93 (MQ=255)

     cacacTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATc                                                                                      >  1:384812/1‑93 (MQ=255)

      acacTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCa                                                                                     <  2:6935951/93‑1 (MQ=255)

       cacTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAg                                                                                    <  2:5555696/93‑1 (MQ=255)

       cacTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAg                                                                                    >  2:1819902/1‑93 (MQ=255)

        acTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAgg                                                                                   <  2:2140427/93‑1 (MQ=255)

        acTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAgg                                                                                   <  2:6857380/93‑1 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGc                                                                                  <  2:6100646/93‑1 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGc                                                                                  <  2:1680702/93‑1 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGc                                                                                  <  2:6746317/93‑1 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGc                                                                                  >  2:3235101/1‑93 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGc                                                                                  >  1:2577030/1‑93 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGc                                                                                  >  1:6540549/1‑93 (MQ=255)

         cTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCAAGCATCAGGc                                                                                  <  1:6035932/93‑1 (MQ=255)

          tGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCg                                                                                 >  1:6319580/1‑93 (MQ=255)

          tGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCg                                                                                 >  2:5509538/1‑93 (MQ=255)

          tGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCg                                                                                 <  1:5655537/93‑1 (MQ=255)

          tGCCATTAACGCCAGCGACATCGAATTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCg                                                                                 >  2:7381281/1‑93 (MQ=255)

           gCCATTACCGCCAGCGACATCGGTTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCGAGCCGCTGTCGGGCATCAGGCgg                                                                                <  1:6189397/93‑1 (MQ=255)

           gCCAATAACGCCAGAAAACTCCGGTTATGGCAGCCGGTTGTGATGTTGATTTAAATGGCCGTGACCTCCGAGCCGCTGTCGGGCATCAGGCgg                                                                                <  1:3332537/93‑1 (MQ=255)

               ttAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCaa                                                                            >  1:1095689/1‑93 (MQ=255)

               ttAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCaa                                                                            >  2:2396403/1‑93 (MQ=255)

               ttAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGGGACCACCGAGCCGCTGTCGGGCATCAGGCGGGGaa                                                                            >  1:5978658/1‑93 (MQ=255)

                tAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCGAGCCGCTGTCGGGCATCAGGCGGGCAAt                                                                           <  1:214823/93‑1 (MQ=255)

                tAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAAt                                                                           >  1:5519019/1‑93 (MQ=255)

                  aCGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCt                                                                         <  1:3058811/93‑1 (MQ=255)

                     ccAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGCCCCCCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGct                                                                      <  1:6409641/93‑1 (MQ=255)

                     ccAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGct                                                                      >  2:4047069/1‑93 (MQ=255)

                     ccAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGct                                                                      <  2:1501072/93‑1 (MQ=255)

                        gCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCtctc                                                                   <  2:5916801/93‑1 (MQ=255)

                        gCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCtctc                                                                   >  2:5415538/1‑93 (MQ=255)

                         cGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCg                                                                  <  2:6579598/93‑1 (MQ=255)

                          gACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGc                                                                 >  1:7019672/1‑93 (MQ=255)

                            cATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCaa                                                               >  2:6670436/1‑93 (MQ=255)

                              tCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATa                                                             >  2:3718482/1‑93 (MQ=255)

                              tCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATa                                                             >  2:4768198/1‑93 (MQ=255)

                               cGAGTTAGGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATAc                                                            <  1:4102300/93‑1 (MQ=255)

                                gAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACt                                                           <  2:6305338/93‑1 (MQ=255)

                                 aGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTc                                                          <  1:945828/93‑1 (MQ=255)

                                   ttATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGc                                                        <  2:2069053/93‑1 (MQ=255)

                                   ttATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGc                                                        <  1:5037347/93‑1 (MQ=255)

                                   ttATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGc                                                        >  2:5268267/1‑93 (MQ=255)

                                      tGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGCGCCGCTGGCGGGCATCAGGCGGACCATCTGCTCCCGCAATACTCGATGa                                                     >  2:294912/1‑93 (MQ=255)

                                      tGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGa                                                     >  2:783052/1‑93 (MQ=255)

                                      tGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGa                                                     >  2:3491164/1‑93 (MQ=255)

                                      tGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGa                                                     <  2:6162931/93‑1 (MQ=255)

                                         cAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCg                                                  <  2:337125/93‑1 (MQ=255)

                                            ccGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCt                                               >  2:5687020/1‑93 (MQ=255)

                                              gCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTaa                                             >  2:4467127/1‑93 (MQ=255)

                                                tGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACa                                           >  1:1735527/1‑93 (MQ=255)

                                                tGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACa                                           >  2:3060242/1‑93 (MQ=255)

                                                 ggTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACaa                                          >  2:5355070/1‑93 (MQ=255)

                                                     atgatgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCa                                      <  1:713415/93‑1 (MQ=255)

                                                     atgatgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCa                                      >  1:1205476/1‑93 (MQ=255)

                                                     ataatgatgTAATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCa                                      <  1:2416346/90‑1 (MQ=255)

                                                      tgatgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCAt                                     >  2:4638575/1‑93 (MQ=255)

                                                       gatgatgTCATCGGCGGTGCCCCCGGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATa                                    <  1:5735103/93‑1 (MQ=255)

                                                       gatgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATa                                    >  2:5149318/1‑93 (MQ=255)

                                                        atgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCCATACTCGCTGACCGGCTAACAAGTCATaa                                   >  2:2679760/1‑93 (MQ=255)

                                                        atgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATaa                                   <  1:5829189/93‑1 (MQ=255)

                                                         tgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAAc                                  <  2:921506/93‑1 (MQ=255)

                                                         tgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAAc                                  <  2:4967914/93‑1 (MQ=255)

                                                         tgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAAc                                  >  1:4281981/1‑93 (MQ=255)

                                                         tgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAAc                                  <  2:227223/93‑1 (MQ=255)

                                                         tgatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAAc                                  >  1:1448338/1‑93 (MQ=255)

                                                          gatgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAAcc                                 <  1:843521/93‑1 (MQ=255)

                                                            tgTCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCga                               <  2:3278108/93‑1 (MQ=255)

                                                              tCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCgaga                             <  2:5792072/93‑1 (MQ=255)

                                                              tCATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCgaga                             >  1:3012214/1‑93 (MQ=255)

                                                               cATCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGaa                            <  2:1636163/93‑1 (MQ=255)

                                                                aTCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAAc                           <  1:1645940/93‑1 (MQ=255)

                                                                 tCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACg                          <  1:7164091/93‑1 (MQ=255)

                                                                 tCGGCGGTGACCACCGAGCCGCTGTCGGGCATCAAGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACg                          >  1:4528847/1‑93 (MQ=255)

                                                                   gggggTGACCCCCGAGCCGCTGTGGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGt                        <  1:7491529/90‑1 (MQ=255)

                                                                    gcggTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGtt                       <  2:7132835/93‑1 (MQ=255)

                                                                    gcggTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGtt                       >  2:3808216/1‑93 (MQ=255)

                                                                      ggTGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTg                     >  1:5569214/1‑93 (MQ=255)

                                                                        tGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCa                   <  1:1013847/93‑1 (MQ=255)

                                                                        tGACCACCGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCa                   <  2:6410583/93‑1 (MQ=255)

                                                                           ccaccGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGAt                <  1:3583964/93‑1 (MQ=255)

                                                                             ccggAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTa              <  1:517564/90‑1 (MQ=255)

                                                                             accGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTa              >  2:6995920/1‑93 (MQ=255)

                                                                             accGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTa              <  2:194664/93‑1 (MQ=255)

                                                                              cggAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAt             <  1:6661919/91‑1 (MQ=255)

                                                                               cGAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtg            <  2:4212567/93‑1 (MQ=255)

                                                                                gAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgc           <  1:2044955/93‑1 (MQ=255)

                                                                                gAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgc           <  2:1236653/93‑1 (MQ=255)

                                                                                gAGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgc           >  1:2312548/1‑93 (MQ=255)

                                                                                gAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgc           >  1:133395/1‑93 (MQ=255)

                                                                                 aGCCGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgct          <  1:6134511/93‑1 (MQ=255)

                                                                                  gCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgctg         <  1:6505139/93‑1 (MQ=255)

                                                                                   ccGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgctgc        <  2:3869474/93‑1 (MQ=255)

                                                                                   ccGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTAtgctgc        >  1:5409819/1‑93 (MQ=255)

                                                                                     gCTGTCGGGCATCAGGCGGGCAATCTGCTCTCACAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGcac      >  1:50975/1‑93 (MQ=255)

                                                                                       tGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGcacca    <  2:2449041/93‑1 (MQ=255)

                                                                                         tCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCaccacc  >  1:3965109/1‑93 (MQ=255)

GTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACC  >  NC_000913/4568573‑4568754

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