BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,370,778	A→G	V145A (GTA→GCA)	nanA ←	N‑acetylneuraminate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,370,778	0	A	G	100.0%	747.3 / NA	232	V145A (GTA→GCA)	nanA	N‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (134/98); total (134/98)

Only 100 of 234 total aligned reads displayed.

TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370688‑3370867

    tCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCAc                                                                                     >  1:304045/1‑93 (MQ=255)

      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  2:95072/1‑93 (MQ=255)

      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  2:5088372/1‑93 (MQ=255)

      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   <  1:4373600/93‑1 (MQ=255)

       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  <  2:4810590/93‑1 (MQ=255)

       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  2:3579414/1‑93 (MQ=255)

       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  1:737427/1‑93 (MQ=255)

       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  1:4578910/1‑93 (MQ=255)

       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  1:251116/1‑93 (MQ=255)

         tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                <  1:6172664/93‑1 (MQ=255)

         tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGGGTTGATCTGATGCAGGGTCAGTTTTGCCCCACTCAgg                                                                                >  2:1953480/1‑93 (MQ=255)

            agaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  1:3128748/1‑88 (MQ=255)

            agaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  <  2:3128748/88‑1 (MQ=255)

             gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            >  1:3311303/1‑93 (MQ=255)

             gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            <  2:4053422/93‑1 (MQ=255)

             gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            >  1:2777148/1‑93 (MQ=255)

              aTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTgg                                                                           >  2:1638494/1‑93 (MQ=255)

               tCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGa                                                                          >  1:2254727/1‑93 (MQ=255)

                 gCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAt                                                                        <  1:1634009/93‑1 (MQ=255)

                  ccAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtg                                                                       >  1:337494/1‑93 (MQ=255)

                  ccAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtg                                                                       >  2:837789/1‑93 (MQ=255)

                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  1:4269810/1‑93 (MQ=255)

                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  1:289158/1‑93 (MQ=255)

                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  2:6096474/1‑93 (MQ=255)

                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    <  2:4713770/93‑1 (MQ=255)

                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    <  1:102456/93‑1 (MQ=255)

                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    >  1:201613/1‑93 (MQ=255)

                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    >  1:1709996/1‑93 (MQ=255)

                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    >  1:368954/1‑93 (MQ=255)

                      agGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttgt                                                                   <  1:5295794/93‑1 (MQ=255)

                       ggTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTa                                                                  <  1:531839/93‑1 (MQ=255)

                       ggTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTa                                                                  >  1:4312269/1‑93 (MQ=255)

                        gTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTac                                                                 >  1:1117952/1‑93 (MQ=255)

                          cTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTacac                                                               >  2:2171070/1‑93 (MQ=255)

                          cTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTacac                                                               >  2:1356712/1‑93 (MQ=255)

                           tGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcacc                                                              <  2:2925351/93‑1 (MQ=255)

                            gTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcacca                                                             >  1:1970371/1‑93 (MQ=255)

                              ttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACaccacc                                                           <  1:1542882/93‑1 (MQ=255)

                               tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                          <  1:4645867/93‑1 (MQ=255)

                               tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                          >  1:4714756/1‑93 (MQ=255)

                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        <  1:3621575/93‑1 (MQ=255)

                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        <  1:4623335/93‑1 (MQ=255)

                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        <  2:1831575/93‑1 (MQ=255)

                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        <  1:2350070/93‑1 (MQ=255)

                                  gcgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCg                                                       <  2:2686136/93‑1 (MQ=255)

                                   cgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCgg                                                      <  1:1839971/93‑1 (MQ=255)

                                    gcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGc                                                     >  2:3901117/1‑93 (MQ=255)

                                      aCCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaa                                                   >  2:3717966/1‑93 (MQ=255)

                                      aCCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaa                                                   >  2:4431315/1‑93 (MQ=255)

                                      aCCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaa                                                   >  1:389148/1‑93 (MQ=255)

                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                  >  2:4668349/1‑93 (MQ=255)

                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                  <  1:5375198/93‑1 (MQ=255)

                                        cTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAc                                                 >  2:3406770/1‑93 (MQ=255)

                                         tACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                                <  2:5877925/93‑1 (MQ=255)

                                           cGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTCTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCAt                                              >  2:464966/1‑93 (MQ=255)

                                            gCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATc                                             <  2:5189477/93‑1 (MQ=255)

                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCAtac                                            >  2:5347571/1‑91 (MQ=255)

                                               aGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgc                                          >  2:4881994/1‑93 (MQ=255)

                                               aGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgc                                          <  1:3799176/93‑1 (MQ=255)

                                                 gCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgccg                                        >  2:434101/1‑93 (MQ=255)

                                                   attGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  2:1139424/3‑93 (MQ=255)

                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  1:341104/1‑93 (MQ=255)

                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      <  2:5391260/93‑1 (MQ=255)

                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  1:5080460/1‑93 (MQ=255)

                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    >  2:6066305/1‑93 (MQ=255)

                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    >  1:4851136/1‑93 (MQ=255)

                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  <  1:2063739/93‑1 (MQ=255)

                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  <  1:884158/93‑1 (MQ=255)

                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  >  2:6012438/1‑93 (MQ=255)

                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  >  2:1426431/1‑93 (MQ=255)

                                                        aaCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                 <  1:2632456/93‑1 (MQ=255)

                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                >  1:1644555/1‑93 (MQ=255)

                                                           aaGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                              <  2:2025304/93‑1 (MQ=255)

                                                             gtgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTg                            >  1:3944614/1‑93 (MQ=255)

                                                               gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                          <  1:399119/93‑1 (MQ=255)

                                                                 gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  2:2386847/1‑93 (MQ=255)

                                                                 gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  1:2808778/1‑93 (MQ=255)

                                                                   tgatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                      <  2:3332446/93‑1 (MQ=255)

                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  2:2034132/93‑1 (MQ=255)

                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     >  1:408961/1‑93 (MQ=255)

                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  1:2444061/93‑1 (MQ=255)

                                                                     atctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAg                    >  1:1169958/1‑93 (MQ=255)

                                                                      tctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGt                   >  1:1951784/1‑93 (MQ=255)

                                                                      tctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGt                   >  2:5619462/1‑93 (MQ=255)

                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                 <  1:3585800/93‑1 (MQ=255)

                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                 <  2:285742/93‑1 (MQ=255)

                                                                         gatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGAt                >  1:3028093/1‑93 (MQ=255)

                                                                         gatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGAt                <  1:5716073/93‑1 (MQ=255)

                                                                          atcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATc               >  2:3221535/1‑93 (MQ=255)

                                                                             cAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa            <  2:5128832/93‑1 (MQ=255)

                                                                             cAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa            <  2:3530152/93‑1 (MQ=255)

                                                                             cAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa            >  1:58263/1‑93 (MQ=255)

                                                                              aGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAg           >  2:360107/1‑93 (MQ=255)

                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  2:4895327/1‑93 (MQ=255)

                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  2:5478098/1‑93 (MQ=255)

                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  2:5774461/1‑93 (MQ=255)

                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  1:1760885/1‑93 (MQ=255)

                                                                                 gTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtgt        <  2:4569265/93‑1 (MQ=255)

                                                                                       tttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCa  <  2:2610814/93‑1 (MQ=255)

                                                                                       tttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCa  <  1:3315839/93‑1 (MQ=255)

TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370688‑3370867

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