BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,074,121	(T)_5→6	coding (82/966 nt)	yggF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,074,116	1	.	T	96.2%	529.2 / 11.4	158	coding (87/966 nt)	yggF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/3); new base T (85/66); total (89/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 190 total aligned reads displayed.

TCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑TTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGACATCATT  >  NC_000913/3074024‑3074203

tCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttt‑a                                                                                           >  2:3571376/1‑93 (MQ=255)

 cAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttt‑a‑t                                                                                         <  1:1373104/93‑1 (MQ=255)

  aaTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCCGGCCATCAatttc‑a‑tt                                                                                        >  1:4190554/1‑93 (MQ=255)

  aaTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttt‑a‑tt                                                                                        >  2:651119/1‑93 (MQ=255)

   aTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttt‑a‑ttt                                                                                       >  1:3281629/1‑93 (MQ=255)

   aTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttt‑a‑ttt                                                                                       >  1:291021/1‑93 (MQ=255)

   aTTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAatttt‑a‑ttt                                                                                       <  1:2084575/93‑1 (MQ=255)

      tctcCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑tttttt                                                                                    >  1:2690477/1‑92 (MQ=255)

       ctcCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑ttttttg                                                                                   <  1:4597737/93‑3 (MQ=255)

        tcCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑ttttttgt                                                                                  >  1:1808093/1‑90 (MQ=255)

        tcCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑ttttttgt                                                                                  >  2:1679218/1‑90 (MQ=255)

         ccTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑ttttttgtc                                                                                 <  2:3092744/93‑5 (MQ=255)

           tccGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGc                                                                               <  2:914173/91‑1 (MQ=255)

            tCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGcc                                                                              >  1:2006577/1‑93 (MQ=255)

              gCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCcac                                                                            >  2:1162840/1‑93 (MQ=255)

               ccGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCcaca                                                                           <  2:1054181/93‑1 (MQ=255)

               ccGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCcaca                                                                           >  1:480969/1‑93 (MQ=255)

               ccGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCcaca                                                                           <  2:2002474/93‑1 (MQ=255)

                cGATACCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACAt                                                                          <  1:4669077/93‑1 (MQ=255)

                cGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACAt                                                                          <  1:3783170/93‑1 (MQ=255)

                cGATAACCACCCGGCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCAGCGCGGTGACCGCCAGCCCATCAATTTT‑ATTTTTT‑GTCGCCACAt                                                                          <  1:3783171/93‑1 (MQ=255)

                 gATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATc                                                                         >  2:401988/1‑93 (MQ=255)

                 gATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATc                                                                         <  2:3460925/93‑1 (MQ=255)

                  aTAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATcc                                                                        <  2:2358395/93‑1 (MQ=255)

                  aTAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑TTTTT‑GTCGCCACATCCg                                                                       <  1:3092089/93‑1 (MQ=255)

                    aaCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCgg                                                                      >  1:672001/1‑93 (MQ=255)

                     aCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGt                                                                     >  1:1226166/1‑93 (MQ=255)

                     aCCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGt                                                                     >  1:3574995/1‑93 (MQ=255)

                      ccACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGtt                                                                    <  1:2841355/93‑1 (MQ=255)

                      ccACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGtt                                                                    <  2:180302/93‑1 (MQ=255)

                       cACTCGTCCATGTAAAGCGTCGTTGTTTAATGCTTGTCGCTTCGTTGTGACCTCTAGGCCATCAATTTT‑ATTTTTT‑GTCGCCATATCCGGttt                                                                   <  2:2014455/93‑1 (MQ=255)

                       cACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTTGT‑TTTTT‑GTCGCCACATCCGGttt                                                                   >  2:3244102/1‑93 (MQ=255)

                        aCTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTg                                                                  >  1:3096983/1‑93 (MQ=255)

                         cTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGTCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGc                                                                 <  2:4190517/93‑1 (MQ=255)

                         cTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGc                                                                 <  1:854097/93‑1 (MQ=255)

                         cTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGc                                                                 >  1:3206985/1‑93 (MQ=255)

                          tCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCg                                                                >  1:3903869/1‑93 (MQ=255)

                          tCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCg                                                                >  2:4353128/1‑93 (MQ=255)

                            gTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGc                                                              <  2:694549/93‑1 (MQ=255)

                             tCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGcc                                                             >  1:1688477/1‑93 (MQ=255)

                              ccACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCa                                                            >  1:2122820/1‑93 (MQ=255)

                              ccACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCa                                                            <  2:2293569/93‑1 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAg                                                           >  2:2121336/1‑93 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAg                                                           <  1:2430733/93‑1 (MQ=255)

                               cACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAg                                                           <  2:883232/93‑1 (MQ=255)

                                aCGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAgg                                                          <  2:510159/93‑1 (MQ=255)

                                 cGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGc                                                         >  1:4703838/1‑93 (MQ=255)

                                 cGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGc                                                         <  2:2004061/93‑1 (MQ=255)

                                   aaaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAg                                                       >  1:760219/1‑93 (MQ=255)

                                    aaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCCTCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCCGc                                                      >  2:1961361/1‑93 (MQ=255)

                                    aaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGGAGc                                                      >  2:56263/1‑93 (MQ=255)

                                    aaaGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGc                                                      >  2:2930264/1‑93 (MQ=255)

                                    aaaGCGACGTCGTTTAATGCCTGGCGCATAGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGc                                                      >  2:4427595/1‑93 (MQ=255)

                                    aaaGCGACGACGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGc                                                      >  2:1243951/1‑93 (MQ=255)

                                      aGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCa                                                    >  1:4374992/1‑93 (MQ=255)

                                      aGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCa                                                    <  2:1411343/93‑1 (MQ=255)

                                      aGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCa                                                    >  2:4504966/1‑93 (MQ=255)

                                       gCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCaa                                                   <  2:578549/93‑1 (MQ=255)

                                       gCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCaa                                                   >  2:2174718/1‑93 (MQ=255)

                                           cgtcgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCt                                               >  1:3426981/1‑93 (MQ=255)

                                             tcgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGc                                             <  2:596881/93‑1 (MQ=255)

                                             tcgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGc                                             <  1:814183/93‑1 (MQ=255)

                                              cgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGcc                                            <  2:738430/93‑1 (MQ=255)

                                              cgtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGcc                                            <  2:175686/93‑1 (MQ=255)

                                               gtTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCt                                           >  2:2431569/1‑93 (MQ=255)

                                                     tGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGTCACATCCGGTTTGCGGCCAGGCAGCTAACGCTGCCCGTTCCg                                     >  2:4217636/1‑93 (MQ=255)

                                                        cTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTaa                                  >  1:3257197/1‑93 (MQ=255)

                                                         tGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAAc                                 <  1:93783/93‑1 (MQ=255)

                                                             gcATCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGa                             >  1:3382845/1‑93 (MQ=255)

                                                               aTCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGaaa                           >  1:2905380/1‑93 (MQ=255)

                                                               aTCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGaaa                           >  2:2826717/1‑93 (MQ=255)

                                                               aTCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGaaa                           <  2:3296922/93‑1 (MQ=255)

                                                                tCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaat                          >  1:3822799/1‑93 (MQ=255)

                                                                tCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaat                          >  1:413065/1‑93 (MQ=255)

                                                                tCGCGGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaat                          >  2:2981920/1‑93 (MQ=255)

                                                                tCGCGGTGACCGCCAGGCCATCAATTTT‑A‑TTTTTCGTCGCCACATCCGGTTTGCGGCCATGCAGCCAACGCTGCCTGGTCCGTCACGCGCaat                          >  2:1206660/1‑93 (MQ=255)

                                                                  gcgGTAACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaataa                        >  2:4709194/1‑93 (MQ=255)

                                                                   cgGTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAaataat                       <  1:100718/93‑1 (MQ=255)

                                                                    ggTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATg                      <  1:1215476/93‑1 (MQ=255)

                                                                    ggTGACCGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATg                      <  2:4140139/93‑1 (MQ=255)

                                                                          cGCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACg                <  1:2232924/93‑1 (MQ=255)

                                                                           gCCAGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGc               <  2:729348/93‑1 (MQ=255)

                                                                              aGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCCAAGCTGCCTGTTCCGTAACGCGCACTCATGGCCACGCCAg            >  1:4073270/1‑93 (MQ=255)

                                                                              aGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAg            >  2:2037832/1‑93 (MQ=255)

                                                                              aGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAg            <  2:2738504/93‑1 (MQ=255)

                                                                              aGGCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAg            >  2:1903335/1‑93 (MQ=255)

                                                                              aGGCCATCAATTTT‑A‑TTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAgg           >  1:1554910/1‑93 (MQ=255)

                                                                               ggCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAgg           >  2:3736393/1‑93 (MQ=255)

                                                                               ggCCATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAgg           >  2:3910616/1‑93 (MQ=255)

                                                                                gCCATCATTTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAggg          <  1:2575694/93‑1 (MQ=255)

                                                                                  cATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAc        >  1:2175511/1‑93 (MQ=255)

                                                                                  cATCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGGTCCGTAACGCGAAATAATGGCCACGACAggggc        >  1:4222362/1‑91 (MQ=255)

                                                                                   aTCAATTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAca       >  2:850760/1‑93 (MQ=255)

                                                                                      aaTTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatca    <  2:3669857/93‑1 (MQ=255)

                                                                                      aaTTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatca    <  1:2057412/93‑1 (MQ=255)

                                                                                      aaTTTT‑ATTTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGACACGCCAGGGAcatca    >  1:1955674/1‑93 (MQ=255)

                                                                                       atttt‑attttTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatcat   <  2:1698604/93‑1 (MQ=255)

                                                                                       atttt‑attttTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatcat   <  1:2962685/93‑1 (MQ=255)

                                                                                       atttt‑attttTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGAcatcat   >  1:2355546/1‑93 (MQ=255)

                                                                                        tttt‑attttTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGTGACATCAtt  <  2:3835997/93‑1 (MQ=255)

TCAATTTCTCCTTCGCCGATAACCACTCGTCCACGAAAAGCGACGTCGTTTAATGCCTGGCGCATCGCGGTGACCGCCAGGCCATCAATTTT‑A‑TTTTT‑GTCGCCACATCCGGTTTGCGGCCAGGCAGCCAACGCTGCCTGTTCCGTAACGCGAAATAATGGCCACGCCAGGGACATCATT  >  NC_000913/3074024‑3074203

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