BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,848,059	T→C	V400A (GTT→GCT)	sppA →	protease IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,848,059	0	T	C	99.5%	643.7 / ‑4.5	187	V400A (GTT→GCT)	sppA	protease IV
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (106/80); total (106/81)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

Only 100 of 191 total aligned reads displayed.

TGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACC  >  NC_000913/1847969‑1848147

tGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTg                                                                                        <  2:3509574/93‑1 (MQ=255)

 gCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGt                                                                                       >  2:340903/1‑93 (MQ=255)

 gCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGt                                                                                       >  2:1305011/1‑93 (MQ=255)

 gCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGt                                                                                       >  1:998278/1‑93 (MQ=255)

    tgtTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATc                                                                                    <  1:3392391/93‑1 (MQ=255)

      ttAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGa                                                                                  <  1:3326359/93‑1 (MQ=255)

      ttAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGa                                                                                  <  2:3572856/93‑1 (MQ=255)

       tAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGAt                                                                                 <  2:6218689/93‑1 (MQ=255)

        aaTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATg                                                                                <  2:6394558/93‑1 (MQ=255)

        aaTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATg                                                                                >  1:2001157/1‑93 (MQ=255)

        aaTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATg                                                                                >  2:3661442/1‑93 (MQ=255)

         aTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATgg                                                                               >  2:5167866/1‑93 (MQ=255)

         aTAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATgg                                                                               >  2:4880308/1‑93 (MQ=255)

           aGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggc                                                                             <  2:2537478/93‑1 (MQ=255)

           aGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGAAGCCCGGGCAGCGGGCAAGCCTGTGGCTGTATCGATgggg                                                                             >  1:4240692/1‑92 (MQ=255)

           aGCCCAGGCGGAAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGACCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggc                                                                             >  1:2152688/1‑93 (MQ=255)

            gCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcg                                                                            >  1:1067923/1‑93 (MQ=255)

              ccAGGCGGCAGCGTTACCTCGTCTGAAGGGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggc                                                                          <  1:699139/93‑1 (MQ=255)

               cAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggca                                                                         <  2:5483884/93‑1 (MQ=255)

               cAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggca                                                                         >  1:3476337/1‑93 (MQ=255)

                aGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcat                                                                        >  1:322771/1‑93 (MQ=255)

                 ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatg                                                                       >  2:2003166/1‑93 (MQ=255)

                 ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatg                                                                       >  2:974069/1‑93 (MQ=255)

                 ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatg                                                                       >  2:3020453/1‑93 (MQ=255)

                 ggcggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatg                                                                       >  2:3304218/1‑93 (MQ=255)

                    ggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcg                                                                    <  1:1854855/93‑1 (MQ=255)

                    ggcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcg                                                                    <  1:552483/93‑1 (MQ=255)

                     gcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcgg                                                                   <  1:1732475/93‑1 (MQ=255)

                     gcAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcgg                                                                   >  2:960414/1‑93 (MQ=255)

                      cAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcggc                                                                  <  2:3191912/93‑1 (MQ=255)

                       aGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcggca                                                                 >  2:4306748/1‑93 (MQ=255)

                        gCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcggcct                                                                >  1:242410/1‑91 (MQ=255)

                        gCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGggcggcatggcggcat                                                                >  1:5464034/1‑93 (MQ=255)

                         cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATc                                                               <  1:1013941/93‑1 (MQ=255)

                         cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATc                                                               >  1:5266474/1‑93 (MQ=255)

                         cGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATc                                                               <  1:4042475/93‑1 (MQ=255)

                          gtgaCCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCt                                                              <  2:4383489/90‑1 (MQ=255)

                          gTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCt                                                              <  1:1004627/93‑1 (MQ=255)

                          gTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCt                                                              >  2:3675727/1‑93 (MQ=255)

                          gTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCt                                                              >  2:2428524/1‑93 (MQ=255)

                           gtACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtg                                                             <  1:5028085/92‑1 (MQ=255)

                             aCCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtggt                                                           >  2:5419778/1‑93 (MQ=255)

                              ccGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtggtg                                                          <  2:3759948/93‑1 (MQ=255)

                               cgcgTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCtggtgg                                                         <  2:149350/93‑1 (MQ=255)

                                   tCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTAc                                                     >  2:3569577/1‑93 (MQ=255)

                                    cTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACt                                                    >  1:2466793/1‑93 (MQ=255)

                                     tGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTg                                                   <  1:6459081/93‑1 (MQ=255)

                                     tGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTg                                                   <  1:4386938/93‑1 (MQ=255)

                                        aGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAt                                                >  2:1349202/1‑93 (MQ=255)

                                         gTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAtt                                               >  2:789128/1‑93 (MQ=255)

                                         gTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAtt                                               <  1:5447793/93‑1 (MQ=255)

                                         gTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGAtt                                               >  2:2644352/1‑93 (MQ=255)

                                              tCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCAc                                          <  1:3273165/93‑1 (MQ=255)

                                                 cgcTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGcc                                       >  1:5814799/1‑93 (MQ=255)

                                                     gAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTTCACGCCAGCt                                   >  1:4256189/1‑93 (MQ=255)

                                                     gAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCt                                   <  1:4288586/93‑1 (MQ=255)

                                                     gAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCt                                   >  1:2763706/1‑93 (MQ=255)

                                                      agcTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTa                                  <  2:5040341/91‑1 (MQ=255)

                                                      aaCTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTa                                  <  2:5919959/93‑1 (MQ=255)

                                                       aCTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTaa                                 >  2:1603096/1‑93 (MQ=255)

                                                         tGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAAtt                               <  1:1391214/93‑1 (MQ=255)

                                                             agcagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt                           >  1:4972203/1‑93 (MQ=255)

                                                             agcagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAt                           <  2:1863965/93‑1 (MQ=255)

                                                              gcagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAtt                          >  2:213885/1‑93 (MQ=255)

                                                              gcagcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACAtt                          >  2:3470615/1‑93 (MQ=255)

                                                                agcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTTCACGCCAGCTAATTACAttct                        >  1:1123019/1‑91 (MQ=255)

                                                                agcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                        <  1:1339184/93‑1 (MQ=255)

                                                                agcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgt                        >  1:4587133/1‑93 (MQ=255)

                                                                 gcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                       >  2:2082217/1‑93 (MQ=255)

                                                                 gcCCGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATtgtg                       >  2:235488/1‑93 (MQ=255)

                                                                  cccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTgg                      >  1:6030051/1‑93 (MQ=255)

                                                                  cccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTgg                      >  2:2868248/1‑93 (MQ=255)

                                                                   ccGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGc                     >  2:530803/1‑93 (MQ=255)

                                                                    cGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCt                    <  2:2901136/93‑1 (MQ=255)

                                                                    cGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCt                    >  1:1482543/1‑93 (MQ=255)

                                                                    cGGGCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCt                    >  2:4671512/1‑93 (MQ=255)

                                                                     gggCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTa                   <  1:2670687/93‑1 (MQ=255)

                                                                     gggCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTa                   >  2:2798051/1‑93 (MQ=255)

                                                                     gggCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTa                   <  1:4622124/93‑1 (MQ=255)

                                                                      ggCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTaa                  >  2:5407584/1‑93 (MQ=255)

                                                                      ggCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTaa                  <  2:4251902/93‑1 (MQ=255)

                                                                       gCAGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAc                 <  1:2196982/93‑1 (MQ=255)

                                                                         aGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAccc               >  1:3600598/1‑93 (MQ=255)

                                                                         aGCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAacc               >  1:1439761/1‑93 (MQ=255)

                                                                          gCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcccc              <  2:5586520/93‑1 (MQ=255)

                                                                          gCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcccc              >  1:5379575/1‑93 (MQ=255)

                                                                          gCGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAAcccc              >  1:3028755/1‑93 (MQ=255)

                                                                           cGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCa             <  1:4198844/93‑1 (MQ=255)

                                                                           cGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCa             <  2:2293847/93‑1 (MQ=255)

                                                                           cGGGTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCa             <  2:4242045/93‑1 (MQ=255)

                                                                            gggTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAg            >  1:5072424/1‑93 (MQ=255)

                                                                             ggTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGc           <  2:567311/93‑1 (MQ=255)

                                                                              gTAAGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCa          <  1:742914/93‑1 (MQ=255)

                                                                                aaGCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCAcc        >  2:3865625/1‑93 (MQ=255)

                                                                                  gCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCt      >  2:5092948/1‑93 (MQ=255)

                                                                                  gCCTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCt      >  2:6462177/1‑93 (MQ=255)

                                                                                  accTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTTGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCt      <  2:4546113/92‑1 (MQ=255)

                                                                                   ccTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTg     <  2:1397632/93‑1 (MQ=255)

                                                                                    cTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGa    >  2:1804662/1‑93 (MQ=255)

                                                                                    cTGTGGCTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGa    >  1:6182622/1‑93 (MQ=255)

TGCGTGTTAATAGCCCAGGCGGCAGCGTTACCGCGTCTGAAGTGATTCGCGCTGAACTGGCAGCAGCCCGGGCAGCGGGTAAGCCTGTGGTTGTATCGATGGGCGGCATGGCGGCATCTGGTGGTTACTGGATTTCCACGCCAGCTAATTACATTGTGGCTAACCCCAGCACCCTGACC  >  NC_000913/1847969‑1848147

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