BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,223,632	A→G	L15P (CTG→CCG)	yohC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,223,632	0	A	G	97.7%	528.0 / ‑0.9	173	L15P (CTG→CCG)	yohC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/2); new base G (103/65); total (105/68)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

GCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTT  >  NC_000913/2223542‑2223721

gccATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGt                                                                                         >  1:2882538/1‑93 (MQ=255)

  cATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcc                                                                                       >  1:1948536/1‑93 (MQ=255)

   aTCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTccc                                                                                      >  1:4523907/1‑93 (MQ=255)

   aTCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTccc                                                                                      >  1:1320063/1‑93 (MQ=255)

   aTCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTccc                                                                                      >  2:1877069/1‑93 (MQ=255)

    tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcccc                                                                                     >  1:688894/1‑93 (MQ=255)

    tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcccc                                                                                     >  1:705108/1‑93 (MQ=255)

    tCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTcccc                                                                                     <  2:4548728/93‑1 (MQ=255)

      aGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAg                                                                                   >  1:98470/1‑93 (MQ=255)

       gCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGa                                                                                  <  1:1137121/93‑1 (MQ=255)

         aaaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACa                                                                                >  1:3883284/1‑93 (MQ=255)

         aaaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACa                                                                                >  2:3939548/1‑93 (MQ=255)

         aaaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATAGGGCTGGGAAACCGGTCCCCAGACa                                                                                >  1:643521/1‑93 (MQ=255)

          aaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACAt                                                                               <  2:1919655/93‑1 (MQ=255)

          aaaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACAt                                                                               >  1:3532984/1‑93 (MQ=255)

           aaCGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATg                                                                              >  2:4327817/1‑93 (MQ=255)

            aCGTGGTGGGTGGAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGACCCCAGACATgg                                                                             >  2:4374560/1‑93 (MQ=255)

              gtggtgGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCt                                                                           >  1:2438298/1‑93 (MQ=255)

              gtggtgGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCt                                                                           >  2:3142624/1‑93 (MQ=255)

               tggtggGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTc                                                                          >  1:2555328/1‑93 (MQ=255)

               tggtggGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTc                                                                          <  2:610597/93‑1 (MQ=255)

               tggtggGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTc                                                                          >  2:1444622/1‑93 (MQ=255)

                   gggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATAc                                                                      <  1:163350/93‑1 (MQ=255)

                    ggtgTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATAcc                                                                     >  1:1090881/1‑93 (MQ=255)

                        tAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccc                                                                 <  1:4096760/93‑1 (MQ=255)

                          aTGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccccc                                                               >  1:3169400/1‑93 (MQ=255)

                          aTGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccccc                                                               >  1:3189676/1‑93 (MQ=255)

                          aTGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccccc                                                               >  1:3676499/1‑93 (MQ=255)

                          aTGGTGAGAAATCGTTTCGTTTTAGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGccccc                                                               >  2:2566025/1‑93 (MQ=255)

                           tGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCt                                                              <  2:4171460/93‑1 (MQ=255)

                           tGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCt                                                              <  2:4682866/93‑1 (MQ=255)

                            ggTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAAAGGAGCCCAGACATGGCTCATACAGCCCCCtt                                                             >  1:733295/1‑93 (MQ=255)

                               gagaAATAGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCttgtt                                                          <  2:4490080/93‑1 (MQ=255)

                                agaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCttgttg                                                         <  1:1849987/93‑1 (MQ=255)

                                 gaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTtgttgt                                                        >  1:3624186/1‑93 (MQ=255)

                                 gaAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTtgttgt                                                        >  1:447830/1‑93 (MQ=255)

                                  aaaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTtgtg                                                       >  2:1217747/1‑93 (MQ=255)

                                   aaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTtgtgg                                                      >  2:1470743/1‑92 (MQ=255)

                                   aaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGa                                                      <  1:4349909/93‑1 (MQ=255)

                                   aaTCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGa                                                      <  2:488000/93‑1 (MQ=255)

                                     tcgtttcgttTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGt                                                    <  2:934877/93‑1 (MQ=255)

                                       gtttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGttt                                                  >  1:599285/1‑92 (MQ=255)

                                       gtttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGtta                                                  >  1:782859/1‑92 (MQ=255)

                                       gtttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTc                                                  <  1:935242/93‑1 (MQ=255)

                                         ttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCAt                                                <  2:884620/93‑1 (MQ=255)

                                         ttcgtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCAt                                                >  1:3790218/1‑93 (MQ=255)

                                            gtttTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGaa                                             <  2:2922053/93‑1 (MQ=255)

                                             ttttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAAc                                            >  1:1232935/1‑93 (MQ=255)

                                             ttttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAAc                                            <  2:423902/93‑1 (MQ=255)

                                             ttttCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCGGGCATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAAc                                            >  1:2687740/1‑93 (MQ=255)

                                                tCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTg                                         >  2:2346333/1‑93 (MQ=255)

                                                  gcgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGAc                                       >  2:852827/1‑93 (MQ=255)

                                                  gcgATTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGAc                                       >  1:2045232/1‑93 (MQ=255)

                                                     aTTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGt                                    >  1:1083464/1‑93 (MQ=255)

                                                     aTTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGt                                    >  2:163551/1‑93 (MQ=255)

                                                     aTTAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGt                                    >  2:1374269/1‑93 (MQ=255)

                                                       tAATCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTc                                  >  2:3582404/1‑93 (MQ=255)

                                                        aaTCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCa                                 <  2:4024047/93‑1 (MQ=255)

                                                          tCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAg                               >  1:4702087/1‑93 (MQ=255)

                                                          tCACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGCTGTGAGTTCATGAACCTGACAGTTCAAg                               >  1:1847714/1‑93 (MQ=255)

                                                           cACCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGt                              <  2:245455/93‑1 (MQ=255)

                                                           cACCTGCATTTCACCATCGGTGTGGGAACACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGt                              <  2:3384942/93‑1 (MQ=255)

                                                            cccTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGta                             <  2:3858203/92‑1 (MQ=255)

                                                            aCCTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAAta                             <  1:1768325/93‑1 (MQ=255)

                                                             ccTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtat                            >  2:2862863/1‑93 (MQ=255)

                                                              cTGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACGTGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGtata                           >  1:3052946/1‑93 (MQ=255)

                                                               tGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATaa                          <  2:4114931/93‑1 (MQ=255)

                                                               tGCATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATaa                          >  1:597190/1‑93 (MQ=255)

                                                                 cATTTCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGt                        <  2:1070045/93‑1 (MQ=255)

                                                                   tttCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCa                      >  2:2124591/1‑93 (MQ=255)

                                                                     tCACGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGc                    <  2:2395059/93‑1 (MQ=255)

                                                                        cGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtg                 <  2:670235/93‑1 (MQ=255)

                                                                        cGATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtg                 >  2:1859338/1‑93 (MQ=255)

                                                                         gATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                >  2:4281428/1‑93 (MQ=255)

                                                                         gATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                >  1:3225691/1‑93 (MQ=255)

                                                                         gATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                <  1:3375977/93‑1 (MQ=255)

                                                                         gATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                >  2:2056077/1‑93 (MQ=255)

                                                                         gATCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgt                <  2:3249752/93‑1 (MQ=255)

                                                                          aTCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTtgtg               >  1:793658/1‑93 (MQ=255)

                                                                           tCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTGAGCTtgtgt              >  2:3820897/1‑92 (MQ=255)

                                                                           tCGGGATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGa              <  1:2634938/93‑1 (MQ=255)

                                                                             gggATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGatt            >  2:2362191/1‑93 (MQ=255)

                                                                              ggATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattt           >  1:892450/1‑92 (MQ=255)

                                                                              ggATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGatta           >  2:2434755/1‑93 (MQ=255)

                                                                               gATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattat          <  2:3453669/93‑1 (MQ=255)

                                                                               gATGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattat          >  1:3599507/1‑93 (MQ=255)

                                                                                aTGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGattatt         <  1:1884662/93‑1 (MQ=255)

                                                                                 tGGGAAAACGGTCCCCAGACATGGCTCATAGCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGGTCAAGTATAAGGCAGCTTGTGATTAttt        >  1:2543219/1‑93 (MQ=255)

                                                                                 tGGGAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttt        <  1:350526/93‑1 (MQ=255)

                                                                                   ggAAAACGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttttt      >  1:235438/1‑93 (MQ=255)

                                                                                   ggAAAACGGTCCCCAGACATGACTCATACCGCCCCCTTGTTCTGAGTTCATGACCCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAttttt      <  1:531722/93‑1 (MQ=255)

                                                                                    gAAAACTGTCCCCAGACATGTCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTAtttttt     <  1:1499505/93‑1 (MQ=255)

                                                                                     aaaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg    <  2:3674235/93‑1 (MQ=255)

                                                                                     aaaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTg    <  1:244621/93‑1 (MQ=255)

                                                                                      aaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGt   <  2:215998/93‑1 (MQ=255)

                                                                                      aaaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGt   >  2:3423443/1‑93 (MQ=255)

                                                                                       aaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  <  2:264765/93‑1 (MQ=255)

                                                                                       aaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  <  1:4307849/93‑1 (MQ=255)

                                                                                       aaCGGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  <  1:120130/93‑1 (MQ=255)

                                                                                       aaCGGTCCCCAGACAGGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGtt  <  2:264764/93‑1 (MQ=255)

GCCATCAGCAAAACGTGGTGGGTGTAATGGTGAGAAATCGTTTCGTTTTCGCGATTAATCACCTGCATTTCACGATCGGGATGGGAAAACAGTCCCCAGACATGGCTCATACCGCCCCCTTGTTGTGAGTTCATGAACCTGACAGTTCAAGTATAAGTCAGCTTGTGATTATTTTTTGTT  >  NC_000913/2223542‑2223721

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