BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,493,567	T→C	F86L (TTT→CTT)	b1423 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,493,567	0	T	C	99.5%	675.6 / ‑5.0	212	F86L (TTT→CTT)	b1423	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (98/113); total (98/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 225 total aligned reads displayed.

TAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGAT  >  NC_000913/1493475‑1493656

  aaggCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGctt                                                                                         <  1:4667148/90‑1 (MQ=255)

  aaTGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGctt                                                                                         >  1:9490000/1‑93 (MQ=255)

   aTGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTg                                                                                        >  2:3089967/1‑93 (MQ=255)

    tGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGc                                                                                       >  1:1163460/1‑93 (MQ=255)

    tGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGc                                                                                       >  2:2219783/1‑93 (MQ=255)

     gCAGACGCGCTGGCGCGACTGAATGTTGAACGTCATGGGGAGATTCGCGGTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGAGGCTTggc                                                                                      >  2:439120/1‑91 (MQ=255)

      cAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCa                                                                                     <  2:9497586/93‑1 (MQ=255)

      cAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCa                                                                                     <  2:7126958/93‑1 (MQ=255)

        gACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCAtt                                                                                   >  2:8148609/1‑93 (MQ=255)

        gACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCAtt                                                                                   >  2:2282093/1‑93 (MQ=255)

         aCGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTa                                                                                  >  2:3250231/1‑93 (MQ=255)

         aCGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTa                                                                                  >  2:171308/1‑93 (MQ=255)

           gAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATg                                                                                >  1:8572292/1‑93 (MQ=255)

             tcTGGCGCGACAGAATTTTGAACGTCGGGGCGCTATGCGCGTTGGGACTG‑‑CAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGa                                                                            <  2:5013531/92‑1 (MQ=255)

             gCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATggg                                                                              <  1:2589533/93‑1 (MQ=255)

              cTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATgggg                                                                             <  2:6827511/93‑1 (MQ=255)

                ggCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGAt                                                                           <  2:301297/93‑1 (MQ=255)

                 gcgcGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATg                                                                          <  1:3592655/93‑1 (MQ=255)

                   gcgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTa                                                                        <  2:9155018/93‑1 (MQ=255)

                    cgACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAc                                                                       <  1:3893398/93‑1 (MQ=255)

                     gACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAcc                                                                      <  1:9106395/93‑1 (MQ=255)

                     gACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAcc                                                                      >  1:9699068/1‑93 (MQ=255)

                     gACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAcc                                                                      >  1:9883604/1‑93 (MQ=255)

                      aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     >  2:3694147/1‑93 (MQ=255)

                      aCTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTAccc                                                                     <  2:6866187/93‑1 (MQ=255)

                       cTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCg                                                                    <  1:3184227/93‑1 (MQ=255)

                       cTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCg                                                                    <  2:641179/93‑1 (MQ=255)

                         gAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGt                                                                  <  2:8714178/93‑1 (MQ=255)

                         gAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGt                                                                  >  2:2606304/1‑93 (MQ=255)

                          aaTGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTg                                                                 <  2:9341951/93‑1 (MQ=255)

                             gTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGc                                                              >  1:8496979/1‑93 (MQ=255)

                              ttGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCt                                                             <  2:3301315/93‑1 (MQ=255)

                               tGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTa                                                            <  1:3664200/93‑1 (MQ=255)

                                   cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        >  1:5654039/1‑93 (MQ=255)

                                   cGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTAc                                                        <  1:3816721/93‑1 (MQ=255)

                                    gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACg                                                       >  1:237483/1‑93 (MQ=255)

                                    gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACg                                                       >  2:6210082/1‑93 (MQ=255)

                                    gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACg                                                       <  1:8266234/93‑1 (MQ=255)

                                    gTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGGGAGCTATTACg                                                       >  2:2078284/1‑93 (MQ=255)

                                     tCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGa                                                      >  1:5747869/1‑93 (MQ=255)

                                     tCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGa                                                      >  2:4180694/1‑93 (MQ=255)

                                      cATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAt                                                     <  2:2748311/93‑1 (MQ=255)

                                       aTGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtc                                                    <  1:6490332/93‑1 (MQ=255)

                                       aTGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtc                                                    <  2:8006504/93‑1 (MQ=255)

                                         ggggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctc                                                  >  1:2747876/1‑93 (MQ=255)

                                          gggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctct                                                 >  2:2147040/1‑93 (MQ=255)

                                          gggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctct                                                 >  2:8493898/1‑93 (MQ=255)

                                          gggCGATTCGCATTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGTGATGTACCCGGTGAGCTATTACGAtctct                                                 <  2:1950233/93‑1 (MQ=255)

                                           ggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctctc                                                >  2:9937119/1‑93 (MQ=255)

                                           ggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctctc                                                >  1:3298943/1‑93 (MQ=255)

                                           ggCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGAtctctc                                                >  1:528824/1‑93 (MQ=255)

                                            gCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCg                                               <  1:8517775/93‑1 (MQ=255)

                                              gATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgca                                             <  1:5188376/93‑1 (MQ=255)

                                              gATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCgca                                             <  1:4807428/93‑1 (MQ=255)

                                                ttCGCGTTGGGACTGCACAAGAGCTTGCTACTCTCTGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGTTATTACGATCTCTCgcagg                                           <  2:5403784/93‑1 (MQ=255)

                                                  cgcgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGCGAGGTACCCGGTGAGCTATCACCATCTCTCgcagcag                                         >  2:617458/1‑90 (MQ=255)

                                                    cgTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGcaggcagg                                       <  1:8382575/93‑1 (MQ=255)

                                                     gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      >  1:8204624/1‑93 (MQ=255)

                                                     gTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAggg                                      <  2:9407023/93‑1 (MQ=255)

                                                      ttGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAgggg                                     >  2:5922022/1‑93 (MQ=255)

                                                        gggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGGGCTATTACGATCTCTCGCAGGCAGGGGTg                                   >  2:8872259/1‑93 (MQ=255)

                                                         ggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGc                                  <  2:9186219/93‑1 (MQ=255)

                                                         ggACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGc                                  <  1:2067397/93‑1 (MQ=255)

                                                          gACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGcc                                 >  2:66009/1‑93 (MQ=255)

                                                          gACTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGcc                                 <  1:8126724/93‑1 (MQ=255)

                                                           aCTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGAGGTGCCg                                >  1:4471906/1‑93 (MQ=255)

                                                            cTGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCgg                               >  2:1380711/1‑93 (MQ=255)

                                                             tGCACAAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGGGCCgga                              >  1:3959222/1‑92 (MQ=255)

                                                               caccAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTAc                            <  2:1835147/93‑1 (MQ=255)

                                                                acaAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACa                           >  2:7665911/1‑93 (MQ=255)

                                                                 caAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAt                          >  1:1473578/1‑93 (MQ=255)

                                                                 caAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAt                          <  2:6123437/93‑1 (MQ=255)

                                                                 caAGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAt                          <  1:6422751/93‑1 (MQ=255)

                                                                  aaGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAtt                         <  1:9953622/93‑1 (MQ=255)

                                                                  aaGAGCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAtt                         >  1:5962296/1‑93 (MQ=255)

                                                                  aaGAGCTTGCTACTCTTCGGCGCATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACAtt                         <  1:2437554/93‑1 (MQ=255)

                                                                   agagCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTc                        <  1:3365083/93‑1 (MQ=255)

                                                                    gagCTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCg                       <  2:7875103/93‑1 (MQ=255)

                                                                       cTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACa                    <  1:1900464/93‑1 (MQ=255)

                                                                       cTTGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACa                    <  2:1439825/93‑1 (MQ=255)

                                                                        ttGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAg                   <  2:7746341/93‑1 (MQ=255)

                                                                         tGCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGc                  <  2:1771821/93‑1 (MQ=255)

                                                                          gCTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCa                 >  2:4072823/1‑93 (MQ=255)

                                                                           cTACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAt                <  1:7547386/93‑1 (MQ=255)

                                                                            tACTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAtt               <  1:6390189/93‑1 (MQ=255)

                                                                             aCTCTTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAttt              <  1:6075724/93‑1 (MQ=255)

                                                                             aCTCGTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCAttt              <  2:4987265/93‑1 (MQ=255)

                                                                              ctctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTc             <  1:5530146/93‑1 (MQ=255)

                                                                              ctctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGGGCCGGTACATTCGACAGCATTTc             >  2:6717885/1‑93 (MQ=255)

                                                                               tctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCg            <  2:3529392/93‑1 (MQ=255)

                                                                                ctTCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCgg           >  1:7173401/1‑93 (MQ=255)

                                                                                 ttCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGc          <  1:8741740/93‑1 (MQ=255)

                                                                                  tCGGCGGATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGcc         <  1:4702001/93‑1 (MQ=255)

                                                                                   cggcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATGTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGccc        >  1:9601011/1‑93 (MQ=255)

                                                                                    ggcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCa       >  1:4140622/1‑93 (MQ=255)

                                                                                     gcggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGAAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCAt      <  1:6314318/93‑1 (MQ=255)

                                                                                     gaggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCACAt      >  2:4803005/3‑93 (MQ=255)

                                                                                      cggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCAtt     <  1:4924568/93‑1 (MQ=255)

                                                                                      cggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCAtt     <  1:6647875/93‑1 (MQ=255)

                                                                                       ggATGCTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATtg    >  2:822033/1‑93 (MQ=255)

TAAATGCAGACGAGCTGGCGCGACTGAATGTTGAACGTCATGGGGCGATTCGCGTTGGGACTGCACAAGAGCTTGCTACTCTTCGGCGGATGTTTGCCATTATGGGGATGTACCCGGTGAGCTATTACGATCTCTCGCAGGCAGGGGTGCCGGTACATTCGACAGCATTTCGGCCCATTGAT  >  NC_000913/1493475‑1493656

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