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%	1269.3 / NA	402	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 (216/186); total (216/186)

Only 100 of 417 total aligned reads displayed.

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837

aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                         <  1:5223920/93‑1 (MQ=255)

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

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

     aTGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTc                                                                                    >  1:5981436/1‑93 (MQ=255)

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

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

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

         tgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGa                                                                                >  1:850368/1‑93 (MQ=255)

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

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

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

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

            cGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgt                                                                             <  1:6686471/93‑1 (MQ=255)

             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            >  1:2173823/1‑93 (MQ=255)

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

                gATGGTCATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         >  1:484426/1‑93 (MQ=255)

                gATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         <  2:3081756/93‑1 (MQ=255)

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

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

                gATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         <  2:7162536/93‑1 (MQ=255)

                  tGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGc                                                                       <  1:1861328/93‑1 (MQ=255)

                  tGGACATCGGAAACGGCGGCGGGTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGc                                                                       >  2:7743242/1‑93 (MQ=255)

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

                     aCATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCAt                                                                    <  2:6473934/93‑1 (MQ=255)

                      cATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATc                                                                   <  2:7587396/93‑1 (MQ=255)

                      cATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATc                                                                   <  2:5767955/93‑1 (MQ=255)

                       aTCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCa                                                                  >  1:4639296/1‑93 (MQ=255)

                       aTCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCa                                                                  >  1:5605747/1‑93 (MQ=255)

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

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

                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGGCAGGATGTTGTGGCCATTaa                                                                 >  1:3211645/1‑93 (MQ=255)

                         cGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAg                                                                <  2:1036859/93‑1 (MQ=255)

                         cGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAg                                                                <  1:696687/93‑1 (MQ=255)

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

                             aaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                            <  2:7527295/93‑1 (MQ=255)

                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                           >  1:3121189/1‑93 (MQ=255)

                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                           >  1:2225273/1‑93 (MQ=255)

                               cGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGGCAGGATGTTGTCGCCATGAAGCAGAAt                                                          >  2:5499505/1‑93 (MQ=255)

                                    tCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     >  2:1290067/1‑93 (MQ=255)

                                    tCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     <  2:3917133/93‑1 (MQ=255)

                                    tCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     <  2:2354818/93‑1 (MQ=255)

                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   <  1:4714701/93‑1 (MQ=255)

                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   <  1:7441714/93‑1 (MQ=255)

                                         ttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTcc                                                <  2:7427300/93‑1 (MQ=255)

                                            tGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                             >  1:7210526/1‑93 (MQ=255)

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

                                             ggAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCAc                                            >  1:7921266/1‑93 (MQ=255)

                                              gAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACg                                           >  1:6430725/1‑93 (MQ=255)

                                                  acCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgc                                       <  1:7645453/93‑1 (MQ=255)

                                                   ccATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgct                                      >  1:4450354/1‑93 (MQ=255)

                                                    cATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCTCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTc                                     <  2:2020797/93‑1 (MQ=255)

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

                                                      tGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCg                                   >  1:169551/1‑93 (MQ=255)

                                                      tGCCAACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCg                                   <  2:3262283/93‑1 (MQ=255)

                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                 <  2:635741/93‑1 (MQ=255)

                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                 <  1:2354631/93‑1 (MQ=255)

                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGtc                               >  1:1411151/1‑92 (MQ=255)

                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               <  1:2065410/93‑1 (MQ=255)

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

                                                           aCTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTAc                              >  1:941979/1‑93 (MQ=255)

                                                           aCTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTAc                              <  1:413334/93‑1 (MQ=255)

                                                           aCTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGGAc                              >  1:3978761/1‑93 (MQ=255)

                                                            cTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                             <  2:2915812/93‑1 (MQ=255)

                                                              ttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGt                           >  1:395322/1‑93 (MQ=255)

                                                              ttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGt                           >  2:5052461/1‑93 (MQ=255)

                                                                  cacgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaa                       <  1:1141021/93‑1 (MQ=255)

                                                                     gcagcaGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACa                    <  1:3776390/93‑1 (MQ=255)

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

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

                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCGGCTCCGGGTACAGTTCAAACAttt                 >  2:4245672/1‑93 (MQ=255)

                                                                         cagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtttt                <  1:3907899/93‑1 (MQ=255)

                                                                          agTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg               <  1:6757039/93‑1 (MQ=255)

                                                                          agTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg               <  2:4884607/93‑1 (MQ=255)

                                                                           gTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCACCGCGTACAGTTCAAACATTttgt              >  2:338293/1‑93 (MQ=255)

                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt             >  1:2759897/1‑93 (MQ=255)

                                                                             gCGATATCCGGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAGTTCACCTTCCGCACGCTGCTACGGGTACAGTTCAAACATTttgttg            <  1:4318570/93‑1 (MQ=255)

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

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

                                                                               gATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa          >  1:549205/1‑93 (MQ=255)

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

                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        <  2:4622790/93‑1 (MQ=255)

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

                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg         <  2:3823566/92‑1 (MQ=255)

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

                                                                                   tCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa      <  1:3348429/93‑1 (MQ=255)

                                                                                   tCCTGAGGGGTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa      >  2:5242839/1‑93 (MQ=255)

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

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

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

                                                                                    ccTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTAc     >  2:5968788/1‑93 (MQ=255)

                                                                                     cTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTCCAGTTCAAACATTTTGTTGAAGGTACg    >  2:61089/1‑93 (MQ=255)

                                                                                     cTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACg    >  2:6441012/1‑93 (MQ=255)

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

                                                                                      tGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGCAGGTACgc   >  1:3267326/1‑93 (MQ=255)

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

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

                                                                                       gaggGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCATGCTGCTCCGGGTACAGTTCAAACATTTTGTAGAAAGTACgca  <  2:7005183/93‑1 (MQ=255)

                                                                                       gaggGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgca  >  2:1524979/1‑93 (MQ=255)

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

                                                                                       gaggGTTGGTCAGGATGTTGTCGACTTCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTTTTGAAGGTACgca  <  2:2638762/93‑1 (MQ=255)

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837

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