BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,824,873	G→A	L679L (CTG→CTA)	recG →	DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,824,873	0	G	A	87.9%	824.9 / 75.4	272	L679L (CTG→CTA)	recG	DNA helicase
Reads supporting (aligned to +/- strand): ref base G (15/18); new base A (126/113); total (141/131)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 278 total aligned reads displayed.

GGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCTTC  >  NC_000913/3824783‑3824962

 gATTTACTGCGCGATCAGGCGATGATCCGGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATa                                                                                        <  2:3363641/93‑1 (MQ=255)

 gATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATa                                                                                        >  1:4244940/1‑93 (MQ=255)

 gATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATa                                                                                        <  1:3520525/93‑1 (MQ=255)

  aTTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAg                                                                                       <  1:1825413/93‑1 (MQ=255)

  aTTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAg                                                                                       <  1:2273326/93‑1 (MQ=255)

  aTTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAg                                                                                       >  2:2857540/1‑93 (MQ=255)

   tttACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGa                                                                                      >  1:3838116/1‑93 (MQ=255)

   tttACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGa                                                                                      >  2:5023079/1‑93 (MQ=255)

    ttACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGaa                                                                                     <  2:6294052/93‑1 (MQ=255)

     tACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAAc                                                                                    <  1:2511884/93‑1 (MQ=255)

     tACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAAc                                                                                    <  2:2618328/93‑1 (MQ=255)

     tACTGCGCGATCAGGCGATGATCCCGGAAGCTCAGCGCCTGGCACCCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAAc                                                                                    <  1:1734206/93‑1 (MQ=255)

         gcgcGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTg                                                                                <  1:6401537/93‑1 (MQ=255)

          cgcgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTgg                                                                               <  2:4010073/93‑1 (MQ=255)

          cgcgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTgg                                                                               <  2:2730855/93‑1 (MQ=255)

           gcgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGa                                                                              >  1:5795391/1‑93 (MQ=255)

           gcgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGa                                                                              <  2:4182006/93‑1 (MQ=255)

           gcgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGa                                                                              <  2:4371455/93‑1 (MQ=255)

            cgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGAt                                                                             <  2:4491349/93‑1 (MQ=255)

            cgATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGAt                                                                             >  1:4268473/1‑93 (MQ=255)

               tCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGcc                                                                          <  1:3091312/93‑1 (MQ=255)

                cAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCg                                                                         >  2:4563910/1‑93 (MQ=255)

                  ggCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGga                                                                       <  2:3661256/93‑1 (MQ=255)

                  ggCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGga                                                                       >  1:1047236/1‑93 (MQ=255)

                  ggCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGga                                                                       >  1:4569383/1‑93 (MQ=255)

                    cGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGgaga                                                                     <  2:2761114/93‑1 (MQ=255)

                    cGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGgaga                                                                     <  2:2084259/93‑1 (MQ=255)

                    cGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGgaga                                                                     <  1:5030798/93‑1 (MQ=255)

                     gatgatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGAc                                                                    <  2:600521/93‑1 (MQ=255)

                     gatgatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGAc                                                                    >  2:1856091/1‑93 (MQ=255)

                      atgatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACg                                                                   >  2:881790/1‑93 (MQ=255)

                      atgatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACg                                                                   >  1:917415/1‑93 (MQ=255)

                       tgatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACgg                                                                  >  2:1270419/1‑93 (MQ=255)

                        gatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGa                                                                 <  2:5584091/93‑1 (MQ=255)

                        gatCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGa                                                                 >  2:4685617/1‑93 (MQ=255)

                          tCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAAc                                                               <  2:6206088/93‑1 (MQ=255)

                          tCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAAc                                                               <  2:305660/93‑1 (MQ=255)

                          tCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAAc                                                               <  1:1509405/93‑1 (MQ=255)

                            ccGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGt                                                             >  1:2600642/1‑93 (MQ=255)

                            ccGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGt                                                             >  1:3969524/1‑93 (MQ=255)

                              ggAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTa                                                           >  1:1546715/1‑93 (MQ=255)

                              ggAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTa                                                           >  2:4420091/1‑93 (MQ=255)

                                aaGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACt                                                         <  2:4275259/93‑1 (MQ=255)

                                 aGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTc                                                        >  2:4108226/1‑93 (MQ=255)

                                 aGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTc                                                        >  1:1308567/1‑93 (MQ=255)

                                 aGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTc                                                        >  1:4853979/1‑93 (MQ=255)

                                 aGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTc                                                        <  1:321893/93‑1 (MQ=255)

                                     cAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAAt                                                    >  2:6630110/1‑93 (MQ=255)

                                       gcgcCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGc                                                  >  1:5914377/1‑93 (MQ=255)

                                       gcgcCTGGCACGACATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATgg                                                  >  1:1583718/1‑92 (MQ=255)

                                         gcCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGt                                                >  1:6193024/1‑93 (MQ=255)

                                          ccTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTa                                               <  2:248786/93‑1 (MQ=255)

                                            tGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTaaa                                             >  2:3896602/1‑93 (MQ=255)

                                             ggCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTaaaa                                            <  1:5139522/93‑1 (MQ=255)

                                              ccACGCCATATTCACGAACGTAACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACTTTACTCGAATGCGTAAAAg                                           <  1:2927430/92‑1 (MQ=255)

                                               cACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAgg                                          >  2:3636778/1‑93 (MQ=255)

                                               cACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAgg                                          <  2:6474309/93‑1 (MQ=255)

                                                 cGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcg                                        <  2:1109339/93‑1 (MQ=255)

                                                 cGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcg                                        <  2:4915549/93‑1 (MQ=255)

                                                  gCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcgg                                       <  2:1957434/93‑1 (MQ=255)

                                                  gCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcgg                                       >  2:6099578/1‑93 (MQ=255)

                                                   ccATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcggc                                      >  1:2140842/1‑93 (MQ=255)

                                                   ccATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcggc                                      <  2:2925321/93‑1 (MQ=255)

                                                   ccATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAggcggc                                      >  1:5653860/1‑93 (MQ=255)

                                                    cATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCa                                     <  1:894260/93‑1 (MQ=255)

                                                     atatTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAg                                    <  1:208857/93‑1 (MQ=255)

                                                     atatTCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAg                                    <  1:6458919/93‑1 (MQ=255)

                                                       atTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTg                                  <  1:5266246/93‑1 (MQ=255)

                                                        ttCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTga                                 <  2:3713814/93‑1 (MQ=255)

                                                        ttCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTga                                 >  1:3670478/1‑93 (MQ=255)

                                                         tCACGAACGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTgag                                >  2:4232144/1‑93 (MQ=255)

                                                          cACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTgaga                               <  2:5426170/93‑1 (MQ=255)

                                                              aaCGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGAc                           >  1:4177946/1‑93 (MQ=255)

                                                              aaCGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGAc                           >  1:3399424/1‑93 (MQ=255)

                                                               aCGTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGAcc                          <  2:5826024/93‑1 (MQ=255)

                                                                 gTTACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGc                        <  2:3290243/93‑1 (MQ=255)

                                                                   tACCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCa                      >  2:5246185/1‑93 (MQ=255)

                                                                    aCCCACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCAt                     >  2:3826120/1‑93 (MQ=255)

                                                                     cccACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCAtt                    >  2:2977741/1‑93 (MQ=255)

                                                                     cccACAACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCAtt                    <  2:1726530/93‑1 (MQ=255)

                                                                       cacaACAGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTc                  <  1:6646798/93‑1 (MQ=255)

                                                                         caacaGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAg                >  2:1092828/1‑93 (MQ=255)

                                                                          aacaGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAgg               >  2:6561830/1‑93 (MQ=255)

                                                                           acaGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGt              <  2:1404342/93‑1 (MQ=255)

                                                                           acaGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGt              <  1:3637208/93‑1 (MQ=255)

                                                                           acaGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGt              >  1:1673321/1‑93 (MQ=255)

                                                                           acaGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGt              <  1:124108/93‑1 (MQ=255)

                                                                            ccagCAAAATCCCCAATAAGACGGAGGCTGCCCGAGACGGAACGTTAACCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGtt             <  1:2564178/90‑1 (MQ=255)

                                                                             aGGCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTa            >  1:5218533/1‑93 (MQ=255)

                                                                              ggCAAAAGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTAc           >  1:548931/1‑93 (MQ=255)

                                                                                 aaaaGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCt        <  2:3544893/93‑1 (MQ=255)

                                                                                   aaGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTAc      <  2:5269803/93‑1 (MQ=255)

                                                                                   aaGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTAc      >  2:5681333/1‑93 (MQ=255)

                                                                                   aaGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTAc      <  1:5383277/93‑1 (MQ=255)

                                                                                    aGCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTAcc     <  1:2425557/93‑1 (MQ=255)

                                                                                     gCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCt    >  2:2391693/1‑93 (MQ=255)

                                                                                     gCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCt    <  2:1024705/93‑1 (MQ=255)

                                                                                     gCCCTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCt    <  1:3310699/93‑1 (MQ=255)

                                                                                      cccTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCtt   >  2:1479312/1‑93 (MQ=255)

                                                                                       ccTAATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCTTc  >  2:3463454/1‑93 (MQ=255)

GGATTTACTGCGCGATCAGGCGATGATCCCGGAAGTTCAGCGCCTGGCACGCCATATTCACGAACGTTACCCACAACAGGCAAAAGCCCTGATAGAACGCTGGATGCCGGAGACGGAACGTTACTCGAATGCGTAAAAGGCGGCAGTGAGAAGACCGCCATTTCAGGTTACCCTACCTTC  >  NC_000913/3824783‑3824962

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