BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	716,988	A→G	G426G (GGT→GGC)	potE ←	putrescine:H(+) symporter/putrescine:ornithine antiporter PotE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	716,988	0	A	G	100.0%	1423.6 / NA	335	G426G (GGT→GGC)	potE	putrescine:H(+) symporter/putrescine:ornithine antiporter PotE
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (188/147); total (188/147)

Only 100 of 335 total aligned reads displayed.

ATTTGCCAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAA  >  NC_000913/716890‑717085

 tttGCCAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCa                                                                                                >  2:1329873/1‑101 (MQ=255)

   tGCCAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATa                                                                                              >  1:1329879/1‑101 (MQ=255)

   tGCCAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATa                                                                                              <  1:1329320/101‑1 (MQ=255)

    gCCAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATAc                                                                                             >  2:1329880/1‑101 (MQ=255)

     ccAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACa                                                                                            >  2:1329883/1‑101 (MQ=255)

         caaacaGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAgtgt                                                                                        >  2:1329888/1‑101 (MQ=255)

          aaacaGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTc                                                                                       <  2:1328941/101‑1 (MQ=255)

          aaacaGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTc                                                                                       >  1:1329889/1‑101 (MQ=255)

           aacaGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTcc                                                                                      >  1:1329893/1‑101 (MQ=255)

             caGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTccaa                                                                                    >  2:1329896/1‑101 (MQ=255)

                 aaGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTccaaccaa                                                                                >  2:1329900/1‑101 (MQ=255)

                    ggCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCaagga                                                                             <  2:1329435/101‑1 (MQ=255)

                     gCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCaaggaa                                                                            >  2:1329903/1‑101 (MQ=255)

                        aTCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGt                                                                         >  2:1329907/1‑101 (MQ=255)

                        aTCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGt                                                                         <  2:1329277/101‑1 (MQ=255)

                          cATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCa                                                                       >  2:1329911/1‑101 (MQ=255)

                          cATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCa                                                                       <  1:1329217/101‑1 (MQ=255)

                           aTTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCAc                                                                      <  2:1329344/101‑1 (MQ=255)

                           aTTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCAc                                                                      >  1:1329915/1‑101 (MQ=255)

                            ttCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACg                                                                     >  1:1329918/1‑101 (MQ=255)

                             tCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGa                                                                    <  1:1329530/101‑1 (MQ=255)

                              cAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGAt                                                                   <  2:1329223/101‑1 (MQ=255)

                              cAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGAt                                                                   >  2:1329922/1‑101 (MQ=255)

                              cAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGAt                                                                   <  1:1329329/101‑1 (MQ=255)

                              cAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGAt                                                                   >  1:1329924/1‑101 (MQ=255)

                               aaTCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATg                                                                  <  1:1329418/101‑1 (MQ=255)

                               aaTCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATg                                                                  <  1:1329399/101‑1 (MQ=255)

                                 tCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGa                                                                <  2:1329690/101‑1 (MQ=255)

                                 tCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGa                                                                <  2:1329061/101‑1 (MQ=255)

                                 tCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGa                                                                <  2:1329130/101‑1 (MQ=255)

                                 tCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGa                                                                <  1:1329436/101‑1 (MQ=255)

                                   gCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATGCAGTGTCCAACCAAGGAAGGTCACGATGGAAc                                                              <  1:1329470/101‑1 (MQ=255)

                                   gCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAAc                                                              <  1:1329378/101‑1 (MQ=255)

                                   gCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAAc                                                              >  2:1329929/1‑101 (MQ=255)

                                    cccTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAAcc                                                             <  1:1329361/101‑1 (MQ=255)

                                     ccTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCg                                                            <  1:1329279/101‑1 (MQ=255)

                                     ccTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCg                                                            >  2:1329932/1‑101 (MQ=255)

                                      cTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGt                                                           >  2:1329938/1‑101 (MQ=255)

                                       ttttttACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTa                                                          >  1:1329942/1‑101 (MQ=255)

                                       ttttttACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTa                                                          <  2:1329323/101‑1 (MQ=255)

                                         ttttACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACa                                                        >  2:1329945/1‑101 (MQ=255)

                                         ttttACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACa                                                        >  1:1329946/1‑101 (MQ=255)

                                          tttACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCa                                                                       >  1:1329951/1‑85 (MQ=255)

                                           ttACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGc                                                      >  2:1329953/1‑101 (MQ=255)

                                              cTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCAt                                                    <  2:1329310/100‑1 (MQ=255)

                                              cTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATg                                                   <  2:1329485/101‑1 (MQ=255)

                                               tttGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATgg                                                  <  2:1329248/101‑1 (MQ=255)

                                                ttGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGc                                                 >  2:1329960/1‑101 (MQ=255)

                                                ttGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGc                                                 >  1:1329959/1‑101 (MQ=255)

                                                   ttttATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGcttc                                              >  2:1329964/1‑101 (MQ=255)

                                                    tttATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGcttct                                             >  1:1329967/1‑101 (MQ=255)

                                                     ttattaACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGcttctt                                            >  1:1329970/1‑101 (MQ=255)

                                                     ttattaACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGctt                                               >  2:1329974/1‑98 (MQ=255)

                                                      tattaACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGct                                                >  1:1329977/1‑96 (MQ=255)

                                                        ttaACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTccc                                         >  2:1329982/1‑101 (MQ=255)

                                                        ttaACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTccc                                         >  1:1329983/1‑101 (MQ=255)

                                                         taACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGaa                                                               <  2:1329447/78‑1 (MQ=255)

                                                          aaCCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCg                                       <  1:1329373/101‑1 (MQ=255)

                                                           aCCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCgg                                      >  2:1329986/1‑101 (MQ=255)

                                                             cGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGAt                                    <  2:1329586/101‑1 (MQ=255)

                                                             cGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGAt                                    >  2:1329991/1‑101 (MQ=255)

                                                              gtgtTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATg                                   >  1:1329992/1‑101 (MQ=255)

                                                               tgtTTATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGa                                  >  2:1329997/1‑101 (MQ=255)

                                                                  ttATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGGTGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTa                               >  2:1330001/1‑101 (MQ=255)

                                                                   tATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTAc                              <  1:1329589/101‑1 (MQ=255)

                                                                   tATTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTAc                              <  2:1329471/101‑1 (MQ=255)

                                                                    aTTTTTCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACa                             <  1:1329768/101‑1 (MQ=255)

                                                                     tttttCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAg                            >  2:1330004/1‑101 (MQ=255)

                                                                      ttttCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAgc                           <  1:1329147/101‑1 (MQ=255)

                                                                      ttttCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAgc                           <  2:1329522/101‑1 (MQ=255)

                                                                      ttttCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAgc                           <  1:1329507/101‑1 (MQ=255)

                                                                        ttCAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAgcgc                         >  1:1330008/1‑101 (MQ=255)

                                                                          cAGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCAt                       >  1:1330011/1‑101 (MQ=255)

                                                                           aGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAgcgc                         >  1:1330017/1‑98 (MQ=255)

                                                                           aGTTCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATa                      >  2:1330015/1‑101 (MQ=255)

                                                                             ttCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGa                                     >  1:1330021/1‑84 (MQ=255)

                                                                             ttCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaaa                    <  2:1329405/101‑1 (MQ=255)

                                                                              tCAAAGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATaaaa                   >  2:1330023/1‑101 (MQ=255)

                                                                                aaaGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGc                 >  2:1330029/1‑101 (MQ=255)

                                                                                aaaGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGc                 >  1:1330028/1‑101 (MQ=255)

                                                                                 aaGCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCt                <  2:1329561/101‑1 (MQ=255)

                                                                                   gCGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCtat              <  1:1329703/101‑1 (MQ=255)

                                                                                    cGTGGTGATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCtata             >  1:1330036/1‑101 (MQ=255)

                                                                                     gtggtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATAc            >  1:1330039/1‑101 (MQ=255)

                                                                                     gtggtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATAc            >  2:1330041/1‑101 (MQ=255)

                                                                                     gtggtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATAc            <  2:1329506/101‑1 (MQ=255)

                                                                                      ttgtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGAt                                    <  1:1329459/74‑1 (MQ=255)

                                                                                      tggtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACa           <  1:1329513/101‑1 (MQ=255)

                                                                                      tggtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACa           <  2:1329375/101‑1 (MQ=255)

                                                                                      tggtgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACa           >  1:1330043/1‑101 (MQ=255)

                                                                                         tgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcg        >  1:1330049/1‑101 (MQ=255)

                                                                                         tgATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcg        >  2:1330050/1‑101 (MQ=255)

                                                                                          gATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcgc       <  1:1329553/101‑1 (MQ=255)

                                                                                          gATACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcgc       <  1:1329614/101‑1 (MQ=255)

                                                                                           aTACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcgcg      <  1:1329480/101‑1 (MQ=255)

                                                                                           aTACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcgcg      >  1:1330053/1‑101 (MQ=255)

                                                                                           aTACCAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATcgcg      <  2:1329441/101‑1 (MQ=255)

                                                                                              ccAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCca   >  1:1330058/1‑101 (MQ=255)

                                                                                               cAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCcaa  <  1:1328864/101‑1 (MQ=255)

                                                                                               cAGGCCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCcaa  <  2:1329618/101‑1 (MQ=255)

ATTTGCCAACAAACAGAAAGGGCGATCATTCAATCGCCCTTTTTTACTTTGTTTTATTAACCGTGTTTATTTTTCAGTTCAAAGCGTGGTGATACCAGACCATACAGTGTCCAACCAAGGAAGGTCACGATGGAACCGTACAGCATGGCTTCTTCCCCGGATGAGTACAGCGCATAAAAGCTATACATCGCGCCAA  >  NC_000913/716890‑717085

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