BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,768,647	T→C	A12A (GCT→GCC)	b1689 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,768,647	0	T	C	98.3%	374.7 / ‑2.2	119	A12A (GCT→GCC)	b1689	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (65/52); total (67/52)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 121 total aligned reads displayed.

GATACTGTGACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCT  >  NC_000913/1768557‑1768735

gATACTGTGACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTa                                                                                        <  2:851162/93‑1 (MQ=255)

gATACTGTGACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTa                                                                                        <  2:4117673/93‑1 (MQ=255)

 aTACTGTGACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTAt                                                                                       <  1:2460562/93‑1 (MQ=255)

  tACTGTGACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATg                                                                                      >  1:3459470/1‑93 (MQ=255)

      gtgACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACt                                                                                  <  1:4187133/93‑1 (MQ=255)

        gACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTcc                                                                                <  2:2176101/93‑1 (MQ=255)

         gctctctTCAAGCCGTAAGCAGGCTGATAAAGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACTATGAACGCCTATGAACTCCa                                                                               <  2:1889373/92‑1 (MQ=255)

         aCTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCa                                                                               <  2:2560422/93‑1 (MQ=255)

         aCTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCa                                                                               >  1:1599105/1‑93 (MQ=255)

            ctctTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGc                                                                            >  2:1592474/1‑93 (MQ=255)

            ctctTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGc                                                                            >  1:1987078/1‑93 (MQ=255)

               ttCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACt                                                                         <  2:1564118/93‑1 (MQ=255)

                tCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTc                                                                        >  1:4610050/1‑93 (MQ=255)

                  aaGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCg                                                                      >  1:1203409/1‑93 (MQ=255)

                    gCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGcc                                                                    >  2:345004/1‑93 (MQ=255)

                     ccGTCAGCAGGCTGATAATGTTATTTATATTTTATAGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCa                                                                   >  2:312286/1‑93 (MQ=255)

                       gTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCata                                                                 <  2:4636740/93‑1 (MQ=255)

                       gTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCata                                                                 >  2:1648677/1‑93 (MQ=255)

                        tCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCatat                                                                >  1:2174051/1‑93 (MQ=255)

                             agGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTg                                                           >  2:3828628/1‑93 (MQ=255)

                              ggCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGc                                                          >  2:189773/1‑93 (MQ=255)

                                cTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCa                                                        >  1:4303055/1‑93 (MQ=255)

                                 tGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCAt                                                       <  1:2214301/93‑1 (MQ=255)

                                  gttAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATg                                                      >  2:1154724/3‑93 (MQ=255)

                                    tAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGAc                                                    <  2:1874341/93‑1 (MQ=255)

                                      aTGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTa                                                  <  1:1649078/93‑1 (MQ=255)

                                       tGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAt                                                 <  2:1326086/93‑1 (MQ=255)

                                        gTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                                <  1:4263889/93‑1 (MQ=255)

                                        gTTATTTATATTTTATTGTCTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTAtt                                                <  1:4616462/93‑1 (MQ=255)

                                         ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                               <  1:4613111/93‑1 (MQ=255)

                                         ttatttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtg                                               >  2:4744285/1‑93 (MQ=255)

                                          tatttatATTTTATTGTTTATTTCATGCAGCCATTACGCGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtga                                              <  2:3375022/93‑1 (MQ=255)

                                           atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                             <  1:1192602/93‑1 (MQ=255)

                                           atttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgat                                             <  1:3496912/93‑1 (MQ=255)

                                             ttatATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATtgatga                                           <  1:1811586/93‑1 (MQ=255)

                                               atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                         >  1:3746975/1‑93 (MQ=255)

                                               atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                         <  2:1128163/93‑1 (MQ=255)

                                               atatTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAAt                                         <  2:2961015/93‑1 (MQ=255)

                                                 atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGACCTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                       <  2:2562602/93‑1 (MQ=255)

                                                 atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCCCTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                       <  2:1146062/93‑1 (MQ=255)

                                                 atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCTTTACTATTGATGAATgc                                       >  2:2301477/1‑93 (MQ=255)

                                                 atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGCCTATTGATGAATgc                                       <  2:2500723/93‑1 (MQ=255)

                                                 atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                       >  1:2584902/1‑93 (MQ=255)

                                                 atTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgc                                       <  2:3403421/93‑1 (MQ=255)

                                                  ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                      >  1:2804362/1‑93 (MQ=255)

                                                  ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                      >  2:2529349/1‑93 (MQ=255)

                                                  ttttATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcg                                      >  2:2462332/1‑93 (MQ=255)

                                                   tttATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTAAGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATgcgc                                     <  2:3760506/93‑1 (MQ=255)

                                                    ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                    >  1:3238127/1‑93 (MQ=255)

                                                    ttATTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGcgcg                                    <  2:4574130/93‑1 (MQ=255)

                                                      aTTGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCTATGACTATTGATGAATGCGCGAc                                  >  1:3708197/1‑93 (MQ=255)

                                                       ttGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                                 >  2:489114/1‑93 (MQ=255)

                                                       ttGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                                 <  1:3885169/93‑1 (MQ=255)

                                                       ttGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGAcc                                 <  2:633638/93‑1 (MQ=255)

                                                        tGTTTATTTCATGCAGCCATTACACGGTGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                >  2:613033/1‑93 (MQ=255)

                                                        tGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATCTTTTTGCCATGACTATTGATGAATGGGCGACCt                                >  2:4039666/1‑93 (MQ=255)

                                                        tGTTTATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCt                                >  2:2804549/1‑93 (MQ=255)

                                                            tATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                            <  2:3843953/93‑1 (MQ=255)

                                                            tATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                            <  1:2870655/93‑1 (MQ=255)

                                                            tATTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAt                            <  1:615239/93‑1 (MQ=255)

                                                             aTTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                           >  2:4026034/1‑93 (MQ=255)

                                                             aTTTCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGAtt                           >  2:3815528/1‑93 (MQ=255)

                                                             aTTTCATGCAGCCATTACACGACGAATGCCTATGAACTCCAGTCACTCCGCCATATTTTTGCCATGACTATTGGTGAATGCGCGACCTGGAtg                           >  1:322817/1‑92 (MQ=255)

                                                              tttCATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTg                          >  2:3412683/1‑93 (MQ=255)

                                                                tCAGGCAGCCATTACACGATGAGCGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGcc                        <  1:1912297/93‑1 (MQ=255)

                                                                 caagCAGCCATTGCACGATGAACGCCTACGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                       <  2:712809/90‑1 (MQ=255)

                                                                 cATGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGccc                       >  1:2652294/1‑93 (MQ=255)

                                                                  aTGCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCa                      <  1:4692130/93‑1 (MQ=255)

                                                                    gCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                    >  1:22293/1‑93 (MQ=255)

                                                                    gCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                    >  1:4512827/1‑93 (MQ=255)

                                                                    gCAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGa                    <  2:3880821/93‑1 (MQ=255)

                                                                     cAGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAc                   <  2:152146/93‑1 (MQ=255)

                                                                      aGCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGAcc                  >  2:587247/1‑93 (MQ=255)

                                                                      aGCCATTACACGATGAACGCCTATGAACTCCACGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCCGAcc                  >  1:4311259/1‑93 (MQ=255)

                                                                       gCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                 >  1:549140/1‑93 (MQ=255)

                                                                       gCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                 >  1:3276465/1‑93 (MQ=255)

                                                                       gCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                 >  1:558467/1‑93 (MQ=255)

                                                                       gCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                 >  1:3276464/1‑93 (MQ=255)

                                                                       gCCATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCg                 >  2:158884/1‑93 (MQ=255)

                                                                        ccATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg                >  2:4000364/1‑93 (MQ=255)

                                                                        ccATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCgg                <  1:2580477/93‑1 (MQ=255)

                                                                         cATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGTTTGCCCAGCCCGGt               >  1:1468608/1‑93 (MQ=255)

                                                                         cATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt               >  1:1468459/1‑93 (MQ=255)

                                                                         cATTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGt               <  2:4059283/93‑1 (MQ=255)

                                                                          aTTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTg              >  2:2379699/1‑93 (MQ=255)

                                                                          aTTACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCACAGACCGGTg              >  2:2379696/1‑93 (MQ=255)

                                                                           ttACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGa             >  1:3218303/1‑93 (MQ=255)

                                                                           ttACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGa             >  1:863826/1‑93 (MQ=255)

                                                                            tACACGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGAc            >  2:2646257/1‑93 (MQ=255)

                                                                             acacGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACa           >  1:2696780/1‑93 (MQ=255)

                                                                              cacGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg          >  2:259642/1‑93 (MQ=255)

                                                                              cacGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAg          >  2:1093641/1‑93 (MQ=255)

                                                                               acGATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGt         >  1:1608193/1‑93 (MQ=255)

                                                                                 gATGAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGa       >  1:4177232/1‑93 (MQ=255)

                                                                                    gAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg    >  2:2289616/1‑93 (MQ=255)

                                                                                    gAACGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAg    <  1:1285191/93‑1 (MQ=255)

                                                                                      aCGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACGTGGATTGCCCAGACCGGTGACAGTGAAAGCt  <  1:231175/93‑1 (MQ=255)

                                                                                      aCGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCt  >  1:2497102/1‑93 (MQ=255)

                                                                                      aCGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCt  <  2:1730975/93‑1 (MQ=255)

                                                                                      aCGCCTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCt  >  1:1958252/1‑93 (MQ=255)

GATACTGTGACTCTCTTCAAGCCGTCAGCAGGCTGATAATGTTATTTATATTTTATTGTTTATTTCATGCAGCCATTACACGATGAACGCTTATGAACTCCAGGCACTCCGCCATATTTTTGCCATGACTATTGATGAATGCGCGACCTGGATTGCCCAGACCGGTGACAGTGAAAGCT  >  NC_000913/1768557‑1768735

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