BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,630,942	T→C	K383K (AAA→AAG)	guaB ←	IMP dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,630,942	0	T	C	99.5%	656.6 / ‑5.1	200	K383K (AAA→AAG)	guaB	IMP dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (108/91); total (109/91)
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.64e-01

Only 100 of 206 total aligned reads displayed.

CACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCC  >  NC_000913/2630852‑2631031

   cAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTa                                                                                      >  1:1092170/1‑93 (MQ=255)

   cAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTa                                                                                      <  2:2970848/93‑1 (MQ=255)

    aGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTaa                                                                                     <  2:2645683/93‑1 (MQ=255)

    aGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTaa                                                                                     <  2:5365705/93‑1 (MQ=255)

    aGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTaa                                                                                     <  1:909744/93‑1 (MQ=255)

     gTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAg                                                                                    <  2:1865751/93‑1 (MQ=255)

     gTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAg                                                                                    <  2:5026370/93‑1 (MQ=255)

      tttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGa                                                                                   >  1:1577902/1‑93 (MQ=255)

      tttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGa                                                                                   >  2:6314056/1‑93 (MQ=255)

      tttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGa                                                                                   <  1:5972428/93‑1 (MQ=255)

       ttGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGaa                                                                                  <  2:6696836/93‑1 (MQ=255)

        tGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAAc                                                                                 <  1:5397918/93‑1 (MQ=255)

        tGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAAc                                                                                 >  2:4543305/1‑93 (MQ=255)

         gTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACg                                                                                <  1:143083/93‑1 (MQ=255)

         gTCGGCAGCGTTAGCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAGCCCATACCACGGTAAGACTTGTAAGAACg                                                                                <  1:5896666/93‑1 (MQ=255)

           cGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGc                                                                              >  1:4081238/1‑93 (MQ=255)

            ggCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGcc                                                                             <  1:4649592/93‑1 (MQ=255)

            ggCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGcc                                                                             >  1:199406/1‑93 (MQ=255)

               aGCGGTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTg                                                                          <  1:3166129/93‑1 (MQ=255)

                 cGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGt                                                                        >  1:1683279/1‑93 (MQ=255)

                 cGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGt                                                                        <  1:2494005/93‑1 (MQ=255)

                  gTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTa                                                                       <  1:1830932/93‑1 (MQ=255)

                   ttATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTag                                                                      <  2:1640590/93‑1 (MQ=255)

                   ttATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTag                                                                      >  2:5111675/1‑93 (MQ=255)

                    tATCGCTCTGGAAATAACGGTCAGAGGCACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTaga                                                                     <  1:3589276/93‑1 (MQ=255)

                     aTCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTagag                                                                    <  2:4899108/93‑1 (MQ=255)

                     aTCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTagag                                                                    >  2:1788447/1‑93 (MQ=255)

                       cGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGtt                                                                  >  2:6166543/1‑93 (MQ=255)

                        gCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTc                                                                 >  1:1290152/1‑93 (MQ=255)

                        gCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTc                                                                 <  1:2207024/93‑1 (MQ=255)

                         ctctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCg                                                                >  2:5946206/1‑93 (MQ=255)

                         ctctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCg                                                                >  1:6168484/1‑93 (MQ=255)

                          tctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGa                                                               >  1:2839166/1‑93 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              <  1:4394798/93‑1 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              >  2:611077/1‑93 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              <  1:6122898/93‑1 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              >  1:5482017/1‑93 (MQ=255)

                           ctGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAt                                                              <  1:2743042/93‑1 (MQ=255)

                            tGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAtt                                                             >  1:6501405/1‑93 (MQ=255)

                             ggAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGAttt                                                            <  1:6690628/93‑1 (MQ=255)

                              gAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTc                                                           >  1:461109/1‑93 (MQ=255)

                               aaaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCa                                                          >  1:3342701/1‑93 (MQ=255)

                                aaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCATGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAc                                                         <  2:4695343/93‑1 (MQ=255)

                                aaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAc                                                         <  2:5160026/93‑1 (MQ=255)

                                aaTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAc                                                         <  2:5803503/93‑1 (MQ=255)

                                 aTAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAcc                                                        >  2:6601364/1‑93 (MQ=255)

                                  tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAccc                                                       <  2:4424505/93‑1 (MQ=255)

                                  tAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAccc                                                       <  2:3375410/93‑1 (MQ=255)

                                    aCGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCgg                                                     >  1:575388/1‑93 (MQ=255)

                                     cGGTCAGAGGAGCCGTTGGACAGCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGga                                                    <  2:3232810/93‑1 (MQ=255)

                                     cGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGga                                                    >  2:1415319/1‑93 (MQ=255)

                                     cGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGga                                                    <  2:5518485/93‑1 (MQ=255)

                                      ggTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGgag                                                   <  1:6014327/93‑1 (MQ=255)

                                       gTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGgaga                                                  >  2:3837844/1‑93 (MQ=255)

                                        tCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAt                                                 >  2:2122773/1‑93 (MQ=255)

                                        tCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAt                                                 <  2:2789261/93‑1 (MQ=255)

                                        tCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAt                                                 <  1:4750021/93‑1 (MQ=255)

                                           gagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                              >  1:4472031/1‑93 (MQ=255)

                                           gagGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGAttct                                              >  1:3479964/1‑93 (MQ=255)

                                              gAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCa                                           <  1:1635249/93‑1 (MQ=255)

                                                aCCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGt                                         >  2:5680596/1‑93 (MQ=255)

                                                  cTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAc                                       >  1:5336746/1‑93 (MQ=255)

                                                   tttGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAcc                                      >  1:775900/1‑93 (MQ=255)

                                                   tttGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTAcc                                      >  1:5254848/1‑93 (MQ=255)

                                                     tGGACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCg                                    <  1:3111016/93‑1 (MQ=255)

                                                      ggACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGc                                   >  1:2380436/1‑93 (MQ=255)

                                                      ggACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGc                                   >  1:426177/1‑93 (MQ=255)

                                                      ggACATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGc                                   >  1:5055788/1‑93 (MQ=255)

                                                        aCATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTGGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCa                                 >  1:1941993/1‑93 (MQ=255)

                                                        aCATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCa                                 <  1:1417447/93‑1 (MQ=255)

                                                         cATCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAg                                >  2:4952987/1‑93 (MQ=255)

                                                           tCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCa                              <  2:5039497/93‑1 (MQ=255)

                                                           tCGCGCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCa                              <  2:850784/93‑1 (MQ=255)

                                                            cgcgCCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCAt                             >  2:1126876/1‑93 (MQ=255)

                                                             gcgcCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATg                            >  2:3949652/1‑93 (MQ=255)

                                                               gcCGAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGa                          <  1:3537497/93‑1 (MQ=255)

                                                               gcCCAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGa                          <  1:5536634/93‑1 (MQ=255)

                                                                  cAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAcc                       <  1:1774870/93‑1 (MQ=255)

                                                                  cAGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGCTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAAcc                       >  2:3561922/1‑93 (MQ=255)

                                                                   aGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCt                      >  1:374675/1‑93 (MQ=255)

                                                                   aGGGAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCt                      <  1:4229713/93‑1 (MQ=255)

                                                                    gggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTa                     >  1:5738424/1‑93 (MQ=255)

                                                                    gggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTa                     >  2:2338924/1‑93 (MQ=255)

                                                                    gggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTa                     >  2:5444627/1‑93 (MQ=255)

                                                                     ggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAc                    <  1:247876/93‑1 (MQ=255)

                                                                     ggAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAc                    >  1:212482/1‑93 (MQ=255)

                                                                      gAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAcc                   >  2:5596849/1‑93 (MQ=255)

                                                                      gAACCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTAcc                   >  2:749740/1‑93 (MQ=255)

                                                                        aCCCATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCAt                 >  2:4235778/1‑93 (MQ=255)

                                                                         cccATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATc                >  2:3062946/1‑93 (MQ=255)

                                                                          ccATACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCa               <  2:4636498/93‑1 (MQ=255)

                                                                            aTACCACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCAcc             <  2:1232697/93‑1 (MQ=255)

                                                                               ccACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACcgcg          >  2:541337/1‑93 (MQ=255)

                                                                                cACGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCgcgc         >  2:2527585/1‑93 (MQ=255)

                                                                                 aCGGTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCt        <  2:6396217/93‑1 (MQ=255)

                                                                                   ggTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTg      >  1:4993429/1‑93 (MQ=255)

                                                                                   ggTAAGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTg      <  2:4307397/93‑1 (MQ=255)

                                                                                       aGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTgagcg  >  2:1729586/1‑92 (MQ=255)

                                                                                       aGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGcc  >  1:2321965/1‑93 (MQ=255)

                                                                                       aGACTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGcc  <  2:2428018/93‑1 (MQ=255)

CACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACCCGGAGATTCTTCAGTACCCGCCAGCATGGAACCTACCATCACCGCGCTTGCGCC  >  NC_000913/2630852‑2631031

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