BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,847,298	A→G	K146K (AAA→AAG)	sppA →	protease IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,847,298	0	A	G	99.1%	323.5 / NA	106	K146K (AAA→AAG)	sppA	protease IV
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (48/57); minor base T (0/1); total (48/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387

tgtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGcc                                                                                         <  1:772100/93‑1 (MQ=255)

 gtgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCg                                                                                        >  2:595710/1‑93 (MQ=255)

  tgATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCgg                                                                                       >  2:486103/1‑93 (MQ=255)

   gATGGATCTGAAAAACTTCGCAGGCGGCGATCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCATGACAGCGGGAATCCGGt                                                                                      <  2:156293/93‑1 (MQ=255)

   gATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGt                                                                                      >  2:4428743/1‑93 (MQ=255)

   gATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGt                                                                                      >  1:2799056/1‑93 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  1:417471/93‑1 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     >  1:118999/1‑93 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  1:4298971/93‑1 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  2:4331001/93‑1 (MQ=255)

    aTGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGtt                                                                                     <  1:369140/93‑1 (MQ=255)

     tGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGttt                                                                                    <  2:4366124/93‑1 (MQ=255)

      taaTCTGAAAAACTTCGCAGGCGGCGACCAACCATCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  1:813832/91‑1 (MQ=255)

      ggATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTa                                                                                   <  2:2638822/93‑1 (MQ=255)

       gATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTAt                                                                                  <  1:3880286/93‑1 (MQ=255)

       gATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTAt                                                                                  <  1:232110/93‑1 (MQ=255)

         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                >  2:752088/1‑93 (MQ=255)

         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                >  1:1567957/1‑93 (MQ=255)

         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                <  2:951913/93‑1 (MQ=255)

         tCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGATTCGTGACAGCGGGAAGCCGGTTTATGc                                                                                >  1:1913474/1‑93 (MQ=255)

          cTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGGAGCCGGTTTATGcc                                                                               >  2:222649/1‑93 (MQ=255)

          cTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGcc                                                                               <  1:1843318/93‑1 (MQ=255)

          cTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGcc                                                                               <  1:4621045/93‑1 (MQ=255)

            gAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGt                                                                             >  2:4542169/1‑93 (MQ=255)

             aaaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGtt                                                                            <  2:448515/93‑1 (MQ=255)

             aaaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGtt                                                                            >  2:1909962/1‑93 (MQ=255)

              aaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTg                                                                           <  1:1111648/93‑1 (MQ=255)

              aaaaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTg                                                                           <  2:4527447/93‑1 (MQ=255)

               aaaCCTCGCAGGCGGCGGCCACCCGTCTGTGCAGTACATCGGCAAAGCTCTGAAGGAGTTTCGTGACAGCGGGAAGGCGGTTTATGCCGTTgg                                                                          <  2:4022872/93‑1 (MQ=255)

                aaCTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGc                                                                         >  1:588608/1‑93 (MQ=255)

                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCga                                                                       <  2:775013/93‑1 (MQ=255)

                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCga                                                                       <  1:14244/93‑1 (MQ=255)

                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCga                                                                       <  1:2929715/93‑1 (MQ=255)

                  cTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCga                                                                       <  2:769372/93‑1 (MQ=255)

                    tCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCgaga                                                                     <  1:3818507/93‑1 (MQ=255)

                     cGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGaa                                                                    <  2:954175/93‑1 (MQ=255)

                      gCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAAc                                                                   >  1:1820345/1‑93 (MQ=255)

                         ggcggcGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTAc                                                                <  2:2449962/93‑1 (MQ=255)

                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  2:2524908/93‑1 (MQ=255)

                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  1:2818056/93‑1 (MQ=255)

                          gcggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACa                                                               <  1:3101290/93‑1 (MQ=255)

                            ggcgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGc                                                             <  1:3636895/93‑1 (MQ=255)

                              cgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCa                                                           <  1:3857076/93‑1 (MQ=255)

                              cgACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCa                                                           <  2:3819441/93‑1 (MQ=255)

                                aCCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAgg                                                         <  2:3455026/93‑1 (MQ=255)

                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        >  1:2120593/1‑93 (MQ=255)

                                 ccAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGCTTATGCCGTTGGCGAGAACTACAGCCAggg                                                        >  1:2120592/1‑93 (MQ=255)

                                   aaCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGc                                                      >  1:3660390/1‑93 (MQ=255)

                                    aCCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCa                                                     >  1:3500405/1‑93 (MQ=255)

                                     ccGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCaa                                                    <  1:3174543/93‑1 (MQ=255)

                                       gTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAata                                                  <  2:3334867/93‑1 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:594100/1‑93 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  1:4174108/1‑93 (MQ=255)

                                        tCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAatat                                                 >  2:2632144/1‑93 (MQ=255)

                                            tGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtc                                             >  2:180185/1‑93 (MQ=255)

                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            >  2:154708/1‑93 (MQ=255)

                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            >  1:383623/1‑93 (MQ=255)

                                             gCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtct                                            >  2:960611/1‑93 (MQ=255)

                                              cAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtctc                                           <  2:1079331/93‑1 (MQ=255)

                                              cAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTAtctc                                           >  1:3730781/1‑93 (MQ=255)

                                               aGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCg                                          <  1:1487143/93‑1 (MQ=255)

                                                gTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         >  1:583902/1‑93 (MQ=255)

                                                atACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGc                                         <  1:413838/92‑1 (MQ=255)

                                                  aCATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCa                                       >  1:1173255/1‑93 (MQ=255)

                                                    aTCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGt                                     >  1:1483010/1‑93 (MQ=255)

                                                    aTCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGt                                     <  1:2016984/93‑1 (MQ=255)

                                                    aTCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGt                                     >  2:3012433/1‑93 (MQ=255)

                                                    aTCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGt                                     >  2:4333822/1‑93 (MQ=255)

                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    >  1:771277/1‑93 (MQ=255)

                                                     tCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGtt                                    >  1:2837176/1‑93 (MQ=255)

                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGttt                                   >  1:2837135/1‑93 (MQ=255)

                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGttt                                   <  1:870036/93‑1 (MQ=255)

                                                      cGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGttt                                   <  2:339005/93‑1 (MQ=255)

                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  >  1:2202490/1‑93 (MQ=255)

                                                       ggCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTc                                  <  2:2815241/93‑1 (MQ=255)

                                                         cAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGc                                <  2:335087/93‑1 (MQ=255)

                                                          aaaGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGcc                               <  2:1736477/93‑1 (MQ=255)

                                                           aaGCTCTGAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCa                              >  1:1281588/1‑93 (MQ=255)

                                                               tctGAAAGAGTTGCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATaa                          <  2:4402258/93‑1 (MQ=255)

                                                                  gAAAGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAt                       <  1:2511709/93‑1 (MQ=255)

                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     >  1:1505458/1‑93 (MQ=255)

                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     >  1:632816/1‑93 (MQ=255)

                                                                    aaGAGTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAAttt                     >  1:3139796/1‑93 (MQ=255)

                                                                      gagTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTgg                   >  1:468767/1‑93 (MQ=255)

                                                                       agTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGc                  >  2:4554488/1‑93 (MQ=255)

                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 <  2:3580719/93‑1 (MQ=255)

                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 >  1:3720851/1‑93 (MQ=255)

                                                                        gTTTCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCt                 >  2:2445522/1‑93 (MQ=255)

                                                                         tttCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTg                >  2:4279691/1‑93 (MQ=255)

                                                                           tCGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtc              >  2:2479292/1‑93 (MQ=255)

                                                                            cGTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtct             >  2:3950771/1‑93 (MQ=255)

                                                                             gTGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGtctc            <  2:2668589/93‑1 (MQ=255)

                                                                              tGACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTcc           <  2:907984/93‑1 (MQ=255)

                                                                               gACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCg          >  1:1534278/1‑93 (MQ=255)

                                                                               gACAGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCg          >  2:921945/1‑93 (MQ=255)

                                                                                  aGCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCaa       <  2:2672253/93‑1 (MQ=255)

                                                                                  aGCGGGAAGCCGGTTTATGCCGTCGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCaa       <  1:3571283/93‑1 (MQ=255)

                                                                                   gCGGGAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAg      <  1:2103052/93‑1 (MQ=255)

                                                                                      ggAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCg   <  1:2788/93‑1 (MQ=255)

                                                                                       gAAGCCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGt  <  2:3024485/93‑1 (MQ=255)

TGTGATGGATCTGAAAAACTTCGCAGGCGGCGACCAACCGTCTATGCAGTACATCGGCAAAGCTCTGAAAGAGTTTCGTGACAGCGGGAAACCGGTTTATGCCGTTGGCGAGAACTACAGCCAGGGGCAATATTATCTCGCCAGTTTCGCCAATAAAATTTGGCTGTCTCCGCAAGGCGT  >  NC_000913/1847208‑1847387

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