BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,021,376	A→G	D580D (GAT→GAC)	b0960 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,021,376	0	A	G	97.6%	390.7 / ‑2.1	126	D580D (GAT→GAC)	b0960	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (57/66); total (60/66)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 131 total aligned reads displayed.

AAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCCTC  >  NC_000913/1021286‑1021465

aaaCGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTc                                                                                         <  1:1358234/93‑1 (MQ=255)

 aaCGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTcg                                                                                        >  2:1733276/1‑93 (MQ=255)

  aCGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTcgc                                                                                       <  1:512056/93‑1 (MQ=255)

   cgcgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGCGTGCGTcgcg                                                                                      >  1:1250368/1‑93 (MQ=255)

    gcgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCgcgg                                                                                     >  2:776330/1‑93 (MQ=255)

    gcgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCgcgg                                                                                     <  2:1180170/93‑1 (MQ=255)

    gcgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCgcgg                                                                                     <  1:2751875/93‑1 (MQ=255)

     cgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCgcggc                                                                                    >  2:3219268/1‑93 (MQ=255)

     cgGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCgcggc                                                                                    <  1:2461281/93‑1 (MQ=255)

      ggCTTCGCGAATTTGCGGGGGAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCgcggcg                                                                                   >  1:4140755/1‑93 (MQ=255)

       gCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGcggcgg                                                                                  <  1:1819219/93‑1 (MQ=255)

       gCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGcggcgg                                                                                  <  1:260232/93‑1 (MQ=255)

          tCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCa                                                                               >  2:2063908/1‑93 (MQ=255)

             cgAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAata                                                                            <  2:296145/93‑1 (MQ=255)

             cgAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAata                                                                            <  1:3545003/93‑1 (MQ=255)

                   tGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAatacgatac                                                                      <  1:1220419/93‑1 (MQ=255)

                      ggggTAACGTTCGGCTCGCTGGACATATTTGATGCCCCCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGcc                                                                   <  2:3383186/93‑1 (MQ=255)

                      ggggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGcc                                                                   <  2:239751/93‑1 (MQ=255)

                      ggggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGcc                                                                   >  2:4365791/1‑93 (MQ=255)

                       gggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCa                                                                  >  2:3767354/1‑93 (MQ=255)

                       gggTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCa                                                                  <  2:1251752/93‑1 (MQ=255)

                         gTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAAACAGATACGCCAGa                                                                >  1:1510018/1‑93 (MQ=255)

                          tAACGTTCGGCTCGCTTGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGAc                                                               <  2:2513768/93‑1 (MQ=255)

                           aaCGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACg                                                              <  1:1346394/93‑1 (MQ=255)

                            aCGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACgg                                                             <  2:524778/93‑1 (MQ=255)

                            aCGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACgg                                                             <  1:763962/93‑1 (MQ=255)

                             cGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGt                                                            <  1:2526314/93‑1 (MQ=255)

                             cGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGt                                                            >  1:1193093/1‑93 (MQ=255)

                              gTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAAGACGATTCGCAAGACGGtt                                                           >  2:2188400/1‑93 (MQ=255)

                              gTTCGGATCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGtt                                                           <  1:2562961/93‑1 (MQ=255)

                               ttCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTa                                                          <  2:3167366/93‑1 (MQ=255)

                               ttCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTa                                                          >  1:4550535/1‑93 (MQ=255)

                                tCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTAt                                                         >  2:1660785/1‑93 (MQ=255)

                                 cGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATc                                                        <  1:3758782/93‑1 (MQ=255)

                                  ggCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCa                                                       <  1:4468228/93‑1 (MQ=255)

                                  ggCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCa                                                       >  2:1879476/1‑93 (MQ=255)

                                  ggCCCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCa                                                       >  2:3548615/1‑93 (MQ=255)

                                   gCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCAc                                                      >  1:1850445/1‑93 (MQ=255)

                                     tCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGAGGGTTATCACGc                                                    >  1:576396/1‑93 (MQ=255)

                                     tCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGc                                                    >  2:4490259/1‑93 (MQ=255)

                                     tCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGc                                                    <  2:4157491/93‑1 (MQ=255)

                                      cGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGcc                                                   <  1:3407822/93‑1 (MQ=255)

                                         tGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCt                                                >  1:802236/1‑93 (MQ=255)

                                         tGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCt                                                <  1:1182506/93‑1 (MQ=255)

                                         tGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCt                                                <  1:4374737/93‑1 (MQ=255)

                                          ggACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTg                                               >  1:3266504/1‑93 (MQ=255)

                                           gACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGa                                              <  1:3354307/93‑1 (MQ=255)

                                            aCATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGAt                                             <  1:2087219/93‑1 (MQ=255)

                                            aCATATTTGATACCACCGACGCCAGCTAAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGAt                                             <  1:3997218/93‑1 (MQ=255)

                                               tatTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGt                                          >  1:3836587/1‑93 (MQ=255)

                                                atTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATggga                                         >  2:3563412/1‑91 (MQ=255)

                                                atTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTa                                         <  1:3526027/93‑1 (MQ=255)

                                                atTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTa                                         <  2:151622/93‑1 (MQ=255)

                                                 tttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAt                                        >  2:2628983/1‑93 (MQ=255)

                                                 tttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAt                                        <  2:667502/93‑1 (MQ=255)

                                                 tttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAt                                        <  1:3437527/93‑1 (MQ=255)

                                                  ttGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTAtt                                       <  2:2940330/93‑1 (MQ=255)

                                                    gATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGc                                     >  2:2535971/1‑93 (MQ=255)

                                                    gATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGc                                     >  2:2698295/1‑93 (MQ=255)

                                                    gATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGc                                     <  1:3780882/93‑1 (MQ=255)

                                                     aTACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGGATTGCt                                    >  2:3752505/1‑93 (MQ=255)

                                                       accaccGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTcc                                  <  2:4074257/93‑1 (MQ=255)

                                                         caccGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAg                                >  1:2636667/1‑93 (MQ=255)

                                                          accGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCGCGCCCCTGATGGTATTGCTCCAGt                               >  1:1925222/1‑93 (MQ=255)

                                                          accGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGt                               <  1:2143141/93‑1 (MQ=255)

                                                           ccGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTa                              <  1:2047359/93‑1 (MQ=255)

                                                           cagACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTTCTCCAGTa                              >  2:3909038/3‑93 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             >  1:3692859/1‑93 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             <  2:3885121/93‑1 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             <  2:2706312/93‑1 (MQ=255)

                                                            cGACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTAt                             <  1:3894009/93‑1 (MQ=255)

                                                             gACGCCAGCTCAGCTTCGCGGTTGTGTGCGTCGCGTCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATg                            <  1:1253788/93‑1 (MQ=255)

                                                             gACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATg                            <  1:1864164/93‑1 (MQ=255)

                                                             gACGCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATg                            <  1:1032863/93‑1 (MQ=255)

                                                                gCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCa                         >  1:1550765/1‑93 (MQ=255)

                                                                gCCAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCa                         >  2:1849299/1‑93 (MQ=255)

                                                                 ccAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCAt                        >  1:1637268/1‑93 (MQ=255)

                                                                 ccAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGGTATCACGCCCCTGATGGTATTGCTCCAGTATGGCac                        >  1:1378479/1‑92 (MQ=255)

                                                                  cAGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATc                       >  1:25051/1‑93 (MQ=255)

                                                                   aGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATATGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCg                      <  1:3032411/93‑1 (MQ=255)

                                                                   aGCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCg                      <  1:3032412/93‑1 (MQ=255)

                                                                   aGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCg                      >  1:4137163/1‑93 (MQ=255)

                                                                    gCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCga                     <  2:755035/93‑1 (MQ=255)

                                                                    gCTCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCga                     >  1:3198317/1‑93 (MQ=255)

                                                                      tCAGCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCgaga                   <  1:2759076/93‑1 (MQ=255)

                                                                         gCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAc                <  2:2097227/93‑1 (MQ=255)

                                                                         gCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAc                >  1:4035355/1‑93 (MQ=255)

                                                                         gCATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAc                >  1:3103091/1‑93 (MQ=255)

                                                                          cATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAcc               >  1:4281045/1‑93 (MQ=255)

                                                                          cATCACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATAcc               >  1:2776030/1‑93 (MQ=255)

                                                                             cACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGAc            <  1:1069563/93‑1 (MQ=255)

                                                                             cACGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGAc            >  2:4143981/1‑93 (MQ=255)

                                                                               cGGTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAg          >  2:3068598/1‑93 (MQ=255)

                                                                                ggTTGTGTGCGTCGCGGCGGGCAATCCGATACGCCAGACGGGTATCACGCCCCTGATGGTCTTGCGCCAGTCTGCAAGCGAGATACCGACAGt         >  1:3536046/1‑93 (MQ=255)

                                                                                ggTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGt         >  1:3536045/1‑93 (MQ=255)

                                                                                ggTTGTGTGCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGt         <  2:2789969/93‑1 (MQ=255)

                                                                                   tgtgtgCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTgg      <  2:559184/93‑1 (MQ=255)

                                                                                     tgtgCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGcc    <  1:579252/93‑1 (MQ=255)

                                                                                      gtgCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCct   >  1:283252/1‑93 (MQ=255)

                                                                                      gtgCGTCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCct   <  2:1977717/93‑1 (MQ=255)

AAACGCGGCTTCGCGAATTTGCGGGGTAACGTTCGGCTCGCTGGACATATTTGATACCACCGACGCCAGCTCAGCATCACGGTTGTGTGCATCGCGGCGGGCAATACGATACGCCAGACGGTTATCACGCCCCTGATGGTATTGCTCCAGTATGGCATCGAGATACCGACAGTTGGCCTC  >  NC_000913/1021286‑1021465

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