BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,021,662	A→G	V485A (GTG→GCG)	b0960 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,021,662	0	A	G	88.1%	674.8 / 77.7	268	V485A (GTG→GCG)	b0960	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (16/16); new base G (134/102); total (150/118)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 268 total aligned reads displayed.

GAGGTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGC  >  NC_000913/1021573‑1021752

gAGGTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGt                                                                                         >  1:2625460/1‑93 (MQ=255)

  ggTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTa                                                                                       <  1:5362473/93‑1 (MQ=255)

   gTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTAc                                                                                      <  2:1420718/93‑1 (MQ=255)

     aaCGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACgg                                                                                    >  1:1459171/1‑93 (MQ=255)

      aCGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGa                                                                                   >  1:3990833/1‑93 (MQ=255)

      aCGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGa                                                                                   <  1:7393321/93‑1 (MQ=255)

      aCGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGa                                                                                   >  2:3153157/1‑93 (MQ=255)

      aCGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGa                                                                                   <  1:4270158/93‑1 (MQ=255)

      aCGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGa                                                                                   <  1:94535/93‑1 (MQ=255)

       cGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGaa                                                                                  >  2:6612592/1‑93 (MQ=255)

       cGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGaa                                                                                  <  2:3294678/93‑1 (MQ=255)

        gCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAg                                                                                 >  1:5447967/1‑93 (MQ=255)

         cTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAgg                                                                                >  2:4785266/1‑93 (MQ=255)

         cTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAgg                                                                                >  1:4876719/1‑93 (MQ=255)

         cTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAgg                                                                                >  2:609733/1‑93 (MQ=255)

           aCTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCa                                                                              >  2:7288528/1‑93 (MQ=255)

           aCTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCa                                                                              >  2:6021311/1‑93 (MQ=255)

           aCTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCa                                                                              >  2:7187572/1‑93 (MQ=255)

            cTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaa                                                                             <  2:352572/93‑1 (MQ=255)

            cTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaa                                                                             >  1:4477600/1‑93 (MQ=255)

             ttCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaaa                                                                            <  2:2695630/93‑1 (MQ=255)

             ttCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaaa                                                                            <  1:186106/93‑1 (MQ=255)

               cAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaaaaa                                                                          <  2:3996627/93‑1 (MQ=255)

                aaaaCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaaaaaa                                                                         >  2:6824700/1‑93 (MQ=255)

                aaaaCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCaaaaaa                                                                         <  2:4969301/93‑1 (MQ=255)

                 aaaCCCTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAag                                                                        <  1:3186357/93‑1 (MQ=255)

                  aaCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAaga                                                                       <  1:2993153/93‑1 (MQ=255)

                  aaCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAaga                                                                       >  1:4645724/1‑93 (MQ=255)

                  aaCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAaga                                                                       <  1:4570649/93‑1 (MQ=255)

                   aCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAagag                                                                      <  1:3517128/93‑1 (MQ=255)

                   aCCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAagag                                                                      >  2:3366283/1‑93 (MQ=255)

                    ccTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGc                                                                     <  1:87651/93‑1 (MQ=255)

                     cTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCa                                                                    <  1:5477391/93‑1 (MQ=255)

                      ttCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCAc                                                                   >  1:5246060/1‑93 (MQ=255)

                      ttCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCAc                                                                   >  2:2563919/1‑93 (MQ=255)

                         aCCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGcc                                                                <  1:66777/93‑1 (MQ=255)

                         aCCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGcc                                                                <  2:5721549/93‑1 (MQ=255)

                          cccAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCg                                                               <  1:3332398/93‑1 (MQ=255)

                          cccAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCg                                                               <  1:3302022/93‑1 (MQ=255)

                           ccAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCgg                                                              >  1:2388763/1‑93 (MQ=255)

                           ccAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCgg                                                              >  2:2002302/1‑93 (MQ=255)

                             aGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGta                                                            >  2:1014476/1‑93 (MQ=255)

                              gTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGtaa                                                           >  1:4962914/1‑93 (MQ=255)

                               tAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGtaat                                                          >  1:2750393/1‑93 (MQ=255)

                                aaGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGtaata                                                         >  1:3460309/1‑93 (MQ=255)

                                 aGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGtaataa                                                        <  1:7108203/93‑1 (MQ=255)

                                  gTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAAc                                                       >  1:1381376/1‑93 (MQ=255)

                                    tAAACCACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGAAAGGCAAAAAAGAGCACGCCGGTAATAACca                                                     <  1:2521325/93‑1 (MQ=255)

                                      aaaCACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACcagc                                                   >  1:3231924/1‑93 (MQ=255)

                                         cacaGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCagcagc                                                <  1:1913181/93‑1 (MQ=255)

                                            aGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCAcc                                             <  1:1183543/93‑1 (MQ=255)

                                             gTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCa                                            <  1:681948/93‑1 (MQ=255)

                                              tAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAg                                           >  2:5875101/1‑93 (MQ=255)

                                              tAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAg                                           >  2:777111/1‑93 (MQ=255)

                                               aGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGc                                          >  1:3810745/1‑93 (MQ=255)

                                               aGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGc                                          <  1:6115402/93‑1 (MQ=255)

                                                gCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCt                                         >  1:1465905/1‑93 (MQ=255)

                                                  aCCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCtnn                                       >  2:2168947/1‑91 (MQ=255)

                                                   ccAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGcc                                      >  2:2309091/1‑93 (MQ=255)

                                                    cAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGccc                                     <  2:1227025/93‑1 (MQ=255)

                                                     aaaaGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCt                                    >  1:5954955/1‑93 (MQ=255)

                                                      aaaGTGTGATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCtt                                   >  1:3541621/1‑93 (MQ=255)

                                                      aaaGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCtt                                   >  1:3706981/1‑93 (MQ=255)

                                                        aGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTcc                                 <  1:6954216/93‑1 (MQ=255)

                                                          tgtgATGAACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAg                               >  1:2137066/1‑93 (MQ=255)

                                                           gtgATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGt                              <  2:6157493/93‑1 (MQ=255)

                                                            tgatgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAggg                             >  1:6238162/1‑91 (MQ=255)

                                                            tgatgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTg                             >  2:6297848/1‑93 (MQ=255)

                                                            tgatgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAAAAACCAGCAGCACCAGCTGCCCTTCCAGTg                             >  2:6499698/1‑93 (MQ=255)

                                                             gatgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGa                            <  2:175254/93‑1 (MQ=255)

                                                             gatgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGa                            >  2:3795729/1‑93 (MQ=255)

                                                               tgaACATCGTTGCATGAGCTTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATg                          <  1:6339573/93‑1 (MQ=255)

                                                               tgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATg                          <  2:1623059/93‑1 (MQ=255)

                                                               tgaACATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATg                          >  2:529439/1‑93 (MQ=255)

                                                                 aaCATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGc                        <  1:1022004/93‑1 (MQ=255)

                                                                 aaCATCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGc                        >  2:5628294/1‑93 (MQ=255)

                                                                    aTCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGcaca                     >  2:4158389/1‑93 (MQ=255)

                                                                     tCGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACaa                    >  2:1663639/1‑93 (MQ=255)

                                                                      cGTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACaaa                   <  2:4647778/93‑1 (MQ=255)

                                                                       gTTGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAAc                  <  1:6712322/93‑1 (MQ=255)

                                                                        ttGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAAcc                 <  1:309435/93‑1 (MQ=255)

                                                                         tGCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACca                <  1:774946/93‑1 (MQ=255)

                                                                          gCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcac               <  1:6450558/93‑1 (MQ=255)

                                                                          gCATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcac               <  2:3236387/93‑1 (MQ=255)

                                                                           cATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcaca              >  2:872909/1‑93 (MQ=255)

                                                                           cATGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcaca              >  2:3213965/1‑93 (MQ=255)

                                                                            aTGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacag             <  1:150768/93‑1 (MQ=255)

                                                                             tGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagc            <  2:6372599/93‑1 (MQ=255)

                                                                             tGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagc            <  2:2901681/93‑1 (MQ=255)

                                                                             tGAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagc            >  2:7029173/1‑93 (MQ=255)

                                                                              nngCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagca           <  2:3570965/91‑1 (MQ=255)

                                                                              gAGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagca           <  2:725452/93‑1 (MQ=255)

                                                                               aGCGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagcac          >  2:1442157/1‑93 (MQ=255)

                                                                                 cGTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACcacagcacag        <  1:6797583/93‑1 (MQ=255)

                                                                                  gTATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAgg       <  2:7111943/93‑1 (MQ=255)

                                                                                   tATTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGa      <  1:4765857/93‑1 (MQ=255)

                                                                                    aTTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCGCAAACCACAGCACAGGaa     >  2:2419698/1‑93 (MQ=255)

                                                                                    aTTGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGaa     <  1:726222/93‑1 (MQ=255)

                                                                                     ttGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAAt    <  1:3028468/93‑1 (MQ=255)

                                                                                      tGCGCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATg   <  1:4119538/93‑1 (MQ=255)

GAGGTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGC  >  NC_000913/1021573‑1021752

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