BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,478,212	A→G	S351G (AGT→GGT)	b1410 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,478,212	0	A	G	99.0%	291.8 / NA	97	S351G (AGT→GGT)	b1410	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (42/54); total (43/54)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCAGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGTTGGCTG  >  NC_000913/1478123‑1478303

cGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGtt                                                                                          >  1:2250882/1‑93 (MQ=255)

 gTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTc                                                                                         <  2:4614300/93‑1 (MQ=255)

 gTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTc                                                                                         >  2:4557847/1‑93 (MQ=255)

 gTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTc                                                                                         >  2:4161250/1‑93 (MQ=255)

  tATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCg                                                                                        >  1:1599352/1‑93 (MQ=255)

   aTCGCTTTATGCGAAGGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGc                                                                                       <  1:3816738/93‑1 (MQ=255)

   aTCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGc                                                                                       <  2:330204/93‑1 (MQ=255)

    tCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCt                                                                                      >  2:1468164/1‑93 (MQ=255)

     cGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTg                                                                                     >  2:4469627/1‑93 (MQ=255)

     cGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTg                                                                                     <  2:335878/93‑1 (MQ=255)

        tttATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGAt                                                                                  >  1:2411922/1‑93 (MQ=255)

        tttATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGCGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGAt                                                                                  >  1:2787592/1‑93 (MQ=255)

         ttATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGAtt                                                                                 >  1:1512551/1‑93 (MQ=255)

           aTGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTAt                                                                               >  1:1478005/1‑93 (MQ=255)

              cGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTc                                                                            <  1:2457035/93‑1 (MQ=255)

                aaaGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCGGGCGGTTCGCTGGATTATGTCTa                                                                          <  1:4377510/93‑1 (MQ=255)

                 aaGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTAt                                                                         <  2:3292021/93‑1 (MQ=255)

                  aGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATc                                                                        <  2:1652820/93‑1 (MQ=255)

                  aGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATc                                                                        <  1:2845177/93‑1 (MQ=255)

                    gCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCaa                                                                      <  2:3923350/93‑1 (MQ=255)

                      gATGAAATTGTTCTGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCaaaa                                                                    <  1:3176795/93‑1 (MQ=255)

                        tGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCATGCGGTTCGCTGGATTATGTCTATCAAAATc                                                                  >  1:1571832/1‑93 (MQ=255)

                        tGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATc                                                                  >  1:2210367/1‑93 (MQ=255)

                         gAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCa                                                                 >  2:998201/1‑93 (MQ=255)

                          aaaTTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCaa                                                                >  2:3937503/1‑93 (MQ=255)

                           aaTTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAAc                                                               <  2:3414351/93‑1 (MQ=255)

                            aTTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAAcc                                                              <  1:2303900/93‑1 (MQ=255)

                            aTTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAAcc                                                              >  1:4020122/1‑93 (MQ=255)

                            aTTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAAcc                                                              <  1:394288/93‑1 (MQ=255)

                             ttGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGAATATGTCTATCAAAATCAACCg                                                             >  2:4078284/1‑93 (MQ=255)

                              tGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGc                                                            <  2:3999298/93‑1 (MQ=255)

                              tGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGc                                                            >  1:2243204/1‑93 (MQ=255)

                              tGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGc                                                            <  1:3023202/93‑1 (MQ=255)

                               gTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCa                                                           >  1:2735144/1‑93 (MQ=255)

                                ttCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCaa                                                          <  2:4056913/93‑1 (MQ=255)

                                 tCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAg                                                         <  1:2287178/93‑1 (MQ=255)

                                 tCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAg                                                         <  1:4164564/93‑1 (MQ=255)

                                  cGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAgg                                                        <  1:277812/93‑1 (MQ=255)

                                   gggACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAgga                                                       >  1:330189/1‑92 (MQ=255)

                                    ggACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGgta                                                      <  1:4661676/93‑1 (MQ=255)

                                    ggACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGgta                                                      >  1:2892193/1‑93 (MQ=255)

                                      aCGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGgtagt                                                    >  2:2254188/1‑93 (MQ=255)

                                        gCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTaa                                                  <  2:1486594/93‑1 (MQ=255)

                                        gCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTaa                                                  >  2:4205875/1‑93 (MQ=255)

                                        gCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTaa                                                  >  1:655352/1‑93 (MQ=255)

                                         cACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAAc                                                 <  2:2104695/93‑1 (MQ=255)

                                          aCTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACg                                                >  1:162380/1‑93 (MQ=255)

                                            tctTCGGGCCTGCATCTCGGAATGAGCACCGTCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCa                                              <  1:613350/93‑1 (MQ=255)

                                            tctTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCa                                              >  1:3184674/1‑93 (MQ=255)

                                            tctTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCa                                              <  2:2058118/93‑1 (MQ=255)

                                             ttgcGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTGCGATGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCAt                                             <  2:4554481/90‑1 (MQ=255)

                                              ttCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCAtt                                            >  1:1528954/1‑93 (MQ=255)

                                              ttCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCAtt                                            <  2:897001/93‑1 (MQ=255)

                                                cGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCg                                          <  2:4122993/93‑1 (MQ=255)

                                                  ggCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCggg                                        <  2:1058831/93‑1 (MQ=255)

                                                  ggCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATAAACCGCAAGGTAGTAACGCACTCggg                                        >  1:1478966/1‑93 (MQ=255)

                                                  ggCCTGCATCTCGGAATGAGAACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCggg                                        >  1:2824486/1‑93 (MQ=255)

                                                   gCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGc                                       >  1:387078/1‑93 (MQ=255)

                                                   gCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGc                                       <  1:3006151/93‑1 (MQ=255)

                                                   gCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGGTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGc                                       >  1:2394347/1‑93 (MQ=255)

                                                      tGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGtt                                    >  2:3878608/1‑93 (MQ=255)

                                                       gCATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGttt                                   <  2:3892321/93‑1 (MQ=255)

                                                        cATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTa                                  <  1:496892/93‑1 (MQ=255)

                                                        cATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTa                                  <  1:1183465/93‑1 (MQ=255)

                                                        cATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTa                                  <  1:84406/93‑1 (MQ=255)

                                                        cATCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTa                                  >  1:3999567/1‑93 (MQ=255)

                                                         aTCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAg                                 >  1:2167079/1‑93 (MQ=255)

                                                         aTCTCGGAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAg                                 <  1:1078262/93‑1 (MQ=255)

                                                               gAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACa                           >  2:1967636/1‑93 (MQ=255)

                                                               gAATGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACa                           <  1:116374/93‑1 (MQ=255)

                                                                 aTGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACaaa                         <  2:3468791/93‑1 (MQ=255)

                                                                  tGAGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACaaaa                        <  2:1692838/93‑1 (MQ=255)

                                                                    aGCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATc                      <  1:2842031/93‑1 (MQ=255)

                                                                     gCACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCt                     <  1:1417127/93‑1 (MQ=255)

                                                                      cACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTa                    <  1:1047218/93‑1 (MQ=255)

                                                                      cACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTa                    <  2:497828/93‑1 (MQ=255)

                                                                      cACCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTa                    <  2:2574863/93‑1 (MQ=255)

                                                                       aCCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTAc                   <  1:3355116/93‑1 (MQ=255)

                                                                       aCCGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTAc                   >  2:1018257/1‑93 (MQ=255)

                                                                        ccGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTAcc                  <  1:245776/93‑1 (MQ=255)

                                                                         cGGCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCt                 >  1:1775990/1‑93 (MQ=255)

                                                                          ggCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTg                >  1:1848015/1‑93 (MQ=255)

                                                                          ggCTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTg                >  1:528982/1‑93 (MQ=255)

                                                                            cTTTGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGaa              <  1:551716/93‑1 (MQ=255)

                                                                             tttGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAAc             >  2:2697517/1‑93 (MQ=255)

                                                                             tttGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAAc             >  1:2687704/1‑93 (MQ=255)

                                                                             tttGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAAc             >  1:2244225/1‑93 (MQ=255)

                                                                             tttGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAAc             <  2:1309852/93‑1 (MQ=255)

                                                                              ttGATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTCTAGTCGACAAAATCTACCTGAACa            <  2:4517983/93‑1 (MQ=255)

                                                                                gATTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAgt          >  1:1876285/1‑93 (MQ=255)

                                                                                 aTTCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAgtg         <  1:205224/93‑1 (MQ=255)

                                                                                  ttCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAgtgt        <  2:331510/93‑1 (MQ=255)

                                                                                   tCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGtt       >  2:1096363/1‑93 (MQ=255)

                                                                                   tCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGtt       <  2:4411792/93‑1 (MQ=255)

                                                                                   tCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGtt       <  1:2921169/93‑1 (MQ=255)

                                                                                   tCAGGCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGtt       <  2:1394507/93‑1 (MQ=255)

                                                                                      ggCGGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGTtggc    <  1:1332848/93‑1 (MQ=255)

                                                                                        cAGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGTtggctg  >  2:3719424/1‑93 (MQ=255)

CGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCTCGGAATGAGCACCGGCTTTGATTCAGGCAGTTCGCTGGATTATGTCTATCAAAATCAACCGCAAGGTAGTAACGCATTCGGGCGTTTAGTCGACAAAATCTACCTGAACAGTGTTGGCTG  >  NC_000913/1478123‑1478303

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