BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,186,123	C→T	R353C (CGC→TGC)	pepT →	putative peptidase T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,186,123	0	C	T	98.3%	630.4 / ‑2.5	179	R353C (CGC→TGC)	pepT	putative peptidase T
Reads supporting (aligned to +/- strand): ref base C (2/1); new base T (78/98); total (80/99)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 187 total aligned reads displayed.

CGCGAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAT  >  NC_000913/1186033‑1186212

cgcgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCtgc                                                                                         <  2:3768599/93‑1 (MQ=255)

 gcgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCtgcg                                                                                        <  1:4743920/93‑1 (MQ=255)

  cgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCgg                                                                                       <  2:3053966/93‑1 (MQ=255)

  cgAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCgg                                                                                       <  1:4352865/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      <  2:1764752/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      <  1:1376150/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      <  2:945238/93‑1 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      >  2:4396511/1‑93 (MQ=255)

   gagaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggt                                                                                      <  2:697875/93‑1 (MQ=255)

     gaAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtgg                                                                                    >  1:604693/1‑93 (MQ=255)

     gaAAGTGGTTCAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtgg                                                                                    <  1:3161531/93‑1 (MQ=255)

      aaaGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtggt                                                                                   <  2:1611638/93‑1 (MQ=255)

      aaaGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCggtggt                                                                                   <  1:2350514/93‑1 (MQ=255)

       aaGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTa                                                                                  >  1:1010589/1‑93 (MQ=255)

        aGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTAc                                                                                 <  1:141270/93‑1 (MQ=255)

            gTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGAc                                                                             >  2:4474113/1‑93 (MQ=255)

            gTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGAc                                                                             >  2:373118/1‑93 (MQ=255)

            gTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGAc                                                                             >  1:3828222/1‑93 (MQ=255)

            gTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGAc                                                                             <  1:2625519/93‑1 (MQ=255)

                 gCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGgcgc                                                                        >  2:433017/1‑93 (MQ=255)

                   aTCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGgcgcgc                                                                      >  2:4092932/1‑93 (MQ=255)

                    tCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCa                                                                     <  2:3453966/93‑1 (MQ=255)

                    tCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCa                                                                     >  2:1696686/1‑93 (MQ=255)

                     ccGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAg                                                                    <  1:3677432/93‑1 (MQ=255)

                       gCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGtt                                                                  >  2:1482654/1‑93 (MQ=255)

                       gCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGtt                                                                  >  2:4298166/1‑93 (MQ=255)

                         atatTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGt                                                                >  2:2621274/1‑93 (MQ=255)

                         atatTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGt                                                                >  1:4281071/1‑93 (MQ=255)

                          tatTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTc                                                               >  1:559821/1‑93 (MQ=255)

                          tatTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTc                                                               >  1:4027889/1‑93 (MQ=255)

                          tatTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTc                                                               <  1:1569062/93‑1 (MQ=255)

                            ttCTCGGTATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGt                                                             >  2:3624039/1‑93 (MQ=255)

                             tctcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGtt                                                            >  1:1039279/1‑93 (MQ=255)

                              ctcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGttt                                                           <  1:1959129/93‑1 (MQ=255)

                               tcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTa                                                          <  1:1763145/93‑1 (MQ=255)

                               tcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTa                                                          <  2:2470487/93‑1 (MQ=255)

                               tcGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTa                                                          >  2:2116967/1‑93 (MQ=255)

                               tcGATAACGCCCAGCAGGCGATGCTCGATTGCGATATTGAACCGGAACTTAAACCGATATGCGGGGGTACCGACGACGCGCAGTTGTCGTTTa                                                          <  2:4674035/93‑1 (MQ=255)

                                cGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTAt                                                         >  2:3377121/1‑93 (MQ=255)

                                 gATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATg                                                        <  1:2502601/93‑1 (MQ=255)

                                    atCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGa                                                     <  2:1139123/93‑1 (MQ=255)

                                    atCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGa                                                     >  2:1993091/1‑93 (MQ=255)

                                     tCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGAt                                                    >  1:3305849/1‑93 (MQ=255)

                                      cgcgCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGAtt                                                   <  1:2792381/93‑1 (MQ=255)

                                      cGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGAtt                                                   <  2:4175867/93‑1 (MQ=255)

                                      cGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGAtt                                                   <  1:2792380/93‑1 (MQ=255)

                                        cccAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTAc                                                 <  1:1560660/93‑1 (MQ=255)

                                         ccAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTAcc                                                >  2:1869955/1‑93 (MQ=255)

                                         ccAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTAcc                                                >  1:3619231/1‑93 (MQ=255)

                                           agcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGt                                              <  1:3542152/93‑1 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             <  1:749727/93‑1 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             >  1:3029082/1‑93 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             <  1:87984/93‑1 (MQ=255)

                                            gcagGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTg                                             >  1:3988767/1‑93 (MQ=255)

                                              agGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGcc                                           >  2:2920782/1‑93 (MQ=255)

                                              agGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGcc                                           <  2:1012598/93‑1 (MQ=255)

                                               ggCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGccc                                          <  2:1677043/93‑1 (MQ=255)

                                               ggCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGccc                                          >  2:789921/1‑93 (MQ=255)

                                                 cGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGa                                        <  1:4259639/93‑1 (MQ=255)

                                                  gATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGaa                                       >  2:2789396/1‑93 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAAcc                                     >  2:2204954/1‑93 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAAcc                                     <  1:1176493/93‑1 (MQ=255)

                                                    tGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAAcc                                     <  1:3345134/93‑1 (MQ=255)

                                                     gcgcGATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCt                                    <  1:4576494/93‑1 (MQ=255)

                                                        cgATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGtt                                 <  1:1013959/93‑1 (MQ=255)

                                                         gATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTc                                >  1:2606084/1‑93 (MQ=255)

                                                         gATTGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTc                                >  2:233738/1‑93 (MQ=255)

                                                           ttGCGATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGAGCCCGAACCTGTTCAc                              <  2:32119/93‑1 (MQ=255)

                                                               gATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGc                          >  2:902364/1‑93 (MQ=255)

                                                               gATATTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGc                          <  2:3917292/93‑1 (MQ=255)

                                                                 tatTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCgg                        >  1:2643421/1‑93 (MQ=255)

                                                                  atTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGt                       <  1:2437418/93‑1 (MQ=255)

                                                                  atTGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGt                       <  1:2590585/93‑1 (MQ=255)

                                                                   ttGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGtt                      <  1:4192948/93‑1 (MQ=255)

                                                                    tGAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTa                     >  2:4382766/1‑93 (MQ=255)

                                                                     gAACCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTAc                    >  1:4631668/1‑93 (MQ=255)

                                                                       aCCGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACaa                  >  1:495997/1‑93 (MQ=255)

                                                                        ccGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAAc                 >  1:2374092/1‑93 (MQ=255)

                                                                         cGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACt                <  1:3588105/93‑1 (MQ=255)

                                                                         cGGAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGACCGAACCTGTTCACTGGCGGCTACACCt                >  1:2055497/1‑93 (MQ=255)

                                                                          ggAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTa               >  2:3168452/1‑93 (MQ=255)

                                                                           gAACTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTAt              >  2:1545565/1‑93 (MQ=255)

                                                                            aaCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATc             <  1:3931989/93‑1 (MQ=255)

                                                                            aaCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATc             <  2:2031285/93‑1 (MQ=255)

                                                                            aaCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATc             <  1:712507/93‑1 (MQ=255)

                                                                            aaCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATc             <  2:1309679/93‑1 (MQ=255)

                                                                             aCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCa            <  1:1612677/93‑1 (MQ=255)

                                                                             aCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCa            >  1:236286/1‑93 (MQ=255)

                                                                             aCTGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCa            <  1:468119/93‑1 (MQ=255)

                                                                               tGAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATg          <  1:870365/93‑1 (MQ=255)

                                                                                gAAACCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATgg         <  2:3078580/93‑1 (MQ=255)

                                                                                 aaaCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGt        <  1:800188/93‑1 (MQ=255)

                                                                                 aaaCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGt        <  1:4153451/93‑1 (MQ=255)

                                                                                  agccGATCTGCGGTGGTAGCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTa       <  2:1675865/91‑1 (MQ=255)

                                                                                   aCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTaa      >  2:2190799/1‑93 (MQ=255)

                                                                                   aCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTaa      <  2:4279790/93‑1 (MQ=255)

                                                                                   aCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTaa      >  1:421942/1‑93 (MQ=255)

                                                                                   aCCGATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTaa      >  1:1812637/1‑93 (MQ=255)

                                                                                      gATCTGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCa   >  1:2170249/1‑93 (MQ=255)

                                                                                       aTCTGCGGTGGAACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAt  >  2:3424481/1‑93 (MQ=255)

CGCGAGAAAGTGGTTGAGCATCCGCATATTCTCGATATCGCCCAGCAGGCGATGCGCGATTGCGATATTGAACCGGAACTGAAACCGATCCGCGGTGGTACCGACGGCGCGCAGTTGTCGTTTATGGGATTACCGTGCCCGAACCTGTTCACTGGCGGTTACAACTATCATGGTAAGCAT  >  NC_000913/1186033‑1186212

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