BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,182,806	C→A	16.2%	P86Q (CCA→CAA)	zupT →	zinc transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,806	0	C	A	16.2%	335.5 / 53.4	130	P86Q (CCA→CAA)	zupT	zinc transporter
Reads supporting (aligned to +/- strand): ref base C (53/56); new base A (13/8); total (66/64)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.25e-01

Only 100 of 139 total aligned reads displayed.

TCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCA  >  NC_000913/3182708‑3182905

tCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCca                                                                                                    >  1:1066122/1‑100 (MQ=255)

tCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCca                                                                                                    >  1:934417/1‑100 (MQ=255)

  cTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcaca                                                                                                  >  1:858896/1‑100 (MQ=255)

     ccGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGc                                                                                               <  2:804179/100‑1 (MQ=255)

        cACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                            <  1:1285750/100‑4 (MQ=255)

        cACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                            <  2:381457/100‑4 (MQ=255)

        cACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                            <  1:237341/100‑4 (MQ=255)

         aCTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                            <  2:548897/99‑4 (MQ=255)

           tAGCAGCTGAAGGAATGTCGCATGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCGTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGGTCATc                                                                                          <  2:3199781/99‑1 (MQ=255)

            agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCg                                                                                        >  1:1341242/1‑100 (MQ=255)

            agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCg                                                                                        <  2:3502984/100‑1 (MQ=255)

            agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggtccag                                                                                        <  2:533913/100‑8 (MQ=255)

            agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggtccag                                                                                        <  1:768727/100‑8 (MQ=255)

             gcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATGTTGGCCTGGACCGCATGTGGCAACATGcggtccagg                                                                                       <  2:1660185/100‑9 (MQ=255)

              tagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCa                                                                                      <  2:287109/99‑1 (MQ=255)

              cagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCa                                                                                      <  2:858890/100‑1 (MQ=255)

              cagcTGAAGGAATGTCGCCTGTATTGGGTTACGGAATCTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCa                                                                                      <  2:2461771/100‑1 (MQ=255)

               agcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCaa                                                                                     >  1:287111/1‑99 (MQ=255)

                  tGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAt                                                                                  >  1:1736938/1‑100 (MQ=255)

                   gAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAtt                                                                                 >  1:3177944/1‑100 (MQ=255)

                    aaGGATTGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGTACCGCATGTTGCCACATGCTCATCCGCAGGAttt                                                                                <  2:3220620/100‑1 (MQ=255)

                       gAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGGTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt                                                                             <  1:225112/100‑1 (MQ=255)

                       gAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt                                                                             <  1:2016634/100‑1 (MQ=255)

                        aaTGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATg                                                                            >  2:1710006/1‑100 (MQ=255)

                         aTGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGc                                                                           <  1:3109018/100‑1 (MQ=255)

                             cGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCaaaa                                                                       <  1:2439742/100‑1 (MQ=255)

                              gCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCaaaaa                                                                      >  1:214498/1‑100 (MQ=255)

                              gCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTAGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGaaaaaa                                                                      <  2:1143766/100‑1 (MQ=255)

                               ccTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCCGGATTCAATGcca                                                                         >  2:3003438/1‑94 (MQ=255)

                               ccTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAAt                                                                     >  1:1776317/1‑100 (MQ=255)

                                 tgtgTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATc                                                                    <  2:2618681/99‑1 (MQ=255)

                                  gtgtTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAAt                                                                     <  2:1974280/97‑1 (MQ=255)

                                  gtgtTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCgg                                                                  >  2:1467531/1‑100 (MQ=255)

                                  gtgaggggTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCCTCCGCAGGATTTAATGGAAAAAACg                                                                   >  2:781963/6‑99 (MQ=255)

                                   tgtTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGt                                                                 >  1:3515647/1‑100 (MQ=255)

                                       gggTTATGGAATTTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAg                                                             <  2:438714/100‑1 (MQ=255)

                                       gggTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAg                                                             >  1:1053788/1‑100 (MQ=255)

                                        ggTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGc                                                            >  1:218504/1‑100 (MQ=255)

                                        ggTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGc                                                            >  1:119086/1‑100 (MQ=255)

                                        ggTTATGGAATGTTTATCTTCGGTCTGCTTGGCTAGTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCTGTGCAg                                                             <  2:1738634/99‑1 (MQ=255)

                                         gTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGAAGGATTTAATGCAAAAATCGGTGCAGcc                                                           >  1:1159315/1‑100 (MQ=255)

                                           tATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGt                                                         >  2:393249/1‑100 (MQ=255)

                                           tATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGt                                                         <  1:2965465/100‑1 (MQ=255)

                                           tATGGAACGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGt                                                         >  1:1917042/1‑100 (MQ=255)

                                               gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGTCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTTCAGACGTTGcc                                                     <  2:1917032/100‑1 (MQ=255)

                                               gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGcc                                                     <  1:1406551/100‑1 (MQ=255)

                                               gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGcc                                                     <  1:392469/100‑1 (MQ=255)

                                               gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGcc                                                     >  2:1227560/1‑100 (MQ=255)

                                                  tGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCaaa                                                  >  2:1465398/1‑100 (MQ=255)

                                                    tttATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAAt                                                >  1:2719037/1‑100 (MQ=255)

                                                        tCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGAt                                            >  1:3129404/1‑100 (MQ=255)

                                                         cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           >  2:2744826/1‑100 (MQ=255)

                                                         cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           <  1:774590/100‑1 (MQ=255)

                                                         cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           <  2:1679836/100‑1 (MQ=255)

                                                         cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGGTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           >  2:1942361/1‑100 (MQ=255)

                                                         cTTCGGTCTGCTGGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTTACCCGCAGGTTTTAATGCCCACCTCGGTGCCGCCGGTCCCAAAATCGATc                                           >  2:638967/1‑100 (MQ=255)

                                                          ttCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                           <  2:2690483/99‑1 (MQ=255)

                                                          ttCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCa                                          <  1:2300346/100‑1 (MQ=255)

                                                           tCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCATCCGTTGCCAAAATCGATCaa                                         >  1:2361018/1‑100 (MQ=255)

                                                            cGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAg                                        <  1:1772402/100‑1 (MQ=255)

                                                             ggTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAgc                                       >  1:2968595/1‑100 (MQ=255)

                                                             ggTCTGCTTGGCTATTTTGGCATGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAGGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAgc                                       <  2:2238566/100‑1 (MQ=255)

                                                              gTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAgcg                                      <  1:470046/100‑1 (MQ=255)

                                                                 tGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGcac                                   <  1:893589/100‑1 (MQ=255)

                                                                  gCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGcaca                                  >  1:3406020/1‑100 (MQ=255)

                                                                  gCTTGGCTATTTTGGCCTAGAGCGCATGTTGCCACATGCTCATCCGAAAAAAAAAAAAAAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGca                                    <  2:13956/98‑1 (MQ=255)

                                                                    ttGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGc                                <  2:1899043/100‑1 (MQ=255)

                                                                         tATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTc                           >  2:297497/1‑100 (MQ=255)

                                                                          aTTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCt                          <  1:2545975/100‑1 (MQ=255)

                                                                          aTTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCt                          <  2:2719019/100‑1 (MQ=255)

                                                                              tGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTc                      <  2:1822070/100‑1 (MQ=255)

                                                                               ggCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCa                     >  2:1323686/1‑100 (MQ=255)

                                                                               ctgcggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCa                     >  1:2550550/11‑100 (MQ=38)

                                                                                tgcggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAc                    >  1:431538/10‑100 (MQ=255)

                                                                                gCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAc                    <  2:1776309/100‑1 (MQ=255)

                                                                                gCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAc                    <  2:1901851/100‑1 (MQ=255)

                                                                                gCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAc                    <  2:2802664/100‑1 (MQ=255)

                                                                                gCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAc                    <  1:2713005/100‑1 (MQ=255)

                                                                                 gcggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAct                   >  2:3494203/9‑100 (MQ=255)

                                                                                  cggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActc                  >  1:2405215/8‑100 (MQ=255)

                                                                                  cggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGATGCCAAAATCGATCAAGCGCACAGCCATTCTCCTCActc                  >  1:3087182/8‑100 (MQ=255)

                                                                                  cTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActc                  >  1:277642/1‑100 (MQ=255)

                                                                                   ggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActct                 >  2:2768088/7‑100 (MQ=255)

                                                                                    gatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACtctc                >  2:2220065/6‑100 (MQ=255)

                                                                                    gatcagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAct                   >  2:843201/6‑97 (MQ=255)

                                                                                      tgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              >  2:1595082/4‑100 (MQ=255)

                                                                                      tgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCAGTTGCCAAAATCGATCAAGCGCACAGCCGTTCTCCTCACTCTCg               >  2:3085977/4‑99 (MQ=255)

                                                                                      aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              <  2:2893904/100‑1 (MQ=255)

                                                                                      aCCGCATGTTGCCACATGCTCATCCACAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCgg              <  1:3096544/100‑1 (MQ=255)

                                                                                       ccGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGc             <  1:863571/100‑1 (MQ=255)

                                                                                        cGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAACATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCa            <  2:1882463/100‑1 (MQ=255)

                                                                                        cGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCa            >  1:3391620/1‑100 (MQ=255)

                                                                                         gCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCAt           >  1:1434140/1‑100 (MQ=255)

                                                                                         gCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCAt           <  1:2089014/100‑1 (MQ=255)

                                                                                          cATGTTGCCACATCGTCACACGCAGGCTTTATGGCACAATTCTGCGCCGCCGAGGCCAAAATCGATCAAGCGCACAGCCACCCTGCTCACTCTCGGCATc          <  2:3225262/100‑1 (MQ=255)

                                                                                            tGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAg        >  2:2009688/1‑100 (MQ=255)

                                                                                            tGTTGCCACATGCTCATCCGCAGGAGTTAGTGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCGCACTCTCGGCATCa         <  2:1675123/99‑1 (MQ=255)

                                                                                               tGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCt     >  2:2677529/1‑100 (MQ=255)

                                                                                                gCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCCAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTg    >  2:1977460/1‑100 (MQ=255)

                                                                                                  cacaTGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCa  <  1:2797621/100‑1 (MQ=255)

TCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCA  >  NC_000913/3182708‑3182905

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