BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,169,304	C→A	13.4%	P73Q (CCG→CAG)	tufA →	protein chain elongation factor EF‑Tu

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_007779	4,169,304	0	C→A	13.4%	29.5 ‑4.0 ‑0.0	104	P73Q (CCG→CAG)	tufA	protein chain elongation factor EF‑Tu
Reads supporting (aligned to +/- strand): new base (2/12): ref base (45/45): total (47/57)
Fisher's exact test strand distribution p-value = 1.87e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.80e-01

Only 100 of 114 total aligned reads displayed.

TACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTT  >  NC_007779/4169204‑4169404

tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                                      <  1:13746992/101‑1 (MQ=32)

tACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACAcccc                                                                                                      <  1:6566669/101‑1 (MQ=32)

 aCGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCg                                                                                                     >  1:1889069/1‑101 (MQ=32)

 aCGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCg                                                                                                     >  1:7478188/1‑101 (MQ=32)

  cggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGa                                                                                                    >  1:8834744/1‑101 (MQ=32)

  cggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGa                                                                                                    <  1:11039692/101‑1 (MQ=32)

   ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                                   >  1:1243737/1‑101 (MQ=32)

   ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                                   >  1:487070/1‑101 (MQ=33)

   ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                                   >  1:8640225/1‑101 (MQ=32)

   ggcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAc                                                                                                   <  1:5276358/101‑1 (MQ=32)

    gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                                  <  1:2737201/101‑1 (MQ=32)

    gcggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGAcc                                                                                                  >  1:1403550/1‑101 (MQ=32)

      ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                                <  1:5843119/101‑1 (MQ=32)

      ggTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCg                                                                                                >  1:4912007/1‑101 (MQ=32)

       gTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGt                                                                                               >  1:3945278/1‑101 (MQ=32)

        tgctgcTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTc                                                                                              <  1:4915180/101‑1 (MQ=32)

        tgctgcTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTc                                                                                              <  1:13181247/101‑1 (MQ=32)

         gctgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAGAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCa                                                                                             <  1:11455975/101‑1 (MQ=17)

          ctgctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTAAc                                                                                            <  1:4101343/101‑1 (MQ=17)

          ctgctCGTGAATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCAc                                                                                            >  1:10130267/1‑101 (MQ=17)

            gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                          >  1:14829498/1‑101 (MQ=32)

            gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                          <  1:4904307/101‑1 (MQ=32)

            gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                          >  1:1225245/1‑101 (MQ=32)

            gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                          >  1:3697553/1‑101 (MQ=32)

            gctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGCTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTa                                                                                          <  1:11858125/101‑1 (MQ=17)

             ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                         >  1:8861442/1‑101 (MQ=32)

             ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                         <  1:10242226/101‑1 (MQ=32)

             ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                         <  1:10234179/101‑1 (MQ=32)

             ctCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTAc                                                                                         <  1:150981/101‑1 (MQ=32)

              tCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACg                                                                                        <  1:1361812/101‑1 (MQ=32)

               cGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGc                                                                                       >  1:11291474/1‑101 (MQ=32)

                gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                                      >  1:15194421/1‑101 (MQ=31)

                gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                                      >  1:14010239/1‑101 (MQ=31)

                gTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGca                                                                                      >  1:5277374/1‑101 (MQ=31)

                 tGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGcac                                                                                     >  1:8683457/1‑101 (MQ=11)

                 tGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGcac                                                                                     <  1:5346277/101‑1 (MQ=11)

                                                                         acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                             <  1:8969986/101‑1 (MQ=11)

                                                                         acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                             >  1:1183981/1‑101 (MQ=11)

                                                                         acacTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTgctgct                             <  1:5916857/101‑1 (MQ=11)

                                                                          cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                            <  1:3284404/94‑1 (MQ=2)

                                                                          cgtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTc                            <  1:10044258/94‑1 (MQ=2)

                                                                           gtctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCa                           <  1:5636874/95‑1 (MQ=2)

                                                                           acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                           >  1:5219838/1‑101 (MQ=32)

                                                                           acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                           <  1:5905181/101‑1 (MQ=32)

                                                                           acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                           >  1:10121127/1‑101 (MQ=32)

                                                                           acTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCa                           >  1:210590/1‑101 (MQ=32)

                                                                            tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                          <  1:7199432/96‑1 (MQ=2)

                                                                            tctcgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAg                          <  1:3633632/96‑1 (MQ=2)

                                                                             ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                         <  1:9747902/101‑1 (MQ=32)

                                                                             ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                         >  1:10896206/1‑101 (MQ=32)

                                                                             ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                         <  1:11780201/101‑1 (MQ=32)

                                                                             ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                         <  1:14000799/101‑1 (MQ=32)

                                                                             ttCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGa                         >  1:6627991/1‑101 (MQ=32)

                                                                              tctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGAt                        <  1:6112270/101‑1 (MQ=32)

                                                                               ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                       <  1:11286160/101‑1 (MQ=32)

                                                                               ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                       <  1:7090720/101‑1 (MQ=32)

                                                                               ctcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATg                       <  1:3902251/101‑1 (MQ=32)

                                                                               cgcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATg                       <  1:6368232/99‑1 (MQ=2)

                                                                                tcACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATgg                      <  1:4860957/101‑1 (MQ=32)

                                                                                gcACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATgg                      >  1:2413898/2‑101 (MQ=2)

                                                                                 cACGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGa                     <  1:14935360/101‑1 (MQ=2)

                                                                                 cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa                     >  1:213703/1‑101 (MQ=32)

                                                                                 cACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGa                     <  1:4649622/101‑1 (MQ=32)

                                                                                  aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc                    <  1:6647398/101‑1 (MQ=2)

                                                                                  aCGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGAc                    <  1:4176891/101‑1 (MQ=2)

                                                                                  aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                    <  1:7557719/101‑1 (MQ=32)

                                                                                  aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                    >  1:8146370/1‑101 (MQ=32)

                                                                                  aCGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGAc                    >  1:4011159/1‑101 (MQ=32)

                                                                                   cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg                   <  1:9330761/101‑1 (MQ=2)

                                                                                   cGTTGAGTACCAGACCGAGAAGCGTCACTACGCTCACGTCGACGCTCCGGGCCACGCCGACTACATCAAGAACATGATCACCGGTGCTGCTCAGATGGACg                   <  1:4351784/101‑1 (MQ=2)

                                                                                   cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                   <  1:1604695/101‑1 (MQ=32)

                                                                                   cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                   <  1:492114/101‑1 (MQ=32)

                                                                                   cGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACg                   >  1:6833539/1‑101 (MQ=32)

                                                                                     ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                 >  1:9242016/1‑101 (MQ=32)

                                                                                     ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                 >  1:12726102/1‑101 (MQ=32)

                                                                                     ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                 <  1:14583839/101‑1 (MQ=32)

                                                                                     ttGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgc                 <  1:419487/101‑1 (MQ=32)

                                                                                      tGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcg                >  1:9836238/1‑101 (MQ=32)

                                                                                       gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc               <  1:5305805/101‑1 (MQ=33)

                                                                                       gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc               <  1:14875435/101‑1 (MQ=33)

                                                                                       gAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGgcgc               <  1:13698253/101‑1 (MQ=33)

                                                                                        aaTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGcgcg              >  1:5617732/1‑101 (MQ=32)

                                                                                         aTACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGa             <  1:9742270/101‑1 (MQ=33)

                                                                                          tACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGAt            <  1:2873872/101‑1 (MQ=33)

                                                                                           aCGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATc           >  1:8941787/1‑101 (MQ=32)

                                                                                            cGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATcc          >  1:7634334/1‑101 (MQ=32)

                                                                                             gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt         >  1:3948897/1‑101 (MQ=32)

                                                                                             gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt         >  1:7505216/1‑101 (MQ=32)

                                                                                             gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt         <  1:11323146/101‑1 (MQ=33)

                                                                                             gACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCt         <  1:6553429/101‑1 (MQ=33)

                                                                                              acacCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg        >  1:9408462/1‑101 (MQ=32)

                                                                                              acacCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTg        >  1:6232983/1‑101 (MQ=32)

                                                                                               cacCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTgg       >  1:6931216/1‑101 (MQ=32)

                                                                                                cccccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGt      <  1:3449316/100‑1 (MQ=18)

                                                                                                acCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGt      >  1:817576/1‑101 (MQ=32)

                                                                                                 ccccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTa     >  1:9563240/1‑101 (MQ=32)

                                                                                                  cccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAg    <  1:14750726/101‑1 (MQ=33)

                                                                                                   ccGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGt   >  1:86558/1‑101 (MQ=32)

                                                                                                    cGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt  <  1:10646092/101‑1 (MQ=33)

                                                                                                    cGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGtt  <  1:6145860/101‑1 (MQ=33)

TACGGCGGTGCTGCTCGTGCATTCGACCAGATCGATAACGCGCCGGAAGAAAAAGCTCGTGGTATCACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACTGCCCGGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCTCAGATGGACGGCGCGATCCTGGTAGTT  >  NC_007779/4169204‑4169404

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.