BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,055,398	A→G	S73P (TCC→CCC)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,055,398	0	A	G	99.0%	727.6 / ‑3.3	209	S73P (TCC→CCC)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (90/117); total (90/119)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

Only 100 of 223 total aligned reads displayed.

GTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGA  >  NC_000913/4055309‑4055487

gTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTTGAGTCGGCGCAGAACGGGATAAACACTGCggtgggtgg                                                                                        >  1:3785719/1‑92 (MQ=255)

gTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGc                                                                                        <  1:2896928/93‑1 (MQ=255)

gTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGAggtgggtgg                                                                                        >  1:1719886/1‑92 (MQ=255)

   aTAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTc                                                                                     <  1:3559578/93‑1 (MQ=255)

   aTAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTc                                                                                     <  2:2943916/93‑1 (MQ=255)

    tAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCt                                                                                    <  2:4846519/93‑1 (MQ=255)

     aGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTg                                                                                   >  2:569590/1‑93 (MQ=255)

     aGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTg                                                                                   >  2:3813452/1‑93 (MQ=255)

      gCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTgg                                                                                  >  1:1152786/1‑93 (MQ=255)

       ccTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGc                                                                                 <  1:3264177/93‑1 (MQ=255)

          tGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATc                                                                              <  2:3965856/93‑1 (MQ=255)

          tGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATc                                                                              <  2:2210629/93‑1 (MQ=255)

           gCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACGGCGGTGGGTGCGTCTGGCATCa                                                                             <  2:118731/93‑1 (MQ=255)

            cAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAg                                                                            <  2:4185212/93‑1 (MQ=255)

            cAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAg                                                                            >  2:857784/1‑93 (MQ=255)

               ggTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCAc                                                                         >  1:6065922/1‑93 (MQ=255)

               ggTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCAc                                                                         <  2:5066283/93‑1 (MQ=255)

                 tGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCa                                                                       <  2:1523769/93‑1 (MQ=255)

                   ccAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATg                                                                     <  2:2149443/93‑1 (MQ=255)

                    cAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGt                                                                    >  2:1620502/1‑93 (MQ=255)

                    cAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGt                                                                    >  2:2619733/1‑93 (MQ=255)

                    cAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGt                                                                    <  2:1158117/93‑1 (MQ=255)

                    cAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGt                                                                    <  2:3430303/93‑1 (MQ=255)

                      ggTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCg                                                                  >  1:517457/1‑93 (MQ=255)

                        ttCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCggc                                                                >  2:5531086/1‑92 (MQ=255)

                        ttCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGa                                                                <  1:1702171/93‑1 (MQ=255)

                         tCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAc                                                               <  1:1881958/93‑1 (MQ=255)

                         tCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAc                                                               <  1:2030566/93‑1 (MQ=255)

                         tCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAc                                                               >  1:2730833/1‑93 (MQ=255)

                         tCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAc                                                               <  1:4451379/93‑1 (MQ=255)

                          cAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACt                                                              >  2:4928907/1‑93 (MQ=255)

                          cAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACt                                                              <  1:4389775/93‑1 (MQ=255)

                          cAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGCCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGAca                                                              >  2:4479825/1‑92 (MQ=255)

                           aaGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGACGAACGGGTCAATCCCTGCGGTGGGTGCGTCTGGCACCAGCACCATGTCGGACTc                                                             >  2:794077/1‑93 (MQ=255)

                           aaGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGACAATCACTGCGGTGGGTGCGTATGGCATCAGCACCATGTCGGACTc                                                             >  1:621800/1‑93 (MQ=255)

                            aGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCg                                                            <  2:1775211/93‑1 (MQ=255)

                             ggATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGt                                                           <  2:2976892/93‑1 (MQ=255)

                                tGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTaa                                                        <  2:189813/93‑1 (MQ=255)

                                tGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTaa                                                        <  2:5221072/93‑1 (MQ=255)

                                tGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTaa                                                        >  1:4156141/1‑93 (MQ=255)

                                 gTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAAt                                                       <  1:3167680/93‑1 (MQ=255)

                                   cGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGc                                                     <  2:323580/93‑1 (MQ=255)

                                    gCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGTATCAGCACCATGTCGGACTCGTTAATGct                                                    >  2:309442/1‑92 (MQ=255)

                                    gCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGcc                                                    >  1:527872/1‑93 (MQ=255)

                                    gCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGcc                                                    >  2:574271/1‑93 (MQ=255)

                                      aaCGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCtt                                                  >  1:1054836/1‑93 (MQ=255)

                                       aCGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCttt                                                 <  1:5223176/93‑1 (MQ=255)

                                       aCGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCttt                                                 <  1:417705/93‑1 (MQ=255)

                                       aCGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCttt                                                 >  1:2591152/1‑93 (MQ=255)

                                         ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTcc                                               <  2:3770775/93‑1 (MQ=255)

                                         ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTcc                                               <  1:1377151/93‑1 (MQ=255)

                                         ggATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTCTcc                                               <  2:3520951/93‑1 (MQ=255)

                                          gATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCa                                              <  2:1290935/93‑1 (MQ=255)

                                            tAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgc                                            <  2:3587828/93‑1 (MQ=255)

                                            tAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgc                                            <  1:1106739/93‑1 (MQ=255)

                                            tAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgc                                            <  1:2712963/93‑1 (MQ=255)

                                             aaTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgcc                                           >  1:5910337/1‑93 (MQ=255)

                                              aTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccg                                          >  1:2763909/1‑93 (MQ=255)

                                              aTCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccg                                          >  2:3869888/1‑93 (MQ=255)

                                               tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgc                                         >  1:3455487/1‑93 (MQ=255)

                                               tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgc                                         >  2:462226/1‑93 (MQ=255)

                                               tCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAgccgc                                         >  1:2272980/1‑93 (MQ=255)

                                                   ggTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAAt                                     >  2:2379430/1‑93 (MQ=255)

                                                   ggTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAAt                                     <  2:1546726/93‑1 (MQ=255)

                                                   ggTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAAt                                     <  2:3367646/93‑1 (MQ=255)

                                                    ttGGAGTCGGCGAAGAACGGGTCAATCACTGCTGTGGGTGCGTCTGGTATCAGCACCATTTCGGACTCGTTAATGTTTTTCCAGCCGCCAATc                                    >  2:5453710/2‑93 (MQ=255)

                                                    gTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATc                                    >  2:4000571/1‑93 (MQ=255)

                                                     tGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCg                                   >  2:5314580/1‑93 (MQ=255)

                                                      ggAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCga                                  <  1:568429/93‑1 (MQ=255)

                                                      ggAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCga                                  <  1:4840643/93‑1 (MQ=255)

                                                       gAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgag                                 <  1:609813/93‑1 (MQ=255)

                                                         gTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgagga                               <  1:839724/93‑1 (MQ=255)

                                                         gTCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgagga                               >  2:4358236/1‑93 (MQ=255)

                                                          tCGGCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCgaggag                              >  2:561916/1‑93 (MQ=255)

                                                            ggCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGcc                            <  2:2082079/93‑1 (MQ=255)

                                                             gCGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCg                           <  1:5743936/93‑1 (MQ=255)

                                                              cGAAGAACGGGTCAATCACTGCGGTGGGTGCGTCTGGCTTCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGt                          >  2:2011669/1‑93 (MQ=255)

                                                               gaagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTc                         >  1:4582907/1‑93 (MQ=255)

                                                                aagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCa                        >  1:2334359/1‑93 (MQ=255)

                                                                aagaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCa                        >  2:4067515/1‑93 (MQ=255)

                                                                  gaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCaaa                      >  1:4270695/1‑93 (MQ=255)

                                                                  gaaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCaaa                      <  1:5296239/93‑1 (MQ=255)

                                                                   aaCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAAc                     <  1:130866/93‑1 (MQ=255)

                                                                    aCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACa                    <  2:2498895/93‑1 (MQ=255)

                                                                    aCGGGTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACa                    <  1:2960635/93‑1 (MQ=255)

                                                                       ggTCAATCCCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGCCTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAttt                 <  2:2743396/93‑1 (MQ=255)

                                                                       ggTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAttt                 <  2:1438796/93‑1 (MQ=255)

                                                                        gTCAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACAtttt                >  2:2534590/1‑93 (MQ=255)

                                                                         gcAATCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAGTGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTg               <  1:5018532/92‑1 (MQ=255)

                                                                           aaTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGcc             >  2:1016542/1‑93 (MQ=255)

                                                                            aTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCct            <  1:550377/93‑1 (MQ=255)

                                                                            aTCACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCct            >  1:5855882/1‑93 (MQ=255)

                                                                              cACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttc          <  2:4007815/93‑1 (MQ=255)

                                                                              cACTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttc          <  2:5931151/93‑1 (MQ=255)

                                                                               aCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttct         >  1:5559355/1‑93 (MQ=255)

                                                                               aCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttct         >  2:306527/1‑93 (MQ=255)

                                                                               aCTGCGGTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCcttct         <  2:677068/93‑1 (MQ=255)

                                                                                    ggtgggtgCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa    <  1:4623879/93‑1 (MQ=255)

                                                                                    ggtgggtgCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaa    <  2:1743610/93‑1 (MQ=255)

                                                                                     gTGGGTGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCgaag   >  2:1822894/1‑93 (MQ=255)

GTCATAGCCTTGCAGGGTGCCAGGTTCAAGGATGTCGCAACGGATAATCAGGGTGGAGTCGGCGAAGAACGGGTCAATCACTGCGGTGGATGCGTCTGGCATCAGCACCATGTCGGACTCGTTAATGCCTTTCCAGCCGCCAATCGAGGAGCCGTCAAACATTTTGCCTTCTTCGAAGA  >  NC_000913/4055309‑4055487

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