BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,830,747	T→C	L99P (CTG→CCG)	xthA →	exonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,830,747	0	T	C	97.8%	407.4 / ‑0.7	134	L99P (CTG→CCG)	xthA	exonuclease III
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (72/59); total (73/61)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

Only 100 of 149 total aligned reads displayed.

CAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACTGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGTTTT  >  NC_000913/1830658‑1830836

cAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgct                                                                                        <  1:3977473/93‑1 (MQ=255)

cAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgcg                                                                                        >  1:1605107/1‑92 (MQ=255)

 aaaGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATTCCCTCCccgctg                                                                                       >  2:1424924/1‑93 (MQ=255)

 aaaGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgctg                                                                                       >  2:3933586/1‑93 (MQ=255)

 aaaGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgctg                                                                                       >  1:1318750/1‑93 (MQ=255)

 aaaGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGACGCAGCGGCGGATTATTATGGCGGAAATCCCCTCAccgctg                                                                                       >  1:1983163/1‑93 (MQ=255)

  aaGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTgg                                                                                      >  2:1137135/1‑93 (MQ=255)

    gagaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGt                                                                                    >  1:1865047/1‑93 (MQ=255)

    gagaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGt                                                                                    >  2:3873847/1‑93 (MQ=255)

      gaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTaa                                                                                  >  1:3906998/1‑93 (MQ=255)

      gaCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGGaa                                                                                  >  1:2886317/1‑93 (MQ=255)

       aCGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAAt                                                                                 <  2:230559/93‑1 (MQ=255)

        cGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATg                                                                                <  1:1593158/93‑1 (MQ=255)

           cGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCa                                                                             >  2:418877/1‑93 (MQ=255)

            gATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCAc                                                                            <  1:4234958/93‑1 (MQ=255)

             aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGGGGAAATCCCCCCACCGCTGGGTAATGTCAcc                                                                           >  2:1470309/1‑93 (MQ=255)

             aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCAcc                                                                           <  2:1103151/93‑1 (MQ=255)

             aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGGAATGTCAcc                                                                           >  2:3410361/1‑93 (MQ=255)

             aTTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAAACCCCTCACCGCTGGGTAATGTCAcc                                                                           >  1:4520778/1‑93 (MQ=255)

               tGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGt                                                                         <  1:4360037/93‑1 (MQ=255)

               tGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGt                                                                         >  1:4438686/1‑93 (MQ=255)

               tGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGt                                                                         <  2:2191664/93‑1 (MQ=255)

               tGCCGTGCGTCGCGGCTGTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGt                                                                         <  1:4360038/93‑1 (MQ=255)

                gCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGCCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTg                                                                        <  2:178450/93‑1 (MQ=255)

                 ccGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGa                                                                       <  1:2287775/93‑1 (MQ=255)

                 ccGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGGAATGTCCCCGTGa                                                                       >  2:111833/1‑93 (MQ=255)

                   gTGCGTCGCGGCTTTCCCGGTGACGAGGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCACCTACCGCTGGGTGATGTCACCGTGATc                                                                     >  2:2197125/1‑93 (MQ=255)

                   gTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATc                                                                     >  2:2732153/1‑93 (MQ=255)

                   gTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATc                                                                     >  2:2721547/1‑93 (MQ=255)

                   gTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATc                                                                     >  2:847530/1‑93 (MQ=255)

                      cGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAAc                                                                  >  1:667953/1‑93 (MQ=255)

                      cGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAAc                                                                  >  2:2583539/1‑93 (MQ=255)

                       gTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACg                                                                 >  1:3550318/1‑93 (MQ=255)

                       gTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACg                                                                 <  1:570893/93‑1 (MQ=255)

                       gTCGCGGCTTTCCCGGCGACGACGAAAAGGCGCAGCGGCGGATTATTATGGCGGAAATACCCTTACCGCTGGGTAATGTCACCGTGATCaaag                                                                 >  2:4069472/1‑91 (MQ=255)

                        tCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACgg                                                                >  2:262674/1‑93 (MQ=255)

                        tCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACgg                                                                >  2:869990/1‑93 (MQ=255)

                           cgGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTa                                                             <  2:3041044/93‑1 (MQ=255)

                            ggCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTAc                                                            <  2:4341269/93‑1 (MQ=255)

                             gCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACt                                                           >  1:58372/1‑93 (MQ=255)

                              cTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACtt                                                          <  1:2882569/93‑1 (MQ=255)

                               tttCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTc                                                         <  1:263395/93‑1 (MQ=255)

                               tttCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTc                                                         >  2:478863/1‑93 (MQ=255)

                                ttCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTcc                                                        <  2:943984/93‑1 (MQ=255)

                                 tCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCCCCGCTGGGTAATGTCACCGTGATCAACGGTTACTTccc                                                       <  1:4146254/93‑1 (MQ=255)

                                  cccGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGTAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCg                                                      <  1:665481/93‑1 (MQ=255)

                                    cGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCa                                                    <  2:755504/93‑1 (MQ=255)

                                    cGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCa                                                    >  1:1423167/1‑93 (MQ=255)

                                     ggTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAg                                                   >  1:3630927/1‑93 (MQ=255)

                                      gTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAgg                                                  >  1:1257124/1‑93 (MQ=255)

                                      gTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAgg                                                  <  2:2461580/93‑1 (MQ=255)

                                       ggacggcGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAggg                                                 <  1:101434/92‑1 (MQ=255)

                                        gacgacGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGt                                                <  2:3634986/93‑1 (MQ=255)

                                         acgacgAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTg                                               >  1:4355161/1‑93 (MQ=255)

                                             cgaAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAg                                           >  1:2274915/1‑93 (MQ=255)

                                               aaGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGcc                                         >  1:334632/1‑93 (MQ=255)

                                                agagGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcg                                        >  1:1690177/1‑93 (MQ=255)

                                                 gagGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcgc                                       >  2:1286199/1‑93 (MQ=255)

                                                 gagGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcgc                                       <  2:762565/93‑1 (MQ=255)

                                                 gagGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcgc                                       <  2:3393232/93‑1 (MQ=255)

                                                  agGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCcgcg                                      >  2:1506484/1‑93 (MQ=255)

                                                   ggCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGa                                     >  1:410443/1‑93 (MQ=255)

                                                    tcgcAGCGGGGGATTATTATGTCGGAAATCCCCTCACCGCTGGGGAATGTCACCGTGATCAACGGTTACTTCCCGCAGGTTGAAAGCcgcgcc                                    <  2:3384765/92‑3 (MQ=255)

                                                    gcgcAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGAc                                    >  1:2521773/1‑93 (MQ=255)

                                                     cgcAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGAcc                                   >  1:1045355/1‑93 (MQ=255)

                                                       cAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCAt                                 <  2:426895/93‑1 (MQ=255)

                                                       cAGCGGCGGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCAt                                 <  2:31356/93‑1 (MQ=255)

                                                         gcggcgGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATcc                               <  1:4698905/93‑1 (MQ=255)

                                                         gcggcgGATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATcc                               >  2:286933/1‑93 (MQ=255)

                                                         gcggcgGATTATTAT‑GCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCg                              >  2:2881750/1‑93 (MQ=255)

                                                           ggcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGa                             <  1:905203/93‑1 (MQ=255)

                                                           ggcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGa                             <  1:639344/93‑1 (MQ=255)

                                                            gcggATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGAt                            <  1:2536468/93‑1 (MQ=255)

                                                               gATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaa                         >  2:4275503/1‑93 (MQ=255)

                                                               gATTATTATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaa                         >  2:302696/1‑93 (MQ=255)

                                                                attattATGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATaaaa                        >  1:2965476/1‑93 (MQ=255)

                                                                 ttattaTGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAAt                       >  1:636559/1‑93 (MQ=255)

                                                                   attatGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTc                     >  1:2471582/1‑93 (MQ=255)

                                                                    ttatGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTcc                    <  2:1788865/93‑1 (MQ=255)

                                                                     tatGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTccc                   <  1:2321234/93‑1 (MQ=255)

                                                                      atGGCGGAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCg                  <  1:978462/93‑1 (MQ=255)

                                                                       tGGCGGAAATCCCCTCACTGCTGGGGAATGGCACCGTGATCAACGGTTACTTCCCGCAGGGGGAAAGCCGCGACCATCCGATAAAATTCCCgg                 <  2:4003922/93‑1 (MQ=255)

                                                                            gAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaaa            <  1:188538/93‑1 (MQ=255)

                                                                            gAAATCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCaaaa            <  2:2067159/93‑1 (MQ=255)

                                                                            gAAACCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCAGCCGATAAAATTCCCGGCaaaa            <  2:3651089/93‑1 (MQ=255)

                                                                             aaaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAg           <  2:3202899/93‑1 (MQ=255)

                                                                              aaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgc          >  2:3879362/1‑93 (MQ=255)

                                                                              aaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgc          <  2:2501229/93‑1 (MQ=255)

                                                                              aaTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgc          <  2:658480/93‑1 (MQ=255)

                                                                               aTCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcg         >  1:2014388/1‑93 (MQ=255)

                                                                                tCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcgc        >  1:4726619/1‑93 (MQ=255)

                                                                                tCCCCTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAgcgc        >  2:847254/1‑93 (MQ=255)

                                                                                 ccccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCa       <  2:4419901/93‑1 (MQ=255)

                                                                                 ccccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCa       <  1:3395106/93‑1 (MQ=255)

                                                                                  cccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAg      >  2:3245948/1‑93 (MQ=255)

                                                                                   ccTCACCGCTGGGTAATGTCCCCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGCCCATCCGATAAAATTCCCGGCAAAAGCGCAGt     <  1:171025/93‑1 (MQ=255)

                                                                                   ccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGt     <  1:3183773/93‑1 (MQ=255)

                                                                                   ccTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGt     <  2:2524733/93‑1 (MQ=255)

                                                                                    cTCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGtt    <  2:1877397/93‑1 (MQ=255)

                                                                                     tCACCGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGttt   >  1:2418842/1‑93 (MQ=255)

CAAAGAGACGCCGATTGCCGTGCGTCGCGGCTTTCCCGGTGACGACGAAGAGGCGCAGCGGCGGATTATTATGGCGGAAATCCCCTCACTGCTGGGTAATGTCACCGTGATCAACGGTTACTTCCCGCAGGGTGAAAGCCGCGACCATCCGATAAAATTCCCGGCAAAAGCGCAGTTTT  >  NC_000913/1830658‑1830836

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