BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,312,334	A→G	L116P (CTG→CCG)	phnP ←	PhnP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,312,334	0	A	G	80.8%	458.5 / 118.7	234	L116P (CTG→CCG)	phnP	PhnP
Reads supporting (aligned to +/- strand): ref base A (16/29); new base G (102/87); total (118/116)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 237 total aligned reads displayed.

AGTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCCAGCGCAGCGGAAAC  >  NC_000913/4312244‑4312426

aGTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCggg                                                                                            >  2:6707022/1‑93 (MQ=255)

 gTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGc                                                                                           <  2:2437531/93‑1 (MQ=255)

 gTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGc                                                                                           >  2:4785523/1‑93 (MQ=255)

  ttttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGcc                                                                                          <  2:5800393/93‑1 (MQ=255)

  ttttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCGGAAATCAAGCGGGcc                                                                                          >  1:3088207/1‑93 (MQ=255)

   tttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGccc                                                                                         <  1:6713280/93‑1 (MQ=255)

   tttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGccc                                                                                         <  2:4224778/93‑1 (MQ=255)

   tttGAGTGGTTGAGCGGCAGGGGCCTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGccc                                                                                         >  2:1925227/1‑93 (MQ=255)

    ttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCg                                                                                        <  1:1703776/93‑1 (MQ=255)

    ttGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCg                                                                                        >  2:2082768/1‑93 (MQ=255)

       aGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGAt                                                                                     >  1:3565977/1‑93 (MQ=255)

       aGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGAt                                                                                     <  1:506698/93‑1 (MQ=255)

       aGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGAt                                                                                     >  2:4958789/1‑93 (MQ=255)

        gTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATg                                                                                    >  1:3381837/1‑93 (MQ=255)

         tGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGt                                                                                   <  2:30709/93‑1 (MQ=255)

         tGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGt                                                                                   >  2:6878179/1‑93 (MQ=255)

                gCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACa                                                                            <  1:1981937/93‑1 (MQ=255)

                gCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACa                                                                            <  2:5447535/93‑1 (MQ=255)

                gCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACa                                                                            <  1:6372977/93‑1 (MQ=255)

                  ggCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGa                                                                          <  2:5033915/93‑1 (MQ=255)

                  ggCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGa                                                                          >  2:6325596/1‑93 (MQ=255)

                  ggCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGa                                                                          <  2:1567669/93‑1 (MQ=255)

                   gCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAt                                                                         >  2:4915485/1‑93 (MQ=255)

                   gCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAt                                                                         <  1:2970788/93‑1 (MQ=255)

                   gCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAt                                                                         <  2:7210089/93‑1 (MQ=255)

                    cAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtc                                                                        <  2:45854/93‑1 (MQ=255)

                    cAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGAtc                                                                        <  1:4651676/93‑1 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgt                                                                      <  2:4050480/93‑1 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgt                                                                      >  1:1095922/1‑93 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgt                                                                      <  1:421088/93‑1 (MQ=255)

                      ggggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGAtcgt                                                                      <  2:389783/93‑1 (MQ=255)

                        ggCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGAtcgtcg                                                                    <  1:5785555/93‑1 (MQ=255)

                            tGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGc                                                                >  1:1188487/1‑93 (MQ=255)

                              aCCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGccc                                                              <  2:3470024/93‑1 (MQ=255)

                              aCCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGccc                                                              >  2:1751564/1‑93 (MQ=255)

                              aCCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGccc                                                              <  1:7388545/93‑1 (MQ=255)

                                 tGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGt                                                           <  1:3160926/93‑1 (MQ=255)

                                 tGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGt                                                           <  1:1469484/93‑1 (MQ=255)

                                  gTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTggt                                                          >  1:2429547/1‑91 (MQ=255)

                                  gTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGtt                                                          <  1:6258252/93‑1 (MQ=255)

                                   tAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTc                                                         >  1:352883/1‑93 (MQ=255)

                                   tAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTc                                                         <  2:4302896/93‑1 (MQ=255)

                                     aCCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCAt                                                       >  1:6634662/1‑93 (MQ=255)

                                     aCCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCAt                                                       >  1:5536014/1‑93 (MQ=255)

                                     aCCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCAt                                                       <  1:2532486/93‑1 (MQ=255)

                                      ccccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATc                                                      >  1:1639118/1‑93 (MQ=255)

                                       cccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATcc                                                     >  2:4134611/1‑93 (MQ=255)

                                       cccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATcc                                                     >  1:4054896/1‑93 (MQ=255)

                                        ccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcg                                                    >  2:6068176/1‑93 (MQ=255)

                                        ccTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcg                                                    >  2:3241547/1‑93 (MQ=255)

                                         cTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcgg                                                   >  2:3924278/1‑93 (MQ=255)

                                          tGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggc                                                  >  2:3582746/1‑93 (MQ=255)

                                          tGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggc                                                  >  2:2565585/1‑93 (MQ=255)

                                          tGCAAATCAAACACCACAAACGGTTCCACCGTGCGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCcggc                                                  >  2:822376/1‑93 (MQ=255)

                                             aaaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCggg                                               >  2:3829592/1‑93 (MQ=255)

                                             aaaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCggg                                               >  2:4729582/1‑93 (MQ=255)

                                              aaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGc                                              >  2:6492041/1‑93 (MQ=255)

                                              aaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGc                                              <  1:5388428/93‑1 (MQ=255)

                                              aaTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGc                                              >  2:7531134/1‑93 (MQ=255)

                                               aTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGcc                                             >  1:2938951/1‑93 (MQ=255)

                                               aTCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGcc                                             >  1:956253/1‑93 (MQ=255)

                                                   aaCACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTaa                                         >  1:6087964/1‑93 (MQ=255)

                                                     caccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAAc                                       <  1:825573/93‑1 (MQ=255)

                                                     caccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAAc                                       <  1:2807017/93‑1 (MQ=255)

                                                      accacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAAcc                                      >  1:3956767/1‑93 (MQ=255)

                                                       ccacAAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCg                                     <  1:5227280/93‑1 (MQ=255)

                                                        cacaAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCgg                                    <  2:3734413/93‑1 (MQ=255)

                                                          caAACGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGgga                                  <  2:2751970/93‑1 (MQ=255)

                                                           aaaCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggat                                 >  2:7375161/1‑93 (MQ=255)

                                                             aCGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcg                               <  2:329406/93‑1 (MQ=255)

                                                             aCGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcg                               >  1:1392182/1‑93 (MQ=255)

                                                             aCGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcg                               <  1:7305749/93‑1 (MQ=255)

                                                              cGGTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcgg                              <  1:1358833/93‑1 (MQ=255)

                                                                gTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggat                            <  2:2938858/93‑1 (MQ=255)

                                                                gTTCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggat                            >  1:6787571/1‑93 (MQ=255)

                                                                gTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggat                            <  1:2083971/93‑1 (MQ=255)

                                                                 ttCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGggatcggatc                           >  1:5817874/1‑93 (MQ=255)

                                                                  tCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGgatcggatcg                          >  2:5584155/1‑93 (MQ=255)

                                                                  tCCACCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGgatcggatcg                          >  1:409368/1‑93 (MQ=255)

                                                                    cACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGcc                        >  2:5064972/1‑93 (MQ=255)

                                                                     aCCGTGTGGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCa                       >  2:53087/1‑93 (MQ=255)

                                                                         tgtgGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGc                   >  1:4444734/1‑93 (MQ=255)

                                                                         tgtgGCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGc                   >  1:6577032/1‑93 (MQ=255)

                                                                          gtgGCTGAAATCAAGCGGGCCCGGATGTTTAAACGGATCGTCGCACCCCTGTTCATCCGGCGGGCCGCAAACCGGGACCGGATCGCCCACAcc                  >  1:5343524/1‑93 (MQ=255)

                                                                          gtgGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGcc                  >  2:2177424/1‑93 (MQ=255)

                                                                            ggCTGAAATCAAGCGGGCCCGGATGTTTAACCAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGcccc                <  2:1639094/93‑1 (MQ=255)

                                                                            ggCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGcccc                <  1:4721855/93‑1 (MQ=255)

                                                                            ggCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGcccc                <  1:1837101/93‑1 (MQ=255)

                                                                             gCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCca               >  2:3995412/1‑93 (MQ=255)

                                                                             gCTGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCca               >  1:7108888/1‑93 (MQ=255)

                                                                             gCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCCACGCCCca               >  1:3810687/1‑93 (MQ=255)

                                                                               tGAAATCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagc             >  1:6986594/1‑93 (MQ=255)

                                                                                 aaaTCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgc           <  2:7101168/93‑1 (MQ=255)

                                                                                  aaTCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgca          <  2:6587548/93‑1 (MQ=255)

                                                                                  aaTCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgca          >  1:1005123/1‑93 (MQ=255)

                                                                                  aaTCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTCCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgca          >  2:3178173/1‑93 (MQ=255)

                                                                                   aTCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgcag         >  2:6192270/1‑93 (MQ=255)

                                                                                   aTCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgcag         <  2:3034192/93‑1 (MQ=255)

                                                                                   aTCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgcag         <  2:1443977/93‑1 (MQ=255)

                                                                                    tCAAGCGGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCcagcgcagc        <  2:6998903/93‑1 (MQ=255)

AGTTTTGAGTGGTTGAGCGGCAGGGGCGTGACCTGTAACCCCTGCAAATCAAACACCACAAACGGTTCCACCGTGTGGCTGAAATCAAGCAGGCCCGGATGTTTAAACAGATCGTCGCAGCCCTGTTCATCCGGCGGGCCGTAAACCGGGATCGGATCGCCAACGCCCCAGCGCAGCGGAAAC  >  NC_000913/4312244‑4312426

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