BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,476,619	A→G	V91A (GTG→GCG)	b4256 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,476,619	0	A	G	99.5%	629.3 / ‑5.2	203	V91A (GTG→GCG)	b4256	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (99/103); total (99/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 206 total aligned reads displayed.

ATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAA  >  NC_000913/4476529‑4476708

  ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTc                                                                                       <  2:1073501/93‑1 (MQ=255)

  ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTc                                                                                       >  1:1550421/1‑93 (MQ=255)

   cGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTca                                                                                      <  1:3263746/93‑1 (MQ=255)

   cGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTca                                                                                      <  1:3524800/93‑1 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                 <  2:3737548/93‑1 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                 >  1:3120417/1‑93 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                 >  2:1560781/1‑93 (MQ=255)

        tAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATc                                                                                 >  2:2670213/1‑93 (MQ=255)

         aaaCAAAACGTACGTGGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCg                                                                                <  1:1175107/93‑1 (MQ=255)

          aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGa                                                                               >  1:147916/1‑93 (MQ=255)

          aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGa                                                                               <  2:3215592/93‑1 (MQ=255)

            cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACa                                                                             >  1:1532047/1‑93 (MQ=255)

            cAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACa                                                                             <  2:4292560/93‑1 (MQ=255)

             aaaaCGTACGTTGTTAATGTTTGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                            >  2:79984/1‑93 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTTGCGCAACATCCGCGCTTCACAATCGACAt                                                                            >  2:547355/1‑93 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                            <  2:1821920/93‑1 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACAt                                                                            <  2:3750465/93‑1 (MQ=255)

             aaaaCGTACGTTGTTAATGTTCGAAGGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTCCGCGCAACATCCGCGCTTCACAATCGACAt                                                                            <  2:761896/93‑1 (MQ=255)

              aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATg                                                                           >  2:413717/1‑93 (MQ=255)

              aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATg                                                                           <  2:4530845/93‑1 (MQ=255)

               aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATgg                                                                          <  1:2292729/93‑1 (MQ=255)

               aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATgg                                                                          <  1:1681920/93‑1 (MQ=255)

                 cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGcc                                                                        <  1:3831023/93‑1 (MQ=255)

                 cgtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGcc                                                                        >  1:2141637/1‑93 (MQ=255)

                  gtacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGccc                                                                       >  2:3097360/1‑93 (MQ=255)

                   tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCTCGCAACATCCGCGCTTCACAATCGACATGGCCCg                                                                      >  2:2095744/1‑93 (MQ=255)

                    acgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGt                                                                     <  2:3318770/93‑1 (MQ=255)

                    acgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGAGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGt                                                                     >  1:2014502/1‑93 (MQ=255)

                     cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTg                                                                    >  2:3075902/1‑93 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                   >  1:3577215/1‑93 (MQ=255)

                      gttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGc                                                                   >  2:3701841/1‑93 (MQ=255)

                       ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCa                                                                  >  1:2250123/1‑93 (MQ=255)

                       ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCa                                                                  >  1:1375074/1‑93 (MQ=255)

                       ttgttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCa                                                                  >  1:1204938/1‑93 (MQ=255)

                         gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGc                                                                <  2:2974718/93‑1 (MQ=255)

                         gttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGc                                                                >  2:1205191/1‑93 (MQ=255)

                           tAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCa                                                              >  1:1008788/1‑93 (MQ=255)

                            aaTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGATTCACAATCGACATGGCCCGTGCAGCCaa                                                             >  2:4604222/1‑93 (MQ=255)

                             aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAg                                                            >  1:2410601/1‑93 (MQ=255)

                             aTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAg                                                            >  1:4379493/1‑93 (MQ=255)

                                ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAgcgc                                                         >  2:362100/1‑93 (MQ=255)

                                 tCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCa                                                        <  2:448291/93‑1 (MQ=255)

                                   gAATGCCCGCGAGTTGAACCCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATa                                                      <  1:870717/93‑1 (MQ=255)

                                   gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGCGCAGCCAAgcgacta                                                      >  1:913582/1‑89 (MQ=255)

                                    aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAg                                                     <  1:4080123/93‑1 (MQ=255)

                                      tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                   <  2:3944652/93‑1 (MQ=255)

                                      tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                   >  1:1811651/1‑93 (MQ=255)

                                        cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                 >  2:3052800/1‑93 (MQ=255)

                                        cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                 <  1:3932989/93‑1 (MQ=255)

                                         ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGat                                                <  1:96956/93‑1 (MQ=255)

                                         ccGCGAGTTAAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGaa                                                >  2:4039878/1‑92 (MQ=255)

                                          cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                               <  1:748747/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                              <  1:2827626/93‑1 (MQ=255)

                                           gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGatat                                              >  1:197745/1‑93 (MQ=255)

                                             gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                            <  1:2580389/93‑1 (MQ=255)

                                             gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                            <  1:452672/93‑1 (MQ=255)

                                              aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                                           <  1:1365777/93‑1 (MQ=255)

                                               gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                                          >  2:2005772/1‑93 (MQ=255)

                                               gTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTc                                          >  2:2827904/1‑93 (MQ=255)

                                                ttGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGGCAAGCGCATAGTCGATATGCTCa                                         >  1:1751569/1‑93 (MQ=255)

                                                ttGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCa                                         <  2:1463357/93‑1 (MQ=255)

                                                 tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCCTAGTGGATATTCTCaa                                        >  2:1414183/1‑93 (MQ=255)

                                                   aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                      <  2:3165217/93‑1 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                     <  2:162920/93‑1 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                     >  1:4415929/1‑93 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                     <  2:1891632/93‑1 (MQ=255)

                                                    acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                     <  2:2047790/93‑1 (MQ=255)

                                                     cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGcc                                    >  1:264744/1‑93 (MQ=255)

                                                         cgGGCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCaaa                                <  2:2998251/93‑1 (MQ=255)

                                                            gCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc                             <  1:3043650/93‑1 (MQ=255)

                                                            gCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc                             <  1:3674065/93‑1 (MQ=255)

                                                            gCATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGc                             <  1:4622314/93‑1 (MQ=255)

                                                             cATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                            <  1:2682489/93‑1 (MQ=255)

                                                             cATTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                            >  2:708745/1‑93 (MQ=255)

                                                              aTTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                           >  1:2170490/1‑93 (MQ=255)

                                                              aTTTTTTTAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                           >  2:4438870/1‑93 (MQ=255)

                                                                  ttttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaa                       <  1:2968901/93‑1 (MQ=255)

                                                                   tttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  2:4428718/93‑1 (MQ=255)

                                                                   tgtAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                      <  2:1795598/91‑1 (MQ=255)

                                                                    ttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAg                     >  1:3403776/1‑93 (MQ=255)

                                                                    ttAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAg                     >  2:1654640/1‑93 (MQ=255)

                                                                     tAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgc                    <  1:4337158/93‑1 (MQ=255)

                                                                     tAGAGTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgc                    <  2:1544921/93‑1 (MQ=255)

                                                                       gagTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgcgc                  >  2:3505020/1‑93 (MQ=255)

                                                                         gTTCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCaa                <  2:4573653/93‑1 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt               <  1:1876797/93‑1 (MQ=255)

                                                                          ttCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAAt               >  1:4365269/1‑93 (MQ=255)

                                                                           tCGCGCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATg              <  1:651704/93‑1 (MQ=255)

                                                                            cgcgCAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATgg             <  1:2817813/93‑1 (MQ=255)

                                                                             gcgcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGc            >  1:2376041/1‑93 (MQ=255)

                                                                              cgcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCt           <  1:3692643/93‑1 (MQ=255)

                                                                               gcAACATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCtt          >  1:3387699/1‑93 (MQ=255)

                                                                                 aaCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTcc        >  1:4032174/1‑93 (MQ=255)

                                                                                 aaCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTcc        <  2:591761/93‑1 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       <  2:2375076/93‑1 (MQ=255)

                                                                                  aCATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt       >  1:1231130/1‑93 (MQ=255)

                                                                                   cATCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCtt      >  2:801432/1‑93 (MQ=255)

                                                                                    aTCCGCGCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCttt     <  2:2410584/93‑1 (MQ=255)

                                                                                       cgcgCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaa  >  1:4154708/1‑93 (MQ=255)

                                                                                       cgcgCTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaa  >  1:4252961/1‑93 (MQ=255)

ATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAA  >  NC_000913/4476529‑4476708

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