BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,374	C→T	D145N (GAC→AAC)	yiePs ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,374	0	C	T	99.7%	1154.6 / ‑5.6	308	D145N (GAC→AAC)	yiePs	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (133/173); total (134/174)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 308 total aligned reads displayed.

GCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGA  >  NC_000913/3938284‑3938463

gcgcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGc                                                                                         >  2:3548727/1‑93 (MQ=255)

gcgcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGc                                                                                         <  2:3930737/93‑1 (MQ=255)

gcgcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGc                                                                                         <  2:6496814/93‑1 (MQ=255)

  gcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTg                                                                                       <  2:1640306/93‑1 (MQ=255)

  gcCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTg                                                                                       >  1:4783895/1‑93 (MQ=255)

   ccTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCCTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTgg                                                                                      <  1:5478186/93‑1 (MQ=255)

    cTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGt                                                                                     >  2:2584536/1‑93 (MQ=255)

    cTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGt                                                                                     <  1:3767854/93‑1 (MQ=255)

    cTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGt                                                                                     >  1:2700075/1‑93 (MQ=255)

    cTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGt                                                                                     >  2:1280856/1‑93 (MQ=255)

     tGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTa                                                                                    >  2:944381/1‑93 (MQ=255)

        cAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAAtt                                                                                 <  2:5715202/93‑1 (MQ=255)

         aaGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTg                                                                                >  2:2825795/1‑93 (MQ=255)

          aGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGa                                                                               <  1:5144842/93‑1 (MQ=255)

          aGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGa                                                                               <  1:4066421/93‑1 (MQ=255)

          aGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGa                                                                               <  2:2235028/93‑1 (MQ=255)

           gCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGAc                                                                              >  2:3340619/1‑93 (MQ=255)

             tttAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGt                                                                            >  1:5601630/1‑93 (MQ=255)

                aaaTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaa                                                                         <  2:3399536/93‑1 (MQ=255)

                  aTGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagt                                                                       >  2:3789847/1‑93 (MQ=255)

                   tGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagta                                                                      >  1:3265429/1‑93 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                     <  2:654500/93‑1 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                     <  1:4706971/93‑1 (MQ=255)

                    gCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaa                                                                     <  2:1772641/93‑1 (MQ=255)

                     ccGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaag                                                                    <  2:6284779/93‑1 (MQ=255)

                     ccGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaag                                                                    <  1:2838925/93‑1 (MQ=255)

                     ccGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGaagtaag                                                                    <  2:4541415/93‑1 (MQ=255)

                       gcgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtg                                                                  <  1:6409333/93‑1 (MQ=255)

                        cgTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAgtgt                                                                 <  1:2940238/93‑1 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGa                                                               <  1:3129815/93‑1 (MQ=255)

                          tCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGa                                                               <  2:2084428/93‑1 (MQ=255)

                           cGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGAt                                                              <  1:6089173/93‑1 (MQ=255)

                            gCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATa                                                             >  2:5543880/1‑93 (MQ=255)

                             ccATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaa                                                            <  2:5814506/93‑1 (MQ=255)

                             ccATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaa                                                            >  2:6594259/1‑93 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                           >  1:4553492/1‑93 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                           >  2:2006535/1‑93 (MQ=255)

                              cATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATaaa                                                           >  2:3098192/1‑93 (MQ=255)

                               ctCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAc                                                          <  1:1028344/92‑1 (MQ=255)

                               aTCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAc                                                          <  2:517507/93‑1 (MQ=255)

                               aTCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAc                                                          >  1:5191914/1‑93 (MQ=255)

                                tCGCTTTGGATAATCGCATCGGCAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                         <  2:4397764/93‑1 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                         >  2:1426225/1‑93 (MQ=255)

                                tCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAAcc                                                         >  1:6219788/1‑93 (MQ=255)

                                 cGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCg                                                        <  2:1184690/93‑1 (MQ=255)

                                 cGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCg                                                        <  1:1375507/93‑1 (MQ=255)

                                  gCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGa                                                       <  1:1840074/93‑1 (MQ=255)

                                   cTTTGGACACGCTCATCGACAAGCGCCGGATGCAGCTCCAGCTTTATCACTGTGTGGCTGGTAATTGACGTGAAGCAAGTGTGACAAACCGaa                                                      <  1:2541260/93‑1 (MQ=255)

                                    tttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAAt                                                     <  2:1865732/93‑1 (MQ=255)

                                     ttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATg                                                    <  1:1608921/93‑1 (MQ=255)

                                     ttGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATg                                                    <  1:309413/93‑1 (MQ=255)

                                         gtAATCGCATCGACAATCGCTTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAAt                                                <  1:4912417/92‑1 (MQ=255)

                                         aTAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAAt                                                <  2:4765483/93‑1 (MQ=255)

                                          tAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATa                                               <  2:5069535/93‑1 (MQ=255)

                                           aaTCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAg                                              <  1:5018507/93‑1 (MQ=255)

                                           aaTCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAg                                              >  2:2382983/1‑93 (MQ=255)

                                            aTCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGc                                             >  2:4818108/1‑93 (MQ=255)

                                             tCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCg                                            >  2:4068297/1‑93 (MQ=255)

                                              cGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGa                                           <  1:2057000/93‑1 (MQ=255)

                                              cGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGAAAAACCGAATGGAATAGCGa                                           <  1:6529432/93‑1 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                          <  2:729206/93‑1 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                          >  2:6068279/1‑93 (MQ=255)

                                               gCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAg                                          >  1:4225568/1‑93 (MQ=255)

                                                 aTCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGc                                        <  1:4145228/93‑1 (MQ=255)

                                                  tCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCg                                       >  1:1917346/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                     >  1:5817479/1‑93 (MQ=255)

                                                    gACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGaa                                     <  1:3942475/93‑1 (MQ=255)

                                                     aCAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAAt                                    <  2:833934/93‑1 (MQ=255)

                                                      cAATCGCCTGATGCCGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGTAGTAAGTGTGATAAACCGAATAGAATAGCGAGGCGAAtt                                   >  1:2299571/1‑92 (MQ=255)

                                                       aaTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGa                                  <  2:2733021/93‑1 (MQ=255)

                                                       aaTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGa                                  <  1:5737267/93‑1 (MQ=255)

                                                        aTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAg                                 >  1:6062321/1‑93 (MQ=255)

                                                        aTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAg                                 <  2:5921360/93‑1 (MQ=255)

                                                        aTCGCCTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAg                                 <  1:497812/93‑1 (MQ=255)

                                                            ccTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCa                             >  1:2792065/1‑93 (MQ=255)

                                                             cTGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaa                            >  2:4243101/1‑93 (MQ=255)

                                                              tGATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaa                           >  1:4041573/1‑93 (MQ=255)

                                                               gATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaaa                          <  1:591416/93‑1 (MQ=255)

                                                               gATGCAGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCaaaa                          <  1:5782846/93‑1 (MQ=255)

                                                                    aGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGa                     >  1:889361/1‑93 (MQ=255)

                                                                    aGATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGa                     >  1:4246391/1‑93 (MQ=255)

                                                                     gATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAt                    <  1:1678804/93‑1 (MQ=255)

                                                                     gATCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAt                    >  2:2428965/1‑93 (MQ=255)

                                                                      aTCCAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGAtt                   <  2:500377/93‑1 (MQ=255)

                                                                         cAGCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgc                <  1:3405444/93‑1 (MQ=255)

                                                                           gCTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGgcgc              >  2:3508965/1‑93 (MQ=255)

                                                                            cTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCt             >  1:668896/1‑93 (MQ=255)

                                                                            cTTTATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCt             >  2:5035168/1‑93 (MQ=255)

                                                                             tttATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTc            >  2:3012860/1‑93 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt          <  1:344485/93‑1 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt          <  1:5541757/93‑1 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt          <  1:2217269/93‑1 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt          >  2:5407493/1‑93 (MQ=255)

                                                                               tATCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAt          <  1:3709237/93‑1 (MQ=255)

                                                                                aTCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAtt         <  2:5787439/93‑1 (MQ=255)

                                                                                 tCACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCAttt        >  1:363527/1‑93 (MQ=255)

                                                                                  cACTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAGTGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTc       <  1:983632/93‑1 (MQ=255)

                                                                                    cTGTGTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCAt     <  2:2297680/93‑1 (MQ=255)

                                                                                     tgtgTTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATa    <  1:2014975/93‑1 (MQ=255)

                                                                                      gtgtTGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAg   >  1:3467268/1‑93 (MQ=255)

GCGCCTGACAAGCTTTAAATGCCGCGTCGCCATCGCTTTGGATAATCGCATCGACAATCGCCTGATGCAGATCCAGCTTTATCACTGTGTCGCTGGTAATTGACGTGAAGTAAGTGTGATAAACCGAATGGAATAGCGAGGCGAATGAGGTCAAAAACGGATTGGCGCTCATTTCATAGA  >  NC_000913/3938284‑3938463

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