BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,011,893	G→A	A123V (GCG→GTG)	b2875 ←	putative synthases

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,011,893	0	G	A	99.5%	766.2 / NA	216	A123V (GCG→GTG)	b2875	putative synthases
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (97/118); minor base C (0/1); total (97/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

GCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCGCCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGAC  >  NC_000913/3011803‑3011982

 cATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcaccg                                                                                        >  2:438550/1‑93 (MQ=255)

  aTGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGcaccgc                                                                                       >  2:2103773/1‑93 (MQ=255)

   tGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCa                                                                                      <  2:4451195/93‑1 (MQ=255)

    ggAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCAc                                                                                     <  2:5706189/93‑1 (MQ=255)

    ggAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCAc                                                                                     <  1:4885173/93‑1 (MQ=255)

     gagaAGCGTGGTTGGTGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCAct                                                                                    >  1:1277015/1‑93 (MQ=255)

      agaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActc                                                                                   <  2:2972109/93‑1 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  <  1:718740/93‑1 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  >  1:1036553/1‑93 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  >  2:2128889/1‑93 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  >  2:4364744/1‑93 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  <  2:4908538/93‑1 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  <  2:5752584/93‑1 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  >  1:1479051/1‑93 (MQ=255)

       gaAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCActct                                                                                  <  2:864521/93‑1 (MQ=255)

         aGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTgg                                                                                <  2:3583504/93‑1 (MQ=255)

          gCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTggg                                                                               <  2:3992175/93‑1 (MQ=255)

          gCGTGGTTGATGGTGGGAATAGTTCACGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTggg                                                                               >  1:5249660/1‑93 (MQ=255)

             tGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTa                                                                            <  1:5815318/93‑1 (MQ=255)

             tGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTa                                                                            <  2:1022815/93‑1 (MQ=255)

             tGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTa                                                                            >  1:3727592/1‑93 (MQ=255)

                 tGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGc                                                                        <  1:2553910/93‑1 (MQ=255)

                 tGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGc                                                                        <  1:1899159/93‑1 (MQ=255)

                  gATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGcc                                                                       <  2:1246790/93‑1 (MQ=255)

                  gATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATAAAATCATAATAGCACGGCACTCTGGGCTATGGcc                                                                       >  2:4224471/1‑93 (MQ=255)

                   aTGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGccc                                                                      <  2:3285567/93‑1 (MQ=255)

                   aTGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGccc                                                                      <  1:1786554/93‑1 (MQ=255)

                     ggtggGAATAGTTCAGGATTGAGGCTCTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCgg                                                                    <  2:5094000/93‑1 (MQ=255)

                      gtggGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGt                                                                   <  1:4284731/93‑1 (MQ=255)

                       tggGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGtt                                                                  >  1:5052891/1‑93 (MQ=255)

                        gggAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTg                                                                 <  1:3912914/93‑1 (MQ=255)

                           aaTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACa                                                              <  2:25008/93‑1 (MQ=255)

                           aaTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACa                                                              >  1:2879427/1‑93 (MQ=255)

                           aaTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACa                                                              <  1:755544/93‑1 (MQ=255)

                            aTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACAt                                                             <  1:2106428/93‑1 (MQ=255)

                            aTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACAt                                                             >  2:2958960/1‑93 (MQ=255)

                            aTAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACAt                                                             <  2:4833336/93‑1 (MQ=255)

                             tAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATg                                                            >  2:3381904/1‑93 (MQ=255)

                             tAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATg                                                            >  1:4544009/1‑93 (MQ=255)

                             tAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATg                                                            >  1:4945826/1‑93 (MQ=255)

                              aGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGc                                                           <  1:4864234/93‑1 (MQ=255)

                               gTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGcc                                                          >  2:1968847/1‑93 (MQ=255)

                                  cAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcg                                                       <  1:5467634/93‑1 (MQ=255)

                                   aGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcgc                                                      <  1:3583740/93‑1 (MQ=255)

                                    ggATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcgcg                                                     >  2:2156797/1‑93 (MQ=255)

                                    ggATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcgcg                                                     <  2:742344/93‑1 (MQ=255)

                                    agATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCcgcg                                                     <  2:2578642/92‑1 (MQ=255)

                                     gATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCgcgc                                                    >  1:2505311/1‑93 (MQ=255)

                                     gATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCgcgc                                                    <  1:2939166/93‑1 (MQ=255)

                                         gAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAAt                                                >  1:2028596/1‑93 (MQ=255)

                                          aGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATc                                               <  2:3848776/93‑1 (MQ=255)

                                             cTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCCCCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAac                                            <  2:23777/93‑1 (MQ=255)

                                               ttCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAacac                                          >  2:181249/1‑93 (MQ=255)

                                                tCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACAcc                                         <  1:4381166/93‑1 (MQ=255)

                                                tCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACAcc                                         <  1:3305926/93‑1 (MQ=255)

                                                tCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACAcc                                         >  1:5188471/1‑93 (MQ=255)

                                                tCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACAcc                                         <  2:5375457/93‑1 (MQ=255)

                                                 tgcgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCa                                        <  1:4602953/92‑1 (MQ=255)

                                                 cgcgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCa                                        <  2:3250255/93‑1 (MQ=255)

                                                  gcgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAg                                       >  2:6242964/1‑93 (MQ=255)

                                                  gcgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAg                                       <  1:599551/93‑1 (MQ=255)

                                                   cgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGa                                      <  1:4477959/93‑1 (MQ=255)

                                                   cgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGa                                      <  1:649041/93‑1 (MQ=255)

                                                   cgATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGa                                      <  1:4218412/93‑1 (MQ=255)

                                                    gATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGAc                                     >  2:4565390/1‑93 (MQ=255)

                                                    gATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGAc                                     >  1:6121992/1‑93 (MQ=255)

                                                    gATAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGAc                                     <  1:5701365/93‑1 (MQ=255)

                                                     aTAAATATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACg                                    <  2:3158635/93‑1 (MQ=255)

                                                       aaaTATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGGCGTg                                  <  1:3858340/93‑1 (MQ=255)

                                                        aaTATCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTgg                                 <  1:5544241/93‑1 (MQ=255)

                                                         atatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGa                                <  1:3364636/93‑1 (MQ=255)

                                                         atatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGa                                <  1:5970719/93‑1 (MQ=255)

                                                          tatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGAc                               >  2:5826/1‑93 (MQ=255)

                                                          tatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGAc                               >  2:2191511/1‑93 (MQ=255)

                                                          tatCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGAc                               >  1:2815017/1‑93 (MQ=255)

                                                           atCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACg                              <  2:3742952/93‑1 (MQ=255)

                                                            tCAGCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGc                             <  1:5620648/93‑1 (MQ=255)

                                                               gCAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATa                          >  2:5980777/1‑93 (MQ=255)

                                                                cAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAc                         <  1:1184921/93‑1 (MQ=255)

                                                                cAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAc                         <  2:602290/93‑1 (MQ=255)

                                                                cAACGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAc                         <  1:6258656/93‑1 (MQ=255)

                                                                 aaCGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAcc                        <  1:2229158/93‑1 (MQ=255)

                                                                 aaCGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATAcc                        >  2:5584922/1‑93 (MQ=255)

                                                                   cGGCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCAt                      <  1:226364/93‑1 (MQ=255)

                                                                    ggCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATGGCCCGGTTGACATGCCCCGCGCCGATCAACAGCAGACGTGGACGCATACCATg                     >  1:3126160/1‑93 (MQ=255)

                                                                    ggCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATg                     >  2:5648902/1‑93 (MQ=255)

                                                                     gCAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATgg                    <  2:4179436/93‑1 (MQ=255)

                                                                      cAATATCAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGa                   <  2:228809/93‑1 (MQ=255)

                                                                            cAAATCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGa             <  2:4312698/93‑1 (MQ=255)

                                                                             aaaTCCTAATGGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGAt            <  2:3322087/93‑1 (MQ=255)

                                                                             aaaTCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGAt            >  2:4705770/1‑93 (MQ=255)

                                                                             aaaTCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGAt            <  2:1375560/93‑1 (MQ=255)

                                                                              aaTCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATa           >  2:4673031/1‑93 (MQ=255)

                                                                               aTCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATaa          >  2:5507638/1‑93 (MQ=255)

                                                                                tCCTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATaaa         >  1:6045569/1‑93 (MQ=255)

                                                                                  cTAATAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATa       <  1:3227553/93‑1 (MQ=255)

                                                                                    aaTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACt     <  1:3005914/93‑1 (MQ=255)

                                                                                     aTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTg    <  1:1653764/93‑1 (MQ=255)

                                                                                     aTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTg    <  2:2740736/93‑1 (MQ=255)

                                                                                     aTAGCACCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTg    <  2:637705/93‑1 (MQ=255)

GCATGGAGAAGCGTGGTTGATGGTGGGAATAGTTCAGGATTGAGGCTTTCGCGATAAATATCAGCAACGGCAATATCAAATCCTAATAGCGCCGCACTCTGGGCTATCGCCCGGTTGACATGCCCCGCGCCAATCAACACCAGACGTGGACGCATACCATGGACGCTGATAAATACTGAC  >  NC_000913/3011803‑3011982

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: `‑`