BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,650,355	T→C	F274L (TTC→CTC)	b1579 →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,650,355	0	T	C	99.6%	910.4 / NA	248	F274L (TTC→CTC)	b1579	putative transposase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (117/130); minor base G (0/1); total (117/131)
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 251 total aligned reads displayed.

GTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAT  >  NC_000913/1650265‑1650444

gTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTc                                                                                         <  1:523470/93‑1 (MQ=255)

gTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTc                                                                                         <  1:3911407/93‑1 (MQ=255)

  cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                                       >  2:1799570/1‑93 (MQ=255)

  cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                                       >  2:4786169/1‑93 (MQ=255)

  cattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCag                                                                                       <  1:5909978/93‑1 (MQ=255)

   attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                      >  2:726036/1‑93 (MQ=255)

   attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                      <  1:484786/93‑1 (MQ=255)

   attcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                      <  1:1615830/93‑1 (MQ=255)

   attcatATGATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATCACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCaga                                                                                      <  1:484784/93‑1 (MQ=255)

     tcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagaga                                                                                    >  1:545228/1‑93 (MQ=255)

      catCTTATTAAGCCAGCAATCGCCGCACTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGCGCATATCATTGATGTTCACCTCagagag                                                                                   <  1:5289660/93‑1 (MQ=255)

       atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                                  <  1:4159671/93‑1 (MQ=255)

       atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                                  <  1:4997380/93‑1 (MQ=255)

       atCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGt                                                                                  >  2:4494496/1‑93 (MQ=255)

        tCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTa                                                                                 <  2:139191/93‑1 (MQ=255)

           tattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATgg                                                                              >  1:967999/1‑93 (MQ=255)

            attaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                                             <  1:1081783/93‑1 (MQ=255)

              taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                           <  1:4542294/93‑1 (MQ=255)

              taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                           >  2:3737225/1‑93 (MQ=255)

                aGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaa                                                                         <  1:3565399/93‑1 (MQ=255)

                  ccAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaaca                                                                       <  2:4391403/93‑1 (MQ=255)

                   cAGCAATCGACGCTCTCCGGAGTCAGATAACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacag                                                                      <  2:3098069/93‑1 (MQ=255)

                    aGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                                     >  1:277643/1‑93 (MQ=255)

                       aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                  <  1:2465939/93‑1 (MQ=255)

                       aaTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                  >  2:6137279/1‑93 (MQ=255)

                        aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                                                 <  2:3219440/93‑1 (MQ=255)

                        aTCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaa                                                                 <  1:1822227/93‑1 (MQ=255)

                         tCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAc                                                                <  2:6120566/93‑1 (MQ=255)

                           gACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaa                                                              >  2:4909918/1‑93 (MQ=255)

                            aCGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaa                                                             >  2:1039685/1‑93 (MQ=255)

                             cGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaa                                                            <  2:910553/93‑1 (MQ=255)

                              gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGCAAAACaaaaa                                                           <  1:4071912/93‑1 (MQ=255)

                              gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaaa                                                           >  1:121951/1‑93 (MQ=255)

                              gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaaa                                                           <  2:5642342/93‑1 (MQ=255)

                               ctctCCGGGGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                          <  1:5129795/93‑1 (MQ=255)

                               ctctCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                          >  2:343375/1‑93 (MQ=255)

                                 ctcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGc                                                        <  2:5248019/93‑1 (MQ=255)

                                  tcCGGAGTCAGATGACATTAGCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                                       <  1:69024/93‑1 (MQ=255)

                                  tcCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCa                                                       <  1:6129855/93‑1 (MQ=255)

                                    cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCAcc                                                     >  1:5824485/1‑93 (MQ=255)

                                     ggAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCt                                                    <  1:4458166/93‑1 (MQ=255)

                                      gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATTCACCtt                                                   <  2:5015470/93‑1 (MQ=255)

                                      gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                   <  2:4525520/93‑1 (MQ=255)

                                      gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                   >  1:3012047/1‑93 (MQ=255)

                                      gAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCtt                                                   <  1:4792955/93‑1 (MQ=255)

                                       aGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCttt                                                  <  2:1916461/93‑1 (MQ=255)

                                       aGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCttt                                                  <  2:5284048/93‑1 (MQ=255)

                                        gTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTg                                                 <  2:1663887/93‑1 (MQ=255)

                                         tCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGt                                                <  2:1634043/93‑1 (MQ=255)

                                          cAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                               >  1:2811712/1‑93 (MQ=255)

                                          cAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtt                                               >  1:1025795/1‑93 (MQ=255)

                                           aGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                              >  1:142190/1‑93 (MQ=255)

                                            gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGtttt                                             >  1:962604/1‑93 (MQ=255)

                                             aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                            <  2:379733/93‑1 (MQ=255)

                                             aTGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttttt                                            <  2:1868387/93‑1 (MQ=255)

                                              tGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTc                                           <  1:3829041/93‑1 (MQ=255)

                                                aCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaa                                         >  1:1294121/1‑93 (MQ=255)

                                                  aTTAACGAGACTGAGCACAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCCCCTTTGTTTTTCAAcc                                       <  1:5485019/93‑1 (MQ=255)

                                                   ttAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                      >  1:5637943/1‑93 (MQ=255)

                                                    tAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                     <  1:1398758/93‑1 (MQ=255)

                                                       cGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAg                                  >  2:5140574/1‑93 (MQ=255)

                                                        gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgt                                 <  2:3631476/93‑1 (MQ=255)

                                                        gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgt                                 <  2:882310/93‑1 (MQ=255)

                                                         agaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtg                                <  1:1693720/93‑1 (MQ=255)

                                                          gaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAACAATGCACCTTTGTTTTTCAACCTGAAgtgt                               >  2:1215861/1‑93 (MQ=255)

                                                           aCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                              <  1:2621745/93‑1 (MQ=255)

                                                           aCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTc                              <  1:3183389/93‑1 (MQ=255)

                                                              gAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCTCCTTTGTTTTTCAACCTGAAGTGTCAGc                           <  2:2688819/93‑1 (MQ=255)

                                                                 cAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgag                        >  1:5402491/1‑93 (MQ=255)

                                                                  aaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga                       <  2:5571550/93‑1 (MQ=255)

                                                                   aaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGagag                      <  2:2165204/93‑1 (MQ=255)

                                                                    agagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGt                     <  1:213001/93‑1 (MQ=255)

                                                                      agCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                   >  2:4323232/1‑93 (MQ=255)

                                                                      agCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                   <  1:5348726/93‑1 (MQ=255)

                                                                      agCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                   >  1:897014/1‑93 (MQ=255)

                                                                      agCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                   >  2:3296001/1‑93 (MQ=255)

                                                                       gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                  >  2:2974971/1‑93 (MQ=255)

                                                                       gCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                  >  2:5451900/1‑93 (MQ=255)

                                                                        cATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                 <  2:5968627/93‑1 (MQ=255)

                                                                        cATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                 <  1:5798824/93‑1 (MQ=255)

                                                                        cATATCATTGACGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaa                 <  2:5968629/93‑1 (MQ=255)

                                                                          tatCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaaa               >  2:4713967/1‑93 (MQ=255)

                                                                            tCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAtt             >  1:2018592/1‑93 (MQ=255)

                                                                             cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTa            <  2:2022670/93‑1 (MQ=255)

                                                                              aTTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTAt           >  1:882803/1‑93 (MQ=255)

                                                                               ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATg          >  2:4015335/1‑93 (MQ=255)

                                                                               ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATg          <  2:260527/93‑1 (MQ=255)

                                                                                tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATgg         <  1:5419206/93‑1 (MQ=255)

                                                                                tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATgg         >  1:568526/1‑93 (MQ=255)

                                                                                 gATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGt        >  1:2255038/1‑93 (MQ=255)

                                                                                 gATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGt        >  2:1508042/1‑93 (MQ=255)

                                                                                  aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg       >  2:2485910/1‑93 (MQ=255)

                                                                                  aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg       <  1:3358146/93‑1 (MQ=255)

                                                                                   tGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGa      >  1:4329822/1‑93 (MQ=255)

                                                                                    gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc     <  1:3834402/93‑1 (MQ=255)

                                                                                    gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc     <  1:501632/93‑1 (MQ=255)

                                                                                    gTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc     >  1:4126950/1‑93 (MQ=255)

                                                                                     ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc    >  2:4106796/1‑93 (MQ=255)

                                                                                     ttCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAcc    >  1:3100881/1‑93 (MQ=255)

                                                                                      tCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCa   >  1:2962777/1‑93 (MQ=255)

GTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAT  >  NC_000913/1650265‑1650444

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