BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,626,669	G→A	A98A (GCG→GCA)	ydfH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,626,669	0	G	A	98.4%	661.1 / ‑2.7	190	A98A (GCG→GCA)	ydfH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/3); new base A (94/93); total (94/96)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 196 total aligned reads displayed.

TCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCGCAAGCA  >  NC_000913/1626579‑1626758

tCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGcacg                                                                                         <  2:1907039/93‑1 (MQ=255)

  cGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTc                                                                                       <  2:316750/93‑1 (MQ=255)

   gCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCg                                                                                      <  2:6473364/93‑1 (MQ=255)

    cAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCgg                                                                                     <  1:6194515/93‑1 (MQ=255)

     aaCGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCggg                                                                                    >  1:4750668/1‑93 (MQ=255)

      aCGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggc                                                                                   <  1:1045515/93‑1 (MQ=255)

      aCGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggc                                                                                   <  2:5855843/93‑1 (MQ=255)

         tGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggcggg                                                                                >  2:1786613/1‑92 (MQ=255)

         tGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggcggc                                                                                >  2:4121548/1‑93 (MQ=255)

         tGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggcggc                                                                                <  1:498938/93‑1 (MQ=255)

         tGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggcggc                                                                                >  2:2491873/1‑93 (MQ=255)

         tGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGggcggc                                                                                <  2:5522195/93‑1 (MQ=255)

             aGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGc                                                                            >  2:2425879/1‑93 (MQ=255)

             aGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGc                                                                            >  2:418385/1‑93 (MQ=255)

             aGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGc                                                                            >  2:4438585/1‑93 (MQ=255)

              gCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCAGCGAGCa                                                                           >  1:6403633/1‑93 (MQ=255)

                tACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATg                                                                         <  1:3737982/93‑1 (MQ=255)

                 aCGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGa                                                                        >  2:3962232/1‑93 (MQ=255)

                  cGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGAt                                                                       <  1:6948339/93‑1 (MQ=255)

                   gTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGAtt                                                                      >  1:7217813/1‑93 (MQ=255)

                   gTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGAtt                                                                      >  1:1995270/1‑93 (MQ=255)

                    tCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTa                                                                     >  2:4271391/1‑93 (MQ=255)

                      aacaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTAcc                                                                   >  1:4210809/1‑93 (MQ=255)

                      aacaaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTAcc                                                                   <  2:7711503/93‑1 (MQ=255)

                        caaAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGa                                                                 <  2:3523853/93‑1 (MQ=255)

                         aaaaTTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTAACGaa                                                                >  2:1432529/1‑93 (MQ=255)

                          aaaTTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGaaa                                                               <  2:7279000/93‑1 (MQ=255)

                          aaaTTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGaaa                                                               <  2:3287960/93‑1 (MQ=255)

                          aaaTTTCCATGGCCCAGGGGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGaaa                                                               <  2:5368844/93‑1 (MQ=255)

                              ttCCGTGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCa                                                           <  1:6535229/93‑1 (MQ=255)

                              ttCCATGGCCCAGGTGCGCAACGTCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCa                                                           <  2:5643121/93‑1 (MQ=255)

                               tCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAg                                                          >  1:7538900/1‑93 (MQ=255)

                               tCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAg                                                          <  2:1239618/93‑1 (MQ=255)

                                    ggCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCtt                                                     >  1:430448/1‑92 (MQ=255)

                                    ggCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTa                                                     <  1:781610/93‑1 (MQ=255)

                                     gCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTAt                                                    <  1:2058505/93‑1 (MQ=255)

                                     gCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTAt                                                    <  1:5034441/93‑1 (MQ=255)

                                      cccAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATc                                                   >  1:2669426/1‑93 (MQ=255)

                                      cccAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATc                                                   >  1:3335264/1‑93 (MQ=255)

                                       ccAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCa                                                  <  2:7591090/93‑1 (MQ=255)

                                        cAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCCTGATTACCGAAAGCCAGTGCTATCaa                                                 >  1:3094426/1‑93 (MQ=255)

                                         aGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAAc                                                <  2:4413382/93‑1 (MQ=255)

                                          ggTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACt                                               >  2:218912/1‑93 (MQ=255)

                                          ggTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACt                                               >  1:6369858/1‑93 (MQ=255)

                                          ggTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACt                                               <  1:5534734/93‑1 (MQ=255)

                                           gTGCGCAACGGCAGTTTTATCNGTCAGGCCATTGAGTGCGCGGTGGCACTTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTg                                              <  1:951530/93‑1 (MQ=255)

                                           gTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTg                                              <  1:6031409/93‑1 (MQ=255)

                                           gTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTg                                              <  1:4440953/93‑1 (MQ=255)

                                           gTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTg                                              <  2:7842902/93‑1 (MQ=255)

                                               gcAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAAc                                          <  1:2200875/93‑1 (MQ=255)

                                                cAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACa                                         >  1:5133194/1‑93 (MQ=255)

                                                 aaCGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaa                                        >  2:6768141/1‑93 (MQ=255)

                                                  aCGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaga                                       >  1:6090742/1‑91 (MQ=255)

                                                  aCGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaaa                                       >  2:3072876/1‑93 (MQ=255)

                                                  aCGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaaa                                       <  2:5737000/93‑1 (MQ=255)

                                                   cGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaaaa                                      <  1:6776189/93‑1 (MQ=255)

                                                    ggCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAAt                                     <  1:1192155/93‑1 (MQ=255)

                                                    ggCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAAt                                     >  1:6469552/1‑93 (MQ=255)

                                                     gCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATc                                    >  2:963339/1‑93 (MQ=255)

                                                      cAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCt                                   >  1:7426472/1‑93 (MQ=255)

                                                       aGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCtt                                  >  2:116931/1‑93 (MQ=255)

                                                       aGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCtt                                  <  2:5836194/93‑1 (MQ=255)

                                                        gTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTc                                 <  1:5261363/93‑1 (MQ=255)

                                                        gTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTc                                 <  2:1957652/93‑1 (MQ=255)

                                                         ttttATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCa                                <  1:3653386/93‑1 (MQ=255)

                                                          tttATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAc                               <  1:1476887/93‑1 (MQ=255)

                                                          tttATCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCAc                               <  1:1356473/93‑1 (MQ=255)

                                                              tCCGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGc                           >  2:6582723/1‑93 (MQ=255)

                                                                cGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCaa                         <  1:465537/93‑1 (MQ=255)

                                                                cGTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCaa                         <  1:1297116/93‑1 (MQ=255)

                                                                 gTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAAc                        >  1:1652030/1‑93 (MQ=255)

                                                                 gTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAAc                        >  1:6112796/1‑93 (MQ=255)

                                                                 gTCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAAc                        >  1:5175295/1‑93 (MQ=255)

                                                                 gTCAGGCCATTGAGAGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAAGCTTCACCAGCAAc                        >  2:5996021/1‑93 (MQ=255)

                                                                  tCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACg                       >  1:2108301/1‑93 (MQ=255)

                                                                  tCAGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACg                       >  2:7499778/1‑93 (MQ=255)

                                                                    aGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCa                     >  2:4650810/1‑93 (MQ=255)

                                                                    aGGCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCa                     <  1:6585783/93‑1 (MQ=255)

                                                                      gCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCAtt                   <  2:6966059/93‑1 (MQ=255)

                                                                      gCCATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCAtt                   <  2:5006423/93‑1 (MQ=255)

                                                                       gcATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTg                  <  1:1699166/92‑1 (MQ=255)

                                                                       ccATTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTg                  >  2:6827259/1‑93 (MQ=255)

                                                                         aTTGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGcc                >  2:2607550/1‑93 (MQ=255)

                                                                          ttGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAACTCTTCACCAGCAACGCATTGCCa               >  2:3935446/1‑93 (MQ=255)

                                                                          ttGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCa               >  1:7148428/1‑93 (MQ=255)

                                                                          ttGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCa               >  2:4500630/1‑93 (MQ=255)

                                                                           tGAGTGCGGGGGGGCACGTCGGGCGGCGAGCATTATTACCGAAAGCCAGTGCTATCAACTGAAACAAAATCTTCACCAGCAACGCATTGCCAt              <  2:2112733/93‑1 (MQ=255)

                                                                           tGAGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCAt              >  2:2894905/1‑93 (MQ=255)

                                                                             aGTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTg            >  1:1769147/1‑93 (MQ=255)

                                                                              gTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGa           >  2:5155570/1‑93 (MQ=255)

                                                                              gTGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGa           >  1:1577551/1‑93 (MQ=255)

                                                                               tGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACaaa                                       >  2:2021289/1‑64 (MQ=255)

                                                                               tGCGCGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAg          >  1:998612/1‑93 (MQ=255)

                                                                                gcgcGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAgc         <  1:5661869/93‑1 (MQ=255)

                                                                                gcgcGGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAgc         >  2:1743151/1‑93 (MQ=255)

                                                                                  gcgGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAgcgc       <  1:7644791/93‑1 (MQ=255)

                                                                                   cgGTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgca      <  1:1140564/93‑1 (MQ=255)

                                                                                    ggTGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaa     <  2:4727382/93‑1 (MQ=255)

                                                                                    ggTGCCACGTCGGGCGGCGAGCATGATTCCCGAAAGCCAGTGCTATCAACGGGAACAAAATCTTCCCCAGCAACGCATTGCCATTGAGCgcaa     <  1:1785538/93‑1 (MQ=255)

                                                                                      tGGCACGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCgcaagc   >  1:5068255/1‑93 (MQ=255)

TCCGCAACGTGGCAGCTACGTCAACAAAATTTCCATGGCCCAGGTGCGCAACGGCAGTTTTATCCGTCAGGCCATTGAGTGCGCGGTGGCGCGTCGGGCGGCGAGCATGATTACCGAAAGCCAGTGCTATCAACTGGAACAAAATCTTCACCAGCAACGCATTGCCATTGAGCGCAAGCA  >  NC_000913/1626579‑1626758

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