BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,340,190	A→G	T93A (ACT→GCT)	yrbI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,340,190	0	A	G	97.7%	667.9 / 2.5	214	T93A (ACT→GCT)	yrbI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/3); new base G (94/115); total (96/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 224 total aligned reads displayed.

GTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAT  >  NC_000913/3340101‑3340279

gTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctc                                                                                        >  2:3285225/1‑93 (MQ=255)

 ttGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctca                                                                                       <  1:3268644/93‑1 (MQ=255)

  tgtgCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGAtcgctcac                                                                                      <  2:230351/93‑1 (MQ=255)

       gcTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTa                                                                                 >  1:2571831/1‑93 (MQ=255)

       gcTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTa                                                                                 >  1:1095485/1‑93 (MQ=255)

        cTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTAt                                                                                <  2:1070571/93‑1 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               >  1:1281769/1‑93 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               >  2:4284413/1‑93 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               >  1:173935/1‑93 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               <  2:4596746/93‑1 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               <  1:2434828/93‑1 (MQ=255)

         tCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATc                                                                               >  2:1578406/1‑93 (MQ=255)

          cACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCa                                                                              <  1:1011841/93‑1 (MQ=255)

           aCCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAg                                                                             <  1:4313883/93‑1 (MQ=255)

           aCCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAg                                                                             >  2:4147668/1‑93 (MQ=255)

           aCCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAg                                                                             >  1:3233406/1‑93 (MQ=255)

            ccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgg                                                                            >  2:3842910/1‑93 (MQ=255)

            ccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgg                                                                            <  2:3103082/93‑1 (MQ=255)

            ccTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgg                                                                            >  1:1831335/1‑93 (MQ=255)

             ctctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAggg                                                                           >  1:974911/1‑93 (MQ=255)

             ctctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAggg                                                                           >  1:1045326/1‑93 (MQ=255)

              tctGATATTTAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  1:687530/1‑93 (MQ=255)

              tctGATATTGAAGTCGCTATCATTACCGGGCGAAATGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  1:633668/1‑93 (MQ=255)

              tctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  1:3884844/1‑93 (MQ=255)

              tctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAgggg                                                                          >  1:37352/1‑93 (MQ=255)

               ctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGc                                                                         <  2:3166323/93‑1 (MQ=255)

               ctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGc                                                                         <  1:522464/93‑1 (MQ=255)

               ctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGc                                                                         >  1:4106061/1‑93 (MQ=255)

               ctGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGc                                                                         >  1:4079347/1‑93 (MQ=255)

                tGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCa                                                                        <  2:3968694/93‑1 (MQ=255)

                tGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCa                                                                        >  2:4184509/1‑93 (MQ=255)

                tGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCa                                                                        >  2:929936/1‑93 (MQ=255)

                tGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCa                                                                        >  1:1303891/1‑93 (MQ=255)

                  atatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGt                                                                      >  2:1720344/1‑93 (MQ=255)

                   tatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTc                                                                     >  1:3895793/1‑93 (MQ=255)

                   tatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTc                                                                     <  2:211522/93‑1 (MQ=255)

                   tatTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTc                                                                     <  2:2721970/93‑1 (MQ=255)

                    atTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTca                                                                    >  2:1843219/1‑93 (MQ=255)

                     ttGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTcta                                                                   >  1:4694626/1‑91 (MQ=255)

                        aaGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaaca                                                                <  2:420849/93‑1 (MQ=255)

                         aGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTcaaacaa                                                               >  2:1968864/1‑93 (MQ=255)

                           tCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCaaacaaac                                                             <  2:3622421/93‑1 (MQ=255)

                           tCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCaaacaaac                                                             <  2:2444849/93‑1 (MQ=255)

                              cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          <  2:580911/93‑1 (MQ=255)

                              cTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGa                                                          >  2:4702492/1‑93 (MQ=255)

                               tATCATTACCGGGCGAAAGGCTATACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         <  2:3226610/93‑1 (MQ=255)

                               tATCATTACCGGGCGAAAGGCTAAACTGGGAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGAt                                                         >  1:2413614/1‑93 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATc                                                        <  2:4432235/93‑1 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATc                                                        <  2:1196468/93‑1 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATc                                                        >  1:2474361/1‑93 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATc                                                        <  2:4128846/93‑1 (MQ=255)

                                aTCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATc                                                        <  1:1127637/93‑1 (MQ=255)

                                 tCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCg                                                       >  1:2913631/1‑93 (MQ=255)

                                 tCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCg                                                       <  1:1997187/93‑1 (MQ=255)

                                  cATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGc                                                      >  1:2083064/1‑93 (MQ=255)

                                   aTTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGTGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGcc                                                     <  2:2521082/93‑1 (MQ=255)

                                    ttACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCt                                                    >  2:3741340/1‑93 (MQ=255)

                                     tgccGGTCGAAACTCTAAACTGGTAGAAGATCGTTTCGCCCGATTGGGGATCGCTCACGTGTATCTGGCGCAGTCAAGCAACCTGATCGCCtt                                                   <  2:1336011/91‑1 (MQ=255)

                                       ccGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTa                                                 <  2:86622/93‑1 (MQ=255)

                                       ccGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTa                                                 <  1:1777324/93‑1 (MQ=255)

                                       ccGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTa                                                 <  2:303638/93‑1 (MQ=255)

                                        cGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAg                                                <  1:2265318/93‑1 (MQ=255)

                                         gggCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGc                                               <  1:2661231/93‑1 (MQ=255)

                                         gggCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATAGCCTTTAGc                                               >  2:340446/1‑93 (MQ=255)

                                           gCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGTGGCAGTCAAACAAACTGATCGCCTTTAGCGa                                             >  2:3867320/1‑93 (MQ=255)

                                            cGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGAt                                            <  1:1073734/93‑1 (MQ=255)

                                              aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          <  2:2044482/93‑1 (MQ=255)

                                              aaaGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATct                                          >  2:3569421/1‑93 (MQ=255)

                                                 agCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgct                                       <  1:493005/92‑1 (MQ=255)

                                                  gCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATctgctg                                      >  2:1532087/1‑93 (MQ=255)

                                                   cTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTgg                                     >  2:265093/1‑93 (MQ=255)

                                                     aaaCTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAACCAAACTGATCGCCTTTAGCGATCTGCTTGaa                                   <  2:636607/93‑1 (MQ=255)

                                                        cTGGTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGaaaaa                                <  2:1821701/93‑1 (MQ=255)

                                                          tgTAGAAGATCGTTGTGCCACATTGGGGATCTCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACt                              <  1:3357259/92‑1 (MQ=255)

                                                          ggTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACt                              <  1:4197435/93‑1 (MQ=255)

                                                          ggTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACt                              <  2:1444573/93‑1 (MQ=255)

                                                           gTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTg                             <  1:1412107/93‑1 (MQ=255)

                                                           gTAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTg                             >  1:4467665/1‑93 (MQ=255)

                                                            tAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTgg                            <  1:4043934/93‑1 (MQ=255)

                                                            tAGAAGATCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTgc                            <  2:515532/93‑2 (MQ=255)

                                                              gaagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCg                          >  1:3694032/1‑93 (MQ=255)

                                                               aagaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGa                         <  1:2027928/93‑1 (MQ=255)

                                                                agaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAt                        >  2:196156/1‑93 (MQ=255)

                                                                 gaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       <  1:893758/93‑1 (MQ=255)

                                                                 gaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       <  1:4542173/93‑1 (MQ=255)

                                                                 gaTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       >  1:2753220/1‑93 (MQ=255)

                                                                 gaTCGTTGTGCCACATTGGGGATCGCTCACTTGCATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGAtt                       <  2:2814337/93‑1 (MQ=255)

                                                                  aTCGTTGTGCCTCATTGGGGATCGCTCGCTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTg                      <  1:2610687/93‑1 (MQ=255)

                                                                  aTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTg                      >  1:3235813/1‑93 (MQ=255)

                                                                  aTCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTg                      <  2:4598059/93‑1 (MQ=255)

                                                                   tCGTTGTGCCACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGc                     >  1:1339745/1‑93 (MQ=255)

                                                                           ccACATTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGaaa             >  2:672581/1‑93 (MQ=255)

                                                                              caTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtg          <  2:724458/93‑1 (MQ=255)

                                                                               aTTGGGGATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAAtgt         <  2:2050116/93‑1 (MQ=255)

                                                                                  ggggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGc      <  1:4242261/93‑1 (MQ=255)

                                                                                   gggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCt     <  2:3491542/93‑1 (MQ=255)

                                                                                   gggATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCt     >  1:3565589/1‑93 (MQ=255)

                                                                                    ggATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCtt    <  2:3916994/93‑1 (MQ=255)

                                                                                     gATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTa   <  1:1926014/93‑1 (MQ=255)

                                                                                     gATCGCTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTa   <  2:1853877/93‑1 (MQ=255)

GTTGTGCGCTCACCTCTGATATTGAAGTCGCTATCATTACCGGGCGAAAGGCTAAACTGGTAGAAGATCGTTGTGCCACATTGGGGATCACTCACTTGTATCAGGGGCAGTCAAACAAACTGATCGCCTTTAGCGATCTGCTGGAAAAACTGGCGATTGCCCCGGAAAATGTGGCTTAT  >  NC_000913/3340101‑3340279

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: `‑`