BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,228,096	A→G	N55S (AAC→AGC)	yjbC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,228,096	0	A	G	86.8%	741.3 / 99.9	302	N55S (AAC→AGC)	yjbC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (21/19); new base G (139/123); total (160/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

Only 100 of 303 total aligned reads displayed.

GGATCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATG  >  NC_000913/4228004‑4228185

  aTCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACg                                                                                         <  2:3943901/93‑1 (MQ=255)

   tCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCgg                                                                                        <  2:5674594/93‑1 (MQ=255)

   tCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCgg                                                                                        >  2:6269600/1‑93 (MQ=255)

   tCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCgg                                                                                        >  2:7054964/1‑93 (MQ=255)

    cGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGt                                                                                       >  1:4686580/1‑93 (MQ=255)

     gCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTc                                                                                      >  1:333930/1‑93 (MQ=255)

      cTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCa                                                                                     >  1:1753291/1‑93 (MQ=255)

      cTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCa                                                                                     <  2:2232427/93‑1 (MQ=255)

       tataTCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAg                                                                                    >  2:5218570/1‑93 (MQ=255)

       tataTCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAg                                                                                    >  2:3779840/1‑93 (MQ=255)

       tataTCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAg                                                                                    >  2:1900728/1‑93 (MQ=255)

        atatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGt                                                                                   >  2:4393891/1‑93 (MQ=255)

         tatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGtt                                                                                  >  2:7389720/1‑93 (MQ=255)

         tatCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGtt                                                                                  >  1:6915906/1‑93 (MQ=255)

          atCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttg                                                                                 >  1:3896063/1‑93 (MQ=255)

           tCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttga                                                                                <  2:5058029/93‑1 (MQ=255)

            cGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgat                                                                               <  2:1815/93‑1 (MQ=255)

            cGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgat                                                                               <  1:6102688/93‑1 (MQ=255)

             gAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgatt                                                                              >  2:230481/1‑93 (MQ=255)

             gAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgatt                                                                              <  2:2626753/93‑1 (MQ=255)

              aGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGttgattg                                                                             <  1:6219767/93‑1 (MQ=255)

                 aaGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAAc                                                                          >  1:6737963/1‑93 (MQ=255)

                  aGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAAcc                                                                         >  1:56997/1‑93 (MQ=255)

                   ggCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCt                                                                        >  2:2107134/1‑93 (MQ=255)

                     cAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCg                                                                      <  1:1242242/93‑1 (MQ=255)

                     cAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCg                                                                      >  1:1456048/1‑93 (MQ=255)

                     cAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCg                                                                      <  1:5229913/93‑1 (MQ=255)

                      aaTGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCgg                                                                     <  2:5515930/93‑1 (MQ=255)

                       aTGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCggg                                                                    >  2:6757074/1‑93 (MQ=255)

                       aTGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCggg                                                                    <  2:3379280/93‑1 (MQ=255)

                        tgtgTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGa                                                                   >  2:2942560/1‑93 (MQ=255)

                        tgtgTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGa                                                                   >  2:6826144/1‑93 (MQ=255)

                         gtgtTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGaa                                                                  <  2:2811506/93‑1 (MQ=255)

                          tgtTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAg                                                                 >  1:3805573/1‑93 (MQ=255)

                            ttCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGcc                                                               >  1:3097762/1‑93 (MQ=255)

                            ttCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGcc                                                               <  1:2286417/93‑1 (MQ=255)

                              ccTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGa                                                             >  1:3293381/1‑93 (MQ=255)

                               cTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGaa                                                            >  1:5049999/1‑93 (MQ=255)

                                ttAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAg                                                           <  2:4484045/93‑1 (MQ=255)

                                ttAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAg                                                           >  1:4573982/1‑93 (MQ=255)

                                ttAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAg                                                           >  2:1217068/1‑93 (MQ=255)

                                 tAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGa                                                          >  1:3674417/1‑93 (MQ=255)

                                 tAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGa                                                          <  2:4763482/93‑1 (MQ=255)

                                   aTGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGAtt                                                        >  1:4807926/1‑93 (MQ=255)

                                    tGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGAttt                                                       >  2:1092080/1‑93 (MQ=255)

                                       cAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGt                                                    >  1:6060156/1‑93 (MQ=255)

                                        aaGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTa                                                   <  1:2260030/93‑1 (MQ=255)

                                        aaGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTa                                                   >  1:7373867/1‑93 (MQ=255)

                                          gCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACt                                                 <  2:6105632/93‑1 (MQ=255)

                                            gAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACttg                                               <  1:4867499/93‑2 (MQ=255)

                                             aGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTAt                                              <  1:2060693/93‑1 (MQ=255)

                                             aGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTAt                                              <  2:2573672/93‑1 (MQ=255)

                                              gCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATc                                             >  2:1772248/1‑93 (MQ=255)

                                              gCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATc                                             >  1:876335/1‑93 (MQ=255)

                                              gCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATc                                             <  2:7184506/93‑1 (MQ=255)

                                               ccaccaTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCg                                            >  2:6529972/1‑93 (MQ=255)

                                                 accaTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGcc                                          >  1:1070692/1‑93 (MQ=255)

                                                   caTTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCt                                        >  2:1518244/1‑93 (MQ=255)

                                                         cGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACaa                                  <  1:3518612/93‑1 (MQ=255)

                                                         cGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACaa                                  <  2:6667177/93‑1 (MQ=255)

                                                          gATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAg                                 <  1:4301335/93‑1 (MQ=255)

                                                           aTCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGc                                >  2:148757/1‑93 (MQ=255)

                                                           aTCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGc                                >  1:5297473/1‑93 (MQ=255)

                                                           aTCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGc                                >  1:819486/1‑93 (MQ=255)

                                                           aTCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGc                                >  1:5825537/1‑93 (MQ=255)

                                                            tCAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGcc                               <  2:5332745/93‑1 (MQ=255)

                                                            tCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGcc                               >  1:2354757/1‑93 (MQ=255)

                                                             cAGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGccc                              <  1:832438/93‑1 (MQ=255)

                                                              aGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCg                             <  2:7487501/93‑1 (MQ=255)

                                                              aGGTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCg                             >  2:4952467/1‑93 (MQ=255)

                                                               ggTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGt                            <  1:397742/93‑1 (MQ=255)

                                                               ggTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGt                            >  2:2406655/1‑93 (MQ=255)

                                                                gTGAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGtt                           >  2:4031273/1‑93 (MQ=255)

                                                                  gAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTgg                         >  2:6034603/1‑93 (MQ=255)

                                                                  gAAACCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTgg                         <  1:5649375/93‑1 (MQ=255)

                                                                  gAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTgg                         >  1:2856970/1‑93 (MQ=255)

                                                                   aaaCCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGt                        >  2:6393460/1‑93 (MQ=255)

                                                                    aaCCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTa                       <  2:4303995/93‑1 (MQ=255)

                                                                     aCCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTAt                      >  1:1613642/1‑93 (MQ=255)

                                                                     aCCCGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTAt                      >  1:6397336/1‑93 (MQ=255)

                                                                      cccGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTAtt                     >  2:5923240/1‑93 (MQ=255)

                                                                       ccGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTg                    >  2:2495498/1‑93 (MQ=255)

                                                                       ccGGCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTg                    <  1:4338264/93‑1 (MQ=255)

                                                                         ggCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTa                  <  1:6004331/93‑1 (MQ=255)

                                                                         ggCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTa                  <  2:4162567/93‑1 (MQ=255)

                                                                          gCGACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTaa                 <  1:7069213/93‑1 (MQ=255)

                                                                            gACGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGc               >  2:4920666/1‑93 (MQ=255)

                                                                              cGTTGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcac             <  2:272488/93‑1 (MQ=255)

                                                                                ttGTGAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcacca           <  2:2059007/93‑1 (MQ=255)

                                                                                 tgtgAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac          <  2:1574564/93‑1 (MQ=255)

                                                                                 tgtgAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGAACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGcaccac          >  1:2971637/1‑93 (MQ=255)

                                                                                  gtgAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCaccacc         >  2:5865966/1‑93 (MQ=255)

                                                                                   tgAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCg        <  1:2946620/93‑1 (MQ=255)

                                                                                    gAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGa       <  1:4633581/93‑1 (MQ=255)

                                                                                    gAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGa       >  1:4141497/1‑93 (MQ=255)

                                                                                    gAAAGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGa       >  2:4698133/1‑93 (MQ=255)

                                                                                       aGTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGa    <  2:5401998/93‑1 (MQ=255)

                                                                                        gTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGAt   <  2:3559824/93‑1 (MQ=255)

                                                                                        gTAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGAt   <  1:2665180/93‑1 (MQ=255)

                                                                                         tAAGCGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATg  >  2:1423380/1‑93 (MQ=255)

GGATCGCTATATCGAGCAAGGCAATGTGTTCCTTAATGGCAAGCGAGCCACCATTGGCGATCAGGTGAAACCCGGCGACGTTGTGAAAGTAAACGGTCAGTTGATTGAACCTCGGGAAGCCGAAGATTTGGTACTTATCGCCCTGAACAAGCCCGTTGGTATTGTAAGCACCACCGAAGATG  >  NC_000913/4228004‑4228185

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