BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,561,521	C→T	S754N (AGC→AAC)	b1489 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,561,521	0	C	T	99.5%	649.2 / NA	187	S754N (AGC→AAC)	b1489	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (87/99); total (87/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

TATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATT  >  NC_000913/1561431‑1561610

tatCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGttt                                                                                         <  2:7706271/93‑1 (MQ=255)

 atCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGtttt                                                                                        <  1:1189103/93‑1 (MQ=255)

 atCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGtttt                                                                                        <  2:3821204/93‑1 (MQ=255)

    ccTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGct                                                                                     >  1:6516473/1‑92 (MQ=255)

    ccTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGcc                                                                                     <  2:2979184/93‑1 (MQ=255)

     cTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGccc                                                                                    <  1:6947684/93‑1 (MQ=255)

     cTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGccc                                                                                    <  2:7477838/93‑1 (MQ=255)

     cTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGccc                                                                                    <  1:1001318/93‑1 (MQ=255)

       gAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCaa                                                                                  <  1:6106881/93‑1 (MQ=255)

       gAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCaa                                                                                  >  2:429229/1‑93 (MQ=255)

        aaTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAAt                                                                                 <  2:5792950/93‑1 (MQ=255)

        aaTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAAt                                                                                 <  2:5094663/93‑1 (MQ=255)

         aTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATg                                                                                <  2:3636780/93‑1 (MQ=255)

         aTAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATg                                                                                >  1:389963/1‑93 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                               >  2:4193226/1‑93 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                               >  1:4455087/1‑93 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                               <  1:327086/93‑1 (MQ=255)

          tAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGc                                                                               <  2:6745555/93‑1 (MQ=255)

           aaCGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCt                                                                              >  2:1509994/1‑93 (MQ=255)

             cgcgACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTgg                                                                            <  1:2372562/93‑1 (MQ=255)

              gcgACAGTGGGTCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAACTTGAGGTTTTGCCCTATGCTGGt                                                                           >  2:5386911/1‑93 (MQ=255)

              gcgACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGt                                                                           >  1:3695536/1‑93 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                         <  2:5634912/93‑1 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                         <  2:3133540/93‑1 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                         >  1:1996068/1‑93 (MQ=255)

                gACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTaa                                                                         >  1:7651470/1‑93 (MQ=255)

                 aCAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAAt                                                                        <  2:3916004/93‑1 (MQ=255)

                       gATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTc                                                                  >  2:2543117/1‑93 (MQ=255)

                       gATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTc                                                                  <  2:6818922/93‑1 (MQ=255)

                       gATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTc                                                                  >  2:4029045/1‑93 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATTGCTTCaa                                                                >  1:3609374/1‑93 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                                <  1:4281148/93‑1 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                                >  2:1712479/1‑93 (MQ=255)

                         tCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCaa                                                                >  1:1823659/1‑93 (MQ=255)

                            tGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTa                                                             <  1:1642901/93‑1 (MQ=255)

                             gCGTAGCATCTCAAATTTCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTaa                                                            <  2:2811031/93‑1 (MQ=255)

                             gCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTaa                                                            >  2:3158421/1‑93 (MQ=255)

                             gCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTaa                                                            <  1:7884760/93‑1 (MQ=255)

                              cGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAg                                                           <  2:6416176/93‑1 (MQ=255)

                               gTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAgg                                                          >  1:1318418/1‑93 (MQ=255)

                                tAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGc                                                         <  2:5196707/93‑1 (MQ=255)

                                 aGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCa                                                        >  1:7516643/1‑93 (MQ=255)

                                 aGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGACGAATGCTGGTAATGGCTTCAACAAAGGCa                                                        >  1:1799597/1‑93 (MQ=255)

                                  gCATCTCAAATCGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCaa                                                       <  2:3237790/93‑1 (MQ=255)

                                   cATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAg                                                      <  1:6628580/93‑1 (MQ=255)

                                   cATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAg                                                      <  2:1704425/93‑1 (MQ=255)

                                    aTCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAgg                                                     >  2:177514/1‑93 (MQ=255)

                                     tctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGa                                                    >  1:6154127/1‑93 (MQ=255)

                                      ctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGAt                                                   >  2:3433099/1‑93 (MQ=255)

                                      ctcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGAt                                                   <  2:4648017/93‑1 (MQ=255)

                                       tcAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATg                                                  >  2:2611237/1‑93 (MQ=255)

                                        cAAATTGCTCTTTGGTTTCGACGCCTGCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGc                                                 <  1:4167427/93‑1 (MQ=255)

                                         aaaTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCg                                                <  1:4867357/93‑1 (MQ=255)

                                          aaTTGTTCTTTGGTTTCGACCCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCAAATGGTGGAAATGGTTTCAAACAAGGCAGGGATGCGt                                               >  1:5471417/1‑93 (MQ=255)

                                           aTTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGtt                                              <  1:637437/93‑1 (MQ=255)

                                             tGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGtttt                                            <  1:4426783/93‑1 (MQ=255)

                                              gCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGttttt                                           >  2:4144279/1‑93 (MQ=255)

                                               ctctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATCCTGGTAATGGCTTCAAGTAAGGCACGGATGCGtttttg                                          >  1:2807594/1‑92 (MQ=255)

                                                tctTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCg                                         >  2:6180498/1‑93 (MQ=255)

                                                  tttGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGt                                       <  1:5288449/93‑1 (MQ=255)

                                                    tGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCa                                     <  2:1229838/93‑1 (MQ=255)

                                                    tGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCa                                     <  1:595218/93‑1 (MQ=255)

                                                     ggTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAg                                    >  1:1366125/1‑93 (MQ=255)

                                                      gTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGa                                   <  1:5809143/93‑1 (MQ=255)

                                                       tttCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGAc                                  >  2:1263485/1‑93 (MQ=255)

                                                        ttCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACa                                 <  1:7055884/93‑1 (MQ=255)

                                                         tCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACaa                                <  1:2916627/93‑1 (MQ=255)

                                                         tCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACaa                                <  1:7080879/93‑1 (MQ=255)

                                                         tCGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACaa                                <  1:2091112/93‑1 (MQ=255)

                                                          cGACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAAc                               >  1:5702048/1‑93 (MQ=255)

                                                           gACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACg                              >  2:4460016/1‑93 (MQ=255)

                                                           gACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACg                              >  1:6310502/1‑93 (MQ=255)

                                                           gACGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACg                              >  1:6705715/1‑93 (MQ=255)

                                                            aCGCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGa                             <  1:7440033/93‑1 (MQ=255)

                                                              gCCTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATc                           <  1:4198011/93‑1 (MQ=255)

                                                              cccTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATc                           <  1:2638813/92‑1 (MQ=255)

                                                                cTTCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGa                         <  1:1689691/93‑1 (MQ=255)

                                                                 ttCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                        >  2:6575910/1‑93 (MQ=255)

                                                                 ttCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                        >  1:5939429/1‑93 (MQ=255)

                                                                 ttCCGCCACGACGGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGAc                        >  2:4813919/1‑93 (MQ=255)

                                                                       cACGACGGTTAAATTGAGGCTGTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACt                  <  2:6785652/93‑1 (MQ=255)

                                                                        acgacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtt                 >  2:4669014/1‑93 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt               <  1:3952777/93‑1 (MQ=255)

                                                                          gacgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACtttt               <  1:3522129/93‑1 (MQ=255)

                                                                            cgGTTAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTgc             >  1:2392735/1‑92 (MQ=255)

                                                                               ttAAATTGAGGTTTTGCCCAATGCTGGCAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCat          >  1:2768885/1‑92 (MQ=255)

                                                                                tAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAt         <  1:3071650/93‑1 (MQ=255)

                                                                                tAAATTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAt         <  1:1380191/93‑1 (MQ=255)

                                                                                 aaaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtt        <  2:7879335/93‑1 (MQ=255)

                                                                                 aaaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtt        <  1:5421032/93‑1 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       >  1:682552/1‑93 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       <  2:240305/93‑1 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       <  1:4649443/93‑1 (MQ=255)

                                                                                  aaTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAACTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAttt       <  2:1921421/93‑1 (MQ=255)

                                                                                  aaTAGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAACCTTTTGTCAAttt       <  2:1002598/93‑1 (MQ=255)

                                                                                   aTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtttt      >  2:897086/1‑93 (MQ=255)

                                                                                   aTTGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAAtttt      >  2:3366580/1‑93 (MQ=255)

                                                                                     tGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGa    <  1:3020917/93‑1 (MQ=255)

                                                                                     tGAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGa    <  1:1186166/93‑1 (MQ=255)

                                                                                      gAGGTTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGAt   >  2:7338243/1‑93 (MQ=255)

TATCCCTGAATAACGCGACAGTGGATCTTGCGTAGCATCTCAAATTGCTCTTTGGTTTCGACGCCTTCCGCCACGACGGTTAAATTGAGGCTTTGCCCAATGCTGGTAATGGCTTCAAGTAAGGCAAGGATGCGTTTTTCGGTCAGACAACGATCGACAAAACTTTTGTCAATTTTGATT  >  NC_000913/1561431‑1561610

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