BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,328,807	A→G	Q243R (CAG→CGG)	proP →	proline permease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,328,807	0	A	G	100.0%	653.5 / NA	208	Q243R (CAG→CGG)	proP	proline permease II
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (102/106); total (102/106)

Only 100 of 215 total aligned reads displayed.

TCGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAAC  >  NC_000913/4328717‑4328895

 cGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGa                                                                                       >  1:3791459/1‑93 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGAt                                                                                      >  2:3750165/1‑93 (MQ=255)

  gggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGAt                                                                                      >  1:2014512/1‑93 (MQ=255)

   ggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATg                                                                                     >  2:708489/1‑93 (MQ=255)

   ggCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATg                                                                                     >  1:3386105/1‑93 (MQ=255)

      tttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGcc                                                                                  >  2:1496721/1‑93 (MQ=255)

      tttACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGcc                                                                                  >  1:643732/1‑93 (MQ=255)

        tACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCg                                                                                <  2:2013968/93‑1 (MQ=255)

        tACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGGTTGCGGGGTGGCCCg                                                                                >  2:41110/1‑93 (MQ=255)

        tACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGGTTGCGGGATGGCCCg                                                                                >  2:4621563/1‑93 (MQ=255)

         aCCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGa                                                                               <  2:1847427/93‑1 (MQ=255)

           cTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGaaa                                                                             <  2:2660555/93‑1 (MQ=255)

            tGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAg                                                                            <  1:1025445/93‑1 (MQ=255)

              cgcCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtc                                                                          <  1:2946524/93‑1 (MQ=255)

              cgcCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtc                                                                          <  1:3698774/93‑1 (MQ=255)

              cgcCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtc                                                                          <  1:4157165/93‑1 (MQ=255)

               gcCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtct                                                                         >  2:2822073/1‑93 (MQ=255)

               gcCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtct                                                                         >  1:1367847/1‑93 (MQ=255)

                ccATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                        >  2:1479077/1‑93 (MQ=255)

                ccATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGtctc                                                                        >  2:3672043/1‑93 (MQ=255)

                 cATGCGCTGGAAGAGACTCCGGCGTTGCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCg                                                                       <  1:2270149/93‑1 (MQ=255)

                 cATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCg                                                                       <  2:350443/93‑1 (MQ=255)

                  aTGCGCTGGAAGAGACTCCGGCGTTGCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGt                                                                      >  1:1947169/1‑93 (MQ=255)

                  aTGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGt                                                                      <  2:1702301/93‑1 (MQ=255)

                   tGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGtt                                                                     <  2:1122260/93‑1 (MQ=255)

                   tGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGtt                                                                     <  2:569118/93‑1 (MQ=255)

                   tGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGtt                                                                     <  2:507938/93‑1 (MQ=255)

                     cgcTGGAAGAGCCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTa                                                                   <  1:908237/93‑1 (MQ=255)

                     cgcTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTa                                                                   <  1:2397983/93‑1 (MQ=255)

                       cTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTaaa                                                                 >  2:3884838/1‑93 (MQ=255)

                        tGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGCTTGCGGGATGGCCCGAGAGTCTCGTTTAAag                                                                >  1:1853955/1‑93 (MQ=255)

                        tGGAAGAAACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAag                                                                >  2:4456985/1‑93 (MQ=255)

                         ggAAGAGACTCCTGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAaga                                                               <  1:773437/93‑1 (MQ=255)

                         ggAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAaga                                                               >  1:1132842/1‑93 (MQ=255)

                         ggAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAaga                                                               >  1:3007065/1‑93 (MQ=255)

                              agaCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGAttc                                                          >  1:3758504/1‑92 (MQ=255)

                                aCTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGcc                                                        <  1:2495477/93‑1 (MQ=255)

                                     ggCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaa                                                   >  2:2811081/1‑93 (MQ=255)

                                      gCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaaa                                                  >  2:3038150/1‑93 (MQ=255)

                                      gCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTaaa                                                  >  2:3496907/1‑93 (MQ=255)

                                       cGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAAt                                                 >  2:4461149/1‑93 (MQ=255)

                                       cGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAAt                                                 <  2:381138/93‑1 (MQ=255)

                                       cGTTCCAGCAGCATGACGATAAACTGAAGCAGGGCGTCCGTGAATTTTTGCGCGATGCCCCGAAAGCCTCGTTTAAAGAGATTGCCACTAAAt                                                 <  2:308065/93‑1 (MQ=255)

                                        gTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATa                                                >  1:1529935/1‑93 (MQ=255)

                                        gTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATa                                                >  1:3521990/1‑93 (MQ=255)

                                        gTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGGTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATa                                                >  2:93044/1‑93 (MQ=255)

                                         ttCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATAc                                               >  1:3236208/1‑93 (MQ=255)

                                          tCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACt                                              <  2:122549/93‑1 (MQ=255)

                                           ccAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTg                                             <  2:3541758/93‑1 (MQ=255)

                                           ccAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTg                                             <  1:4093251/93‑1 (MQ=255)

                                            cagcagCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTgg                                            <  2:108581/93‑1 (MQ=255)

                                            cagcagCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTgg                                            <  2:831792/93‑1 (MQ=255)

                                             agcagcATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgc                                           >  2:1264955/1‑93 (MQ=255)

                                              gcagcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgcg                                          >  2:874918/1‑93 (MQ=255)

                                              gcagcaTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGgcg                                          >  1:4523309/1‑93 (MQ=255)

                                                   aTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGGTGGCCCGAGAGTCTCGTTTAAAGAGATTGACACTAAATACTGGCGCAGcc                                     >  2:1154114/1‑93 (MQ=255)

                                                   aTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGcc                                     <  2:1374765/93‑1 (MQ=255)

                                                   aTGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGcc                                     <  2:3185560/93‑1 (MQ=255)

                                                      tCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGt                                  >  2:3580934/1‑93 (MQ=255)

                                                      tCGATAAACTGGAACAGGGCGACCGGGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCATCCTGt                                  <  2:1944184/93‑1 (MQ=255)

                                                       cGATAGACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGtt                                 <  2:753444/93‑1 (MQ=255)

                                                       cGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGtt                                 <  2:3275237/93‑1 (MQ=255)

                                                          tAAACTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGAc                              >  1:2725936/1‑93 (MQ=255)

                                                            aaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACAt                            <  1:2526406/93‑1 (MQ=255)

                                                            aaCTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGGCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACAt                            >  1:1650914/1‑93 (MQ=255)

                                                              cTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGt                          <  1:4661680/93‑1 (MQ=255)

                                                              cTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGt                          <  1:3311886/93‑1 (MQ=255)

                                                              cTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGt                          <  1:2605534/93‑1 (MQ=255)

                                                              cTGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGt                          <  2:232981/93‑1 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTa                         >  2:2501627/1‑93 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTa                         <  1:348569/93‑1 (MQ=255)

                                                               tGGAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTa                         <  2:4654262/93‑1 (MQ=255)

                                                                ggAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTAt                        >  2:316485/1‑93 (MQ=255)

                                                                ggAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTAt                        >  2:2440067/1‑93 (MQ=255)

                                                                 gAACAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTAtt                       <  1:4185659/93‑1 (MQ=255)

                                                                  aaCAGGGCGACCGTGAGGGTTTGCGGGATGACCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTg                      <  2:2972308/93‑1 (MQ=255)

                                                                   aCAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTgg                     >  1:3397210/1‑93 (MQ=255)

                                                                    cAGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGt                    <  1:3996072/93‑1 (MQ=255)

                                                                     aGGGCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTc                   >  1:1289637/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                  >  1:3941227/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                  >  2:2905891/1‑93 (MQ=255)

                                                                      gggCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCt                  <  2:2487629/93‑1 (MQ=255)

                                                                        gCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTgg                <  1:1043740/93‑1 (MQ=255)

                                                                        gCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTgg                <  1:3127876/93‑1 (MQ=255)

                                                                        gCGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTgg                <  2:1915569/93‑1 (MQ=255)

                                                                         cGACCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGt               <  2:3520180/93‑1 (MQ=255)

                                                                           aCCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTaa             <  1:264185/93‑1 (MQ=255)

                                                                           aCCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTaa             >  2:4493332/1‑93 (MQ=255)

                                                                           aCCGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTaa             <  1:1927000/93‑1 (MQ=255)

                                                                             cGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAAtt           >  2:3834011/1‑93 (MQ=255)

                                                                             cGTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAAtt           >  1:2225594/1‑93 (MQ=255)

                                                                              gTGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTg          >  2:4250565/1‑93 (MQ=255)

                                                                               tGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGc         >  1:737703/1‑93 (MQ=255)

                                                                               tGAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGc         >  1:995192/1‑93 (MQ=255)

                                                                                gAAGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGcc        >  1:3814231/1‑93 (MQ=255)

                                                                                 aaGGTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGcca       <  2:1857566/93‑1 (MQ=255)

                                                                                   ggTTTGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGccacc     <  1:3170988/93‑1 (MQ=255)

                                                                                     tttGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCcaacaa   >  2:649152/1‑93 (MQ=255)

                                                                                     tttGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCaa   <  1:4639363/93‑1 (MQ=255)

                                                                                      ttGCGGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAAc  <  2:3803509/93‑1 (MQ=255)

TCGGGCTTTACCTGCGCCATGCGCTGGAAGAGACTCCGGCGTTCCAGCAGCATGTCGATAAACTGGAACAGGGCGACCGTGAAGGTTTGCAGGATGGCCCGAAAGTCTCGTTTAAAGAGATTGCCACTAAATACTGGCGCAGCCTGTTGACATGTATTGGTCTGGTAATTGCCACCAAC  >  NC_000913/4328717‑4328895

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