BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,781,134	T→C	D32G (GAC→GGC)	b2648 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,781,134	0	T	C	93.8%	651.3 / 22.0	209	D32G (GAC→GGC)	b2648	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (5/8); new base C (119/77); total (124/85)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 218 total aligned reads displayed.

TTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGC  >  NC_000913/2781042‑2781223

ttACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGt                                                                                           <  2:4623407/93‑1 (MQ=255)

  aCATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGccc                                                                                         <  2:3408721/93‑1 (MQ=255)

   cATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGcccc                                                                                        >  1:3236098/1‑93 (MQ=255)

    aTAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCcccc                                                                                       >  1:3679679/1‑93 (MQ=255)

    aTAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCcccc                                                                                       <  1:2021106/93‑1 (MQ=255)

      aaTCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCcccccc                                                                                     >  1:4302900/1‑93 (MQ=255)

      aaTCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCcccccc                                                                                     >  2:5518918/1‑93 (MQ=255)

      aaTCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCcccccc                                                                                     >  1:5449976/1‑93 (MQ=255)

      aaTCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCcccccc                                                                                     >  2:896328/1‑93 (MQ=255)

       aTCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCa                                                                                    >  1:5996430/1‑93 (MQ=255)

        tCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCAt                                                                                   >  1:4208030/1‑93 (MQ=255)

           tttCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCAttgt                                                                                >  2:874000/1‑93 (MQ=255)

           tttCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCAttgt                                                                                >  2:3352734/1‑93 (MQ=255)

           tttCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCAttgt                                                                                >  2:3654614/1‑93 (MQ=255)

           tttCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCAttgt                                                                                >  2:4915974/1‑93 (MQ=255)

             tCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCAttgttg                                                                              <  1:4801131/93‑1 (MQ=255)

              cTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTgg                                                                             <  2:1934914/93‑1 (MQ=255)

               tgatgaATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTggg                                                                            >  2:2852270/1‑93 (MQ=255)

                   gaATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCaaa                                                                        >  2:4044782/1‑93 (MQ=255)

                   gaATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCaaa                                                                        >  2:5017166/1‑93 (MQ=255)

                   gaATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCaaa                                                                        >  2:3103260/1‑93 (MQ=255)

                    tatatCTTCTGCCCACATAAAAATCACACAATAACTTTGAGATCGTAGATTGTGTTACTTTTACAGCATTCGCCCTCCCATTGTTGGGCAAat                                                                       >  1:5494358/2‑93 (MQ=37)

                    aaTATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAat                                                                       >  2:1522413/1‑93 (MQ=255)

                      tatCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAatat                                                                     <  1:1049881/93‑1 (MQ=255)

                      tatCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAatat                                                                     >  1:2484959/1‑93 (MQ=255)

                      tatCTTCTGCTCACATAAAAATCACACAATAACATTGTGATCGCCGATTGTTTTACTTTTACAGCATTCGCCCCCCCCTTGTTGGGCAAataa                                                                     >  1:4451068/1‑92 (MQ=37)

                       atCTTCTGCTCTCATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAAtata                                                                    >  2:781388/1‑93 (MQ=255)

                       atCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAAtata                                                                    <  1:3322333/93‑1 (MQ=255)

                       atCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAAtata                                                                    <  2:845640/93‑1 (MQ=255)

                       atCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAAtata                                                                    >  2:5297277/1‑93 (MQ=255)

                         cttctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGa                                                                  <  2:2742428/93‑1 (MQ=255)

                         cttctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGa                                                                  <  2:1460262/93‑1 (MQ=255)

                         cttctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGa                                                                  <  2:5695374/93‑1 (MQ=255)

                         cttctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGa                                                                  <  2:3974328/93‑1 (MQ=255)

                          ttctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGAt                                                                 >  2:3468928/1‑93 (MQ=255)

                           tctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGAtt                                                                <  1:2708540/93‑1 (MQ=255)

                            ttGCTCACATAAAAATCACACAATAACTTTGAGCTCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTg                                                               <  1:4704759/92‑1 (MQ=255)

                            ctGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTg                                                               <  2:5432946/93‑1 (MQ=255)

                             tGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTgg                                                              >  2:2154672/1‑93 (MQ=255)

                             tGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTgg                                                              >  1:5774690/1‑93 (MQ=255)

                             tGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTgg                                                              >  1:3418248/1‑93 (MQ=255)

                              gCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTggg                                                             >  2:2646047/1‑93 (MQ=255)

                               cTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGc                                                            >  1:4724304/1‑93 (MQ=255)

                                tCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGcc                                                           >  2:1729130/1‑93 (MQ=255)

                                 cacaTAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCa                                                          >  2:793408/1‑93 (MQ=255)

                                  acaTAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCag                                                         >  2:1284838/1‑93 (MQ=255)

                                  acaTAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCag                                                         <  2:3260842/93‑1 (MQ=255)

                                   caTAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCaga                                                        <  2:4562097/93‑1 (MQ=255)

                                    aTAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCagag                                                       >  1:2923872/1‑93 (MQ=255)

                                    aTAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCagag                                                       >  1:3704705/1‑93 (MQ=255)

                                     tAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGc                                                      <  1:2239289/93‑1 (MQ=255)

                                     tAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGc                                                      >  2:5220014/1‑93 (MQ=255)

                                        aaaTCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACg                                                   >  1:5222998/1‑93 (MQ=255)

                                           tCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGaaa                                                <  1:4077519/93‑1 (MQ=255)

                                           tCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGaaa                                                <  2:2574601/93‑1 (MQ=255)

                                           tCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGaaa                                                <  1:121453/93‑1 (MQ=255)

                                           tCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGaaa                                                <  2:169860/93‑1 (MQ=255)

                                           tCACACAATAACTCTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGaaa                                                <  1:3348597/93‑1 (MQ=255)

                                            cacacaATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAg                                               <  2:6077577/93‑1 (MQ=255)

                                             acacaATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGt                                              <  1:2122660/93‑1 (MQ=255)

                                                caATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTaa                                           <  1:1245709/93‑1 (MQ=255)

                                                caATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTaa                                           >  2:2290094/1‑93 (MQ=255)

                                                caATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTaa                                           <  1:1358439/93‑1 (MQ=255)

                                                caATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTaa                                           >  1:776848/1‑93 (MQ=255)

                                                caATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTaa                                           >  2:5928620/1‑93 (MQ=255)

                                                caATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTaa                                           >  2:5558063/1‑93 (MQ=255)

                                                     aCTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCa                                      >  1:6082905/1‑93 (MQ=255)

                                                      cTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCAc                                     <  1:4025902/93‑1 (MQ=255)

                                                       tttGAGATCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACg                                    >  2:4267057/1‑93 (MQ=255)

                                                        ttGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGt                                   <  1:1694382/93‑1 (MQ=255)

                                                            gaTCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGc                               >  1:656010/1‑93 (MQ=255)

                                                            gaTCGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGc                               >  2:3732395/1‑93 (MQ=255)

                                                               cGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGcaca                            >  2:1163115/1‑93 (MQ=255)

                                                               cGCAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGcaca                            <  1:180211/93‑1 (MQ=255)

                                                                gCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCCCACAg                           >  2:4394479/1‑93 (MQ=255)

                                                                 cAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGc                          <  1:5082455/93‑1 (MQ=255)

                                                                 cAGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGc                          >  2:4526375/1‑93 (MQ=255)

                                                                 cAGATTGTTTTACTTTGACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGc                          <  1:12820/93‑1 (MQ=255)

                                                                  aGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGct                         >  2:1590650/1‑93 (MQ=255)

                                                                  aGATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGct                         >  1:5929848/1‑93 (MQ=255)

                                                                   gATTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGctc                        >  1:1985498/1‑93 (MQ=255)

                                                                    aTTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGctcc                       >  1:900922/1‑93 (MQ=255)

                                                                    aTTGTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGctcc                       <  1:3538180/93‑1 (MQ=255)

                                                                       gTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCtcctcc                    >  2:3859347/1‑93 (MQ=255)

                                                                       gTTTTACTTTTACAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCtcctcc                    >  1:253384/1‑93 (MQ=255)

                                                                        ttttacttttacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGATCGCACAGCTACTCCa                   >  1:4866007/1‑93 (MQ=255)

                                                                              ttttacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCa             >  1:5700585/1‑93 (MQ=255)

                                                                               tttacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCAc            >  2:3758815/1‑93 (MQ=255)

                                                                               tttacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCAc            >  2:1971623/1‑93 (MQ=255)

                                                                                 tacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGa          >  1:3156550/1‑93 (MQ=255)

                                                                                 tacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGa          >  1:5592691/1‑93 (MQ=255)

                                                                                 tacAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGa          <  2:1224982/93‑1 (MQ=255)

                                                                                  acAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGAc         <  1:926758/93‑1 (MQ=255)

                                                                                   cAGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACa        >  1:2662450/1‑93 (MQ=255)

                                                                                    aGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACaa       >  2:3197957/1‑93 (MQ=255)

                                                                                    aGCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACaa       >  1:435030/1‑93 (MQ=255)

                                                                                     gCATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAg      >  1:1466231/1‑93 (MQ=255)

                                                                                      cATTCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGa     <  1:988975/93‑1 (MQ=255)

                                                                                        ttCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATg   <  2:6103949/93‑1 (MQ=255)

                                                                                         tCGCCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGc  >  2:2955144/1‑93 (MQ=255)

TTACATAATCGTTTCTGATGAATATCTTCTGCTCACATAAAAATCACACAATAACTTTGAGATCGCAGATTGTTTTACTTTTACAGCATTCGTCCCCCCATTGTTGGGCAAATATAGATTGGGCCAGAGCACGAAAGTTAATACCACGTTCGCACAGCTCCTCCAACAGCACGACAAGATGC  >  NC_000913/2781042‑2781223

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