BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	147,141	A→G	G58G (GGA→GGG)	yadD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	147,141	0	A	G	99.3%	1035.5 / ‑2.9	301	G58G (GGA→GGG)	yadD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (158/141); total (159/142)
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 318 total aligned reads displayed.

GATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAA  >  NC_000913/147051‑147230

gaTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgca                                                                                         <  2:3440087/93‑1 (MQ=255)

gaTACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGgca                                                                                         <  1:7179375/93‑1 (MQ=255)

  tACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGcaca                                                                                       <  2:1277869/93‑1 (MQ=255)

    cATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAATGGGCACAGc                                                                                     <  1:2322123/93‑1 (MQ=255)

    cATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGc                                                                                     >  2:6951633/1‑93 (MQ=255)

     aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGTGCACAGCa                                                                                    >  1:1342429/1‑93 (MQ=255)

     aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCa                                                                                    <  2:3027957/93‑1 (MQ=255)

     aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCa                                                                                    <  2:3873236/93‑1 (MQ=255)

     aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCa                                                                                    <  2:2121212/93‑1 (MQ=255)

     aTTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCa                                                                                    <  2:7981783/93‑1 (MQ=255)

          ccAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGAc                                                                               <  1:296882/93‑1 (MQ=255)

          ccAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGAc                                                                               <  2:3611405/93‑1 (MQ=255)

           cAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACg                                                                              <  1:998177/93‑1 (MQ=255)

           cAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACg                                                                              >  1:6480983/1‑93 (MQ=255)

            aGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGt                                                                             <  1:4697493/93‑1 (MQ=255)

             gTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTg                                                                            >  2:1945272/1‑93 (MQ=255)

             gTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTg                                                                            <  2:3927944/93‑1 (MQ=255)

              tGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGc                                                                           >  2:8661430/1‑93 (MQ=255)

               ggAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGct                                                                          >  2:8723235/1‑93 (MQ=255)

               ggAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGct                                                                          >  2:7679079/1‑93 (MQ=255)

                  aTTACGCGAACTTTGTGGCCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTa                                                                       <  1:2888119/93‑1 (MQ=255)

                  aTTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCtata                                                                       >  2:6894960/1‑93 (MQ=255)

                  aTTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTa                                                                       <  1:3345156/93‑1 (MQ=255)

                   ttACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTAt                                                                      <  2:6324913/93‑1 (MQ=255)

                    tACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTAtt                                                                     >  2:234830/1‑93 (MQ=255)

                      cgcgAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTcc                                                                   >  2:204487/1‑93 (MQ=255)

                         gAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTg                                                                >  1:2702677/1‑93 (MQ=255)

                         gAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAAAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTg                                                                >  1:6176050/1‑93 (MQ=255)

                          aaCTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGc                                                               >  1:5638711/1‑93 (MQ=255)

                             tttGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCaaa                                                            <  2:2810372/93‑1 (MQ=255)

                                gtgACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGc                                                         <  1:6643133/93‑1 (MQ=255)

                                  gACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                                       >  1:1260315/1‑93 (MQ=255)

                                  gACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAg                                                       <  1:5110273/93‑1 (MQ=255)

                                   aCCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAgg                                                      <  1:8264315/93‑1 (MQ=255)

                                     cTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGc                                                    >  1:3584495/1‑93 (MQ=255)

                                     cTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGc                                                    <  2:6375168/93‑1 (MQ=255)

                                      tCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACATCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCa                                                   >  1:4395556/1‑93 (MQ=255)

                                      tCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCCCGGACGTGCTCTATTCTGCGCAAATGCAGGGCa                                                   >  1:8580413/1‑93 (MQ=255)

                                         acacGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATc                                                >  2:7546680/1‑93 (MQ=255)

                                          cacGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATcc                                               <  1:7386276/93‑1 (MQ=255)

                                           acGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATccc                                              <  1:7438715/93‑1 (MQ=255)

                                            cGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCg                                             <  1:4659402/93‑1 (MQ=255)

                                             gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATTCAGGGCAATCCCgg                                            >  1:686309/1‑93 (MQ=255)

                                             gCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCgg                                            >  1:8376240/1‑93 (MQ=255)

                                              cTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGt                                           <  2:6728723/93‑1 (MQ=255)

                                              cTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGt                                           <  1:3183855/93‑1 (MQ=255)

                                               ttCATTTAGAGTCGGGGAGTTTCATTGAAGAGAACCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGtt                                          >  2:2491673/1‑93 (MQ=255)

                                                tCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                                         <  1:873319/93‑1 (MQ=255)

                                                tCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGGGATCTATTCCGTGCAAATGCAGGGCAATCCCGGTTa                                         >  1:919809/1‑93 (MQ=255)

                                                 cATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTAt                                        >  2:5699177/1‑93 (MQ=255)

                                                  aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                       >  2:658676/1‑93 (MQ=255)

                                                  aTTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATc                                       >  1:748187/1‑93 (MQ=255)

                                                      agagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCa                                   >  1:498859/1‑93 (MQ=255)

                                                       gagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAt                                  <  1:1037742/93‑1 (MQ=255)

                                                       gagTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGCTATCTGCAt                                  >  2:4443779/1‑93 (MQ=255)

                                                          tCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgtt                               <  1:4108648/93‑1 (MQ=255)

                                                           cGGGGAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttg                              >  2:465194/1‑93 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                             >  2:5723832/1‑93 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                             >  2:4282761/1‑93 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                             >  1:8167722/1‑93 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgg                             >  2:6016311/1‑92 (MQ=255)

                                                            ggggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGCGCAAATGCAGGGCAATCCCGGTTATCTGCAtgttgt                             <  1:5344484/93‑1 (MQ=255)

                                                             gggAGTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTtgtg                            >  1:548637/1‑93 (MQ=255)

                                                                 gTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTTTGATTg                        <  2:6261258/93‑1 (MQ=255)

                                                                 gTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTg                        >  2:7858995/1‑93 (MQ=255)

                                                                 gTTTCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTg                        >  1:8516304/1‑93 (MQ=255)

                                                                  tttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGa                       <  1:4209892/93‑1 (MQ=255)

                                                                   ttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGaa                      <  1:3730363/93‑1 (MQ=255)

                                                                   ttCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCACGGTTATCTGCATGTTGTGATTGaa                      >  2:1689602/1‑93 (MQ=255)

                                                                    tCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAac                     >  2:8424174/1‑93 (MQ=255)

                                                                    tCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAac                     >  1:18684/1‑93 (MQ=255)

                                                                    tCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAac                     <  1:8084735/93‑1 (MQ=255)

                                                                    tCATTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAac                     <  2:227326/93‑1 (MQ=255)

                                                                      aTTGAAGAGAGCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAacac                   >  1:3445405/1‑93 (MQ=255)

                                                                       ttGAAGAGAGCCTGAAAGGGCACAGCAAGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCAGGTTATCTGCCTGTTGTGATTGAACAcc                  >  1:6129543/1‑93 (MQ=255)

                                                                         gAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCaa                <  1:5229442/93‑1 (MQ=255)

                                                                           agagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaag              <  1:5950845/93‑1 (MQ=255)

                                                                           agagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaag              <  2:7389998/93‑1 (MQ=255)

                                                                           agagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaag              <  2:2897268/93‑1 (MQ=255)

                                                                           agagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaag              >  2:4925252/1‑93 (MQ=255)

                                                                            gagagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagc             <  2:6820665/93‑1 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           <  1:301907/93‑1 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           <  2:545281/93‑1 (MQ=255)

                                                                              gagCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaa           <  2:8296928/93‑1 (MQ=255)

                                                                               agCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaag          <  2:3373832/93‑1 (MQ=255)

                                                                                gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc         <  2:2509692/93‑1 (MQ=255)

                                                                                gCCTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAaagcaagc         >  2:2655793/1‑93 (MQ=255)

                                                                                  cTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCg       <  2:4198867/93‑1 (MQ=255)

                                                                                  cTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCg       <  2:8219215/93‑1 (MQ=255)

                                                                                  cTGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCg       <  1:783574/93‑1 (MQ=255)

                                                                                   tGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCgg      >  1:7938411/1‑93 (MQ=255)

                                                                                   tGAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACAAAAGCAAGCCgg      <  2:112612/93‑1 (MQ=255)

                                                                                    gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa     <  1:1657826/93‑1 (MQ=255)

                                                                                    gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa     >  2:731642/1‑93 (MQ=255)

                                                                                    gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa     >  2:7119663/1‑93 (MQ=255)

                                                                                    gAAAGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGa     >  1:1483610/1‑93 (MQ=255)

                                                                                     aaaGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGAt    <  2:7026564/93‑1 (MQ=255)

                                                                                      aaGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATa   >  2:4064254/1‑93 (MQ=255)

                                                                                       aGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATaa  >  1:1595585/1‑93 (MQ=255)

                                                                                       aGGGCACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGGGATTGAACACCAAAGCAAGCCGGATag  >  1:4582863/1‑92 (MQ=255)

GATACATTTGCCAGTGGAATTACGCGAACTTTGTGACCTCAACACGCTTCATTTAGAGTCGGGGAGTTTCATTGAAGAGAGCCTGAAAGGACACAGCACGGACGTGCTCTATTCCGTGCAAATGCAGGGCAATCCCGGTTATCTGCATGTTGTGATTGAACACCAAAGCAAGCCGGATAA  >  NC_000913/147051‑147230

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