BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,041,927	A→G	Y50C (TAC→TGC)	b3861 →	putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,041,927	0	A	G	100.0%	601.9 / NA	182	Y50C (TAC→TGC)	b3861	putative GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (114/68); total (114/68)

Only 100 of 182 total aligned reads displayed.

GCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCG  >  NC_000913/4041837‑4042016

gCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCg                                                                                         >  2:1027350/1‑93 (MQ=255)

 cacaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGa                                                                                        >  1:3438651/1‑93 (MQ=255)

  acaTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGat                                                                                       <  2:1662528/93‑1 (MQ=255)

   caTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGata                                                                                      <  1:1352199/93‑1 (MQ=255)

   caTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGata                                                                                      <  1:2159752/93‑1 (MQ=255)

    aTGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGatat                                                                                     <  2:592836/93‑1 (MQ=255)

     tGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATg                                                                                    <  2:1292175/93‑1 (MQ=255)

        tACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACg                                                                                 <  1:1114088/93‑1 (MQ=255)

         acacAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGt                                                                                <  1:957646/93‑1 (MQ=255)

         acacAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGt                                                                                <  2:592911/93‑1 (MQ=255)

                 aCGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTAAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTAt                                                                        >  2:3358709/1‑93 (MQ=255)

                  cGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATc                                                                       <  2:2259847/93‑1 (MQ=255)

                     aaaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaa                                                                    >  2:500513/1‑93 (MQ=255)

                      aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   >  1:2305220/1‑93 (MQ=255)

                      aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   >  1:4174421/1‑93 (MQ=255)

                      aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   <  2:4065425/93‑1 (MQ=255)

                      aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   >  2:4494069/1‑93 (MQ=255)

                      aaTTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCaaaa                                                                   >  1:144791/1‑93 (MQ=255)

                          cAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAAtttt                                                               <  2:1809525/93‑1 (MQ=255)

                           aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTg                                                              >  1:2491782/1‑93 (MQ=255)

                           aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTg                                                              >  1:4379759/1‑93 (MQ=255)

                           aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTg                                                              >  2:2847561/1‑93 (MQ=255)

                           aGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGACATGACGTTAAGTTATCAAAATTTTg                                                              <  2:1678263/93‑1 (MQ=255)

                            gCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGa                                                             >  1:1060147/1‑93 (MQ=255)

                              cTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACa                                                           <  1:47995/93‑1 (MQ=255)

                               tGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACaa                                                          >  2:26005/1‑93 (MQ=255)

                               tGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACaa                                                          >  1:1462026/1‑93 (MQ=255)

                                 gtgtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAg                                                        <  1:4250199/93‑1 (MQ=255)

                                  tgtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAgg                                                       >  1:2934329/1‑93 (MQ=255)

                                   gtCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGc                                                      >  2:4612355/1‑93 (MQ=255)

                                     cGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtt                                                    >  1:1065888/1‑93 (MQ=255)

                                      gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttt                                                   >  2:1484650/1‑93 (MQ=255)

                                      gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttt                                                   >  2:3598146/1‑93 (MQ=255)

                                      gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttt                                                   >  2:4366012/1‑93 (MQ=255)

                                      gAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCttt                                                   <  1:2455503/93‑1 (MQ=255)

                                       aaaaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttt                                                  >  1:3130098/1‑93 (MQ=255)

                                         aaaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCtttttt                                                <  1:3672270/93‑1 (MQ=255)

                                          aaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTa                                               <  1:1145319/93‑1 (MQ=255)

                                          aaTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTa                                               >  1:51490/1‑93 (MQ=255)

                                           aTTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTaa                                              >  2:1936650/1‑93 (MQ=255)

                                            tttATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAAt                                             >  2:3843263/1‑93 (MQ=255)

                                              tattaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTc                                           >  1:3773718/1‑93 (MQ=255)

                                              tattaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTc                                           >  2:3736935/1‑93 (MQ=255)

                                               attaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTcc                                          <  2:81824/93‑1 (MQ=255)

                                               attaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTcc                                          >  2:757277/1‑93 (MQ=255)

                                               attaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTcc                                          <  1:2733577/93‑1 (MQ=255)

                                               attaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTcc                                          <  1:3246574/93‑1 (MQ=255)

                                                ttaAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTccc                                         <  1:4099533/93‑1 (MQ=255)

                                                 taAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCg                                        <  1:4510690/93‑1 (MQ=255)

                                                  aaaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGt                                       >  1:3593067/1‑93 (MQ=255)

                                                   aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGtt                                      >  1:684069/1‑93 (MQ=255)

                                                   aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGtt                                      >  2:3253422/1‑93 (MQ=255)

                                                   aaGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGtt                                      <  2:4255226/93‑1 (MQ=255)

                                                    aGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGttt                                     >  2:1336570/1‑93 (MQ=255)

                                                    aGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGttt                                     <  1:2541244/93‑1 (MQ=255)

                                                        tttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCTCGTTTTcaa                                 >  2:2656770/1‑93 (MQ=255)

                                                         ttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaa                                >  1:2600526/1‑93 (MQ=255)

                                                         ttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaa                                >  2:4453973/1‑93 (MQ=255)

                                                         ttAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATAT‑ACGGTAAGGTATCAAAATTTTTACAAAGGCTTTTTTTATTTCCGGTTTtaaat                               >  2:4742390/1‑93 (MQ=255)

                                                          tAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaat                               >  1:238353/1‑93 (MQ=255)

                                                          tAACGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaat                               <  2:2237192/93‑1 (MQ=255)

                                                           aaCGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatg                              >  1:2237976/1‑93 (MQ=255)

                                                            aCGATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgc                             >  1:2483144/1‑93 (MQ=255)

                                                              gATGCTAAAAAGAAAGGTGAACATGCCTGGGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaa                           >  2:2905161/1‑93 (MQ=255)

                                                              gATGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaa                           >  2:4666495/1‑93 (MQ=255)

                                                               aTGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaa                          >  1:534940/1‑93 (MQ=255)

                                                               aTGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaa                          <  2:2670786/93‑1 (MQ=255)

                                                                tGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaat                         >  1:3311627/1‑93 (MQ=255)

                                                                tGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaat                         <  1:777915/93‑1 (MQ=255)

                                                                tGCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaat                         <  2:3926879/93‑1 (MQ=255)

                                                                 gCTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTcaaatgcaaatg                        >  1:897953/1‑93 (MQ=255)

                                                                  cTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGa                       <  1:2358935/93‑1 (MQ=255)

                                                                  cTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGa                       <  2:1605987/93‑1 (MQ=255)

                                                                  cTAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGa                       >  1:1795768/1‑93 (MQ=255)

                                                                   tAAAAAGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGAc                      <  1:2876694/93‑1 (MQ=255)

                                                                    aaaaaGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACa                     <  2:2994456/93‑1 (MQ=255)

                                                                     aaaaGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAt                    >  2:2318675/1‑93 (MQ=255)

                                                                     aaaaGAAAGGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAt                    >  1:3874749/1‑93 (MQ=255)

                                                                      aaagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAtt                   >  2:1497912/1‑93 (MQ=255)

                                                                      aaagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAtt                   <  2:697846/93‑1 (MQ=255)

                                                                      aaagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACAtt                   <  1:2744912/93‑1 (MQ=255)

                                                                       aagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTc                  >  2:2807769/1‑93 (MQ=255)

                                                                       aagaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTc                  >  1:612962/1‑93 (MQ=255)

                                                                        agaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCTAATGCAAATGACATTCg                 >  1:2444067/1‑93 (MQ=255)

                                                                        agaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCg                 >  1:3286303/1‑93 (MQ=255)

                                                                         gaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                >  2:3169886/1‑93 (MQ=255)

                                                                         gaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                >  2:785222/1‑93 (MQ=255)

                                                                         gaaagGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGa                >  1:182328/1‑93 (MQ=255)

                                                                            agGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATaa             <  2:2103081/93‑1 (MQ=255)

                                                                            agGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATaa             <  1:3564053/93‑1 (MQ=255)

                                                                            agGTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATaa             <  2:1144004/93‑1 (MQ=255)

                                                                             ggTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAAc            <  1:573728/93‑1 (MQ=255)

                                                                              gTGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACg           <  2:187836/93‑1 (MQ=255)

                                                                               tGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACgg          >  1:3854352/1‑93 (MQ=255)

                                                                               tGAACATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACgg          >  1:840184/1‑93 (MQ=255)

                                                                                   cATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa      <  1:3621679/93‑1 (MQ=255)

                                                                                   cATGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCa      >  1:2948723/1‑93 (MQ=255)

                                                                                    aTGCCTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCAc     >  2:4065349/1‑93 (MQ=255)

                                                                                       ccTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTTCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCg  >  2:4426287/1‑93 (MQ=255)

                                                                                       ccTGCGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCg  >  1:3977356/1‑93 (MQ=255)

GCACATGGTACACAGGTACGCAAATTCAGCCTGGTGTCGAAAAATTTATTAAAGATTTTAACGATGCTAAAAAGAAAGGTGAACATGCCTACGATATGACGTTAAGTTATCAAAATTTTGACAAAGGCTTTTTTAATTCCCGTTTTCAAATGCAAATGACATTCGATAACGGTGCACCCG  >  NC_000913/4041837‑4042016

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