BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,905,748	A→G	S83P (TCT→CCT)	pstB ←	ATP‑binding component of high‑affinity phosphate‑specific transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,905,748	0	A	G	100.0%	705.0 / NA	224	S83P (TCT→CCT)	pstB	ATP‑binding component of high‑affinity phosphate‑specific transport system
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (109/115); total (109/115)

Only 100 of 237 total aligned reads displayed.

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837

aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCATCAGTGCGATATCCTGAGGGt                                                                                         >  1:3127191/1‑93 (MQ=255)

aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                         >  2:1487773/1‑93 (MQ=255)

 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                        <  2:1507290/93‑1 (MQ=255)

 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                        <  2:1956787/93‑1 (MQ=255)

  gCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTg                                                                                       <  2:3663220/93‑1 (MQ=255)

  gCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTg                                                                                       <  1:3689613/93‑1 (MQ=255)

    gATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGt                                                                                     >  2:2508732/1‑93 (MQ=255)

    gATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGt                                                                                     >  2:511128/1‑93 (MQ=255)

    gATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGt                                                                                     <  1:1067102/93‑1 (MQ=255)

      tgttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCa                                                                                   <  2:372510/93‑1 (MQ=255)

      tgttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCa                                                                                   <  2:2027912/93‑1 (MQ=255)

      tgttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCa                                                                                   >  1:4090151/1‑93 (MQ=255)

        ttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAgg                                                                                 >  1:2888714/1‑93 (MQ=255)

         tgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGa                                                                                <  1:1642910/93‑1 (MQ=255)

          gtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAt                                                                               <  1:1985314/93‑1 (MQ=255)

           tCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtg                                                                              <  1:3806063/93‑1 (MQ=255)

           tCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtg                                                                              >  1:2615582/1‑93 (MQ=255)

             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            <  2:2651348/93‑1 (MQ=255)

              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           >  1:2004859/1‑93 (MQ=255)

              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           <  2:1693750/93‑1 (MQ=255)

              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           >  1:3905719/1‑93 (MQ=255)

              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           <  1:1300331/93‑1 (MQ=255)

              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           >  2:929147/1‑93 (MQ=255)

                gATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         >  1:27989/1‑93 (MQ=255)

                  tGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGc                                                                       >  1:1348543/1‑93 (MQ=255)

                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGcc                                                                      >  2:2623653/1‑93 (MQ=255)

                    gACATCGGAAGCGGCGTCGGTGTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCa                                                                     <  1:1231382/93‑1 (MQ=255)

                    gACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCa                                                                     <  2:2175823/93‑1 (MQ=255)

                    gACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCa                                                                     >  2:3955298/1‑93 (MQ=255)

                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 <  1:267050/93‑1 (MQ=255)

                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 >  2:1095921/1‑93 (MQ=255)

                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 <  1:2306210/93‑1 (MQ=255)

                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                               <  2:3133838/93‑1 (MQ=255)

                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                               >  1:4603624/1‑93 (MQ=255)

                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                               <  2:1942275/93‑1 (MQ=255)

                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                              >  2:1014500/1‑93 (MQ=255)

                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                              >  1:3675519/1‑93 (MQ=255)

                            aaaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAg                                                             >  1:3795542/1‑93 (MQ=255)

                             aaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                            >  2:2417490/1‑93 (MQ=255)

                             aaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                            >  1:1820669/1‑93 (MQ=255)

                             aaCGGCGGCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCACTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                            >  2:3072064/1‑93 (MQ=255)

                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                           <  2:1538709/93‑1 (MQ=255)

                               cGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAt                                                          >  2:3907044/1‑93 (MQ=255)

                                ggCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                                         <  1:2429702/93‑1 (MQ=255)

                                 gCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAttt                                                        >  2:740228/1‑93 (MQ=255)

                                  cGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTc                                                       >  1:3722424/1‑93 (MQ=255)

                                   gTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCa                                                      >  2:4159786/1‑93 (MQ=255)

                                   gTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCa                                                      >  1:3182349/1‑93 (MQ=255)

                                    gcGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     <  1:1898838/92‑1 (MQ=255)

                                     cGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                                    <  2:3249445/93‑1 (MQ=255)

                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   >  2:1421659/1‑93 (MQ=255)

                                        tttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTc                                                 >  2:4303916/1‑93 (MQ=255)

                                          tCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCg                                               <  2:1463014/93‑1 (MQ=255)

                                          tCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCg                                               <  2:4052975/93‑1 (MQ=255)

                                             ggAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCAc                                            <  1:4433865/93‑1 (MQ=255)

                                                acacCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                                         >  1:3089854/1‑93 (MQ=255)

                                                    cATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTc                                     >  2:4663207/1‑93 (MQ=255)

                                                    cATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTc                                     <  1:1597826/93‑1 (MQ=255)

                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                    <  2:1654095/93‑1 (MQ=255)

                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                    <  2:582340/93‑1 (MQ=255)

                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                    >  2:3675812/1‑93 (MQ=255)

                                                       gCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCgg                                  <  2:347695/93‑1 (MQ=255)

                                                         ccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGt                                <  1:1977927/93‑1 (MQ=255)

                                                         ccACTTGCGCACGCAGCAGTGCGATGTCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCGCGCTGCTCCGGGt                                <  2:3186140/93‑1 (MQ=255)

                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               >  1:3257166/1‑93 (MQ=255)

                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGATa                               >  2:4499252/1‑93 (MQ=255)

                                                           aCTTTAGCACGCAGCAGTGCGATATCCTGAGGGTTGGCCAGGATGTTGTCGCCATAAAGCAGAATATCATCTTCCGCACGCTGATCCggggcc                              >  2:138990/1‑90 (MQ=255)

                                                            cTTTCGCACGCATCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                             <  1:4285479/93‑1 (MQ=255)

                                                             tttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAg                            >  2:3422679/1‑93 (MQ=255)

                                                              ttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGt                           <  1:819987/93‑1 (MQ=255)

                                                                cgcacgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTc                         >  1:3407856/1‑93 (MQ=255)

                                                                 gcacgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCa                        <  2:4084055/93‑1 (MQ=255)

                                                                    cgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAAc                     <  1:3347242/93‑1 (MQ=255)

                                                                    cgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAAc                     <  1:4352278/93‑1 (MQ=255)

                                                                     gcagcaGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACa                    >  1:539971/1‑93 (MQ=255)

                                                                      cagcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAt                   >  1:4568154/1‑93 (MQ=255)

                                                                      cagcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAt                   >  2:3553285/1‑93 (MQ=255)

                                                                       agcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtt                  >  2:2953852/1‑93 (MQ=255)

                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 >  1:959518/1‑93 (MQ=255)

                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 >  1:4389487/1‑93 (MQ=255)

                                                                          agTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg               >  2:4425574/1‑93 (MQ=255)

                                                                           gTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgt              >  1:4482628/1‑93 (MQ=255)

                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt             <  1:1205654/93‑1 (MQ=255)

                                                                             gCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgttg            <  1:2113085/93‑1 (MQ=255)

                                                                              cGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGa           <  1:2910309/93‑1 (MQ=255)

                                                                               gATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa          <  2:853914/93‑1 (MQ=255)

                                                                               gATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTCaa          >  2:2264862/1‑93 (MQ=255)

                                                                               gATATCCTGAGGGTTGGTCAAGATGTTGTCGCCCTCAAGCAGAATTTTACCTTTCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa          <  1:2387588/93‑1 (MQ=255)

                                                                                atatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg         <  2:1565529/93‑1 (MQ=255)

                                                                                atatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg         >  2:2559985/1‑93 (MQ=255)

                                                                                atatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg         <  2:2822261/93‑1 (MQ=255)

                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        >  2:1608374/1‑93 (MQ=255)

                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       >  2:707800/1‑93 (MQ=255)

                                                                                   tCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa      >  1:4200251/1‑93 (MQ=255)

                                                                                    ccTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTAc     <  1:3052013/93‑1 (MQ=255)

                                                                                    ccTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTAc     <  2:2309021/93‑1 (MQ=255)

                                                                                     cTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACg    <  1:3095590/93‑1 (MQ=255)

                                                                                      tGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgc   >  2:4123852/1‑93 (MQ=255)

                                                                                       gaggGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgca  <  1:1490602/93‑1 (MQ=255)

                                                                                       gaggGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgca  <  2:3311451/93‑1 (MQ=255)

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837

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