BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,049,149	G→A	G75S (GGC→AGC)	yihI →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,049,149	0	G	A	99.6%	956.4 / ‑5.6	254	G75S (GGC→AGC)	yihI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (120/132); total (122/132)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 257 total aligned reads displayed.

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

cgcgCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGc                                                                                         <  2:4052515/93‑1 (MQ=255)

cgcgCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGc                                                                                         <  1:5247099/93‑1 (MQ=255)

  cgcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGt                                                                                       <  1:2603540/93‑1 (MQ=255)

   gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                                      >  2:3106506/1‑93 (MQ=255)

   gcAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTg                                                                                      >  2:597362/1‑93 (MQ=255)

     aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGAc                                                                                    <  2:1716882/93‑1 (MQ=255)

     aGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGAc                                                                                    <  1:5504257/93‑1 (MQ=255)

      gcgggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACt                                                                                   >  1:1111510/1‑93 (MQ=255)

      gcgggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACt                                                                                   >  2:3523259/1‑93 (MQ=255)

      gcgggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACt                                                                                   >  1:3488727/1‑93 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 >  1:4840147/1‑93 (MQ=255)

        gggcggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGa                                                                                 <  1:1916772/93‑1 (MQ=255)

           cggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaa                                                                              <  1:4095490/93‑1 (MQ=255)

            ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                             <  1:3905369/93‑1 (MQ=255)

            ggAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGaaaaa                                                                             <  1:1349293/93‑1 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGCTCCACGTATTGGCAGTAAAACACCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            >  2:1609992/1‑93 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            >  1:2020818/1‑93 (MQ=255)

             gAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAg                                                                            <  1:5788809/93‑1 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         >  2:424369/1‑93 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         >  1:3951001/1‑93 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         >  1:310659/1‑93 (MQ=255)

                acacCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCa                                                                         <  1:6141702/93‑1 (MQ=255)

                  accacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCAcc                                                                       <  1:2396376/93‑1 (MQ=255)

                   ccacGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCa                                                                      >  1:6016942/1‑93 (MQ=255)

                     acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                                    >  1:439834/1‑93 (MQ=255)

                     acGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCaaa                                                                    >  1:226472/1‑93 (MQ=255)

                      cGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAAc                                                                   <  1:4651727/93‑1 (MQ=255)

                      cGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAAc                                                                   <  1:6104493/93‑1 (MQ=255)

                       gTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACa                                                                  >  1:2233691/1‑93 (MQ=255)

                        tCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAg                                                                 <  2:773510/93‑1 (MQ=255)

                         cAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGc                                                                <  1:2311736/93‑1 (MQ=255)

                          aGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGca                                                               <  2:1010232/93‑1 (MQ=255)

                           ggCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGcac                                                              <  2:2331436/93‑1 (MQ=255)

                             cagcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACaa                                                            <  2:4314947/93‑1 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          >  2:2299995/1‑93 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          >  1:2871996/1‑93 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAAc                                                          >  1:6194798/1‑93 (MQ=255)

                               gcaAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCCGTTCCATTGAGCGTGAATGAAAAAGTCACCAAACAGCACAAAc                                                          >  1:5927333/1‑93 (MQ=255)

                                 aaaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCg                                                        >  2:5822223/1‑93 (MQ=255)

                                  aaGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGa                                                       >  1:3759337/1‑93 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      >  2:3789257/1‑93 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      <  2:1437579/93‑1 (MQ=255)

                                   aGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGaa                                                      <  1:2650295/93‑1 (MQ=255)

                                    ggCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAag                                                     <  2:4872789/93‑1 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    <  2:4381403/93‑1 (MQ=255)

                                     gCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAaga                                                    >  2:5360945/1‑93 (MQ=255)

                                      ccAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAagag                                                   >  1:1736730/1‑93 (MQ=255)

                                       cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                                  >  2:694469/1‑93 (MQ=255)

                                       cAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGt                                                  >  1:3604846/1‑93 (MQ=255)

                                        aGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTg                                                 >  1:5238383/1‑93 (MQ=255)

                                          aaCGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTgag                                               <  2:4502859/93‑1 (MQ=255)

                                             gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                            <  2:332697/93‑1 (MQ=255)

                                             gCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGaaa                                            >  2:2419759/1‑93 (MQ=255)

                                              cACCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAc                                           >  2:4819239/1‑93 (MQ=255)

                                               aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                          <  2:410912/93‑1 (MQ=255)

                                               aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                          <  2:499981/93‑1 (MQ=255)

                                               aCCAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAAcc                                          >  2:5033788/1‑93 (MQ=255)

                                                ccAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCt                                         >  2:5742483/1‑93 (MQ=255)

                                                 cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                                        <  1:4002103/93‑1 (MQ=255)

                                                 cAAAAGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTa                                        <  1:1791345/93‑1 (MQ=255)

                                                  aaaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTAt                                       <  1:143988/93‑1 (MQ=255)

                                                   aaaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATg                                      <  1:907610/93‑1 (MQ=255)

                                                    aaGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGc                                     <  1:5664821/93‑1 (MQ=255)

                                                     aGATCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCt                                    <  1:1579946/93‑1 (MQ=255)

                                                       aTCCACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCttt                                  <  1:4234706/93‑1 (MQ=255)

                                                          cACGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCAc                               >  2:3032301/1‑93 (MQ=255)

                                                            cGTATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCg                             <  1:3788913/93‑1 (MQ=255)

                                                              tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                           <  2:618593/93‑1 (MQ=255)

                                                              tATTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCa                           >  2:2526028/1‑93 (MQ=255)

                                                               aTTGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAg                          >  2:4128903/1‑93 (MQ=255)

                                                                ttGGCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAgg                         <  1:6075281/93‑1 (MQ=255)

                                                                  ggCAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCg                       <  2:4538551/93‑1 (MQ=255)

                                                                    cAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGa                     >  1:2152415/1‑93 (MQ=255)

                                                                    cAGTAAAACCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGa                     >  2:3607415/1‑93 (MQ=255)

                                                                        aaaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTg                 <  2:5291919/93‑1 (MQ=255)

                                                                         aaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTgg                <  1:1513197/93‑1 (MQ=255)

                                                                         aaaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTgg                <  1:2532127/93‑1 (MQ=255)

                                                                          aaCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGa               >  2:4611536/1‑93 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCTAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  1:3389559/93‑1 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  1:105696/93‑1 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  1:5619119/93‑1 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  1:2908340/93‑1 (MQ=255)

                                                                           aCCCCTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGaa              <  2:2806173/93‑1 (MQ=255)

                                                                              ccTGTTCCATTGCGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa           >  2:1925764/1‑93 (MQ=255)

                                                                              ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa           >  2:5511144/1‑93 (MQ=255)

                                                                              ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa           >  2:4767634/1‑93 (MQ=255)

                                                                              ccTGTTCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTa           >  2:5192235/1‑93 (MQ=255)

                                                                               ctgttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTAc          <  1:5042462/90‑1 (MQ=255)

                                                                                tgttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACt         <  2:2691503/91‑1 (MQ=255)

                                                                                 gttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTg        >  1:5302497/2‑93 (MQ=255)

                                                                                 gttCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTg        >  2:5571527/2‑93 (MQ=255)

                                                                                  ttCCATCGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTgg       <  2:1249460/93‑1 (MQ=255)

                                                                                   tCCATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGa      <  2:4652939/93‑1 (MQ=255)

                                                                                     cATTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGaaa    >  2:271351/1‑93 (MQ=255)

                                                                                      aTTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAAc   >  2:1059578/1‑93 (MQ=255)

                                                                                      aTTGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAAc   <  2:872190/93‑1 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACg  >  2:5906794/1‑93 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACg  >  1:1711648/1‑93 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACg  >  1:4081895/1‑93 (MQ=255)

                                                                                       ttGAGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACg  <  2:3666987/93‑1 (MQ=255)

CGCGCAGCGGGCGGAAACACCACGTCAGGCAGCAAAGGCCAGAACGCACCAAAAGATCCACGTATTGGCAGTAAAACCCCTATTCCATTGGGCGTGACTGAAAAAGTCACCAAACAGCACAAACCGAAGAGTGAGAAACCTATGCTTTCACCGCAGGCGGAGTTGGAATTACTGGAAACG  >  NC_000913/4049059‑4049238

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