BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	411,562	T→C	A76A (GCA→GCG)	araJ ←	AraJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	411,562	0	T	C	99.4%	969.3 / ‑4.1	309	A76A (GCA→GCG)	araJ	AraJ
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (164/143); total (165/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 320 total aligned reads displayed.

GCACAACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGTGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAATT  >  NC_000913/411472‑411650

  acaACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCat                                                                                      <  2:1126959/93‑1 (MQ=255)

   caACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCata                                                                                     <  1:1073142/93‑1 (MQ=255)

   caACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCata                                                                                     >  1:3213065/1‑93 (MQ=255)

   caACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCata                                                                                     <  2:2986164/93‑1 (MQ=255)

    aaCGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCCGCGCataa                                                                                    >  1:6748468/1‑93 (MQ=255)

     aCGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCataat                                                                                   >  1:10031798/1‑93 (MQ=255)

      cGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCataata                                                                                  <  2:9171911/93‑1 (MQ=255)

      cGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCataata                                                                                  <  1:4284485/93‑1 (MQ=255)

        ccaccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACg                                                                                <  2:2758254/93‑1 (MQ=255)

        ccaccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACg                                                                                >  1:7660266/1‑93 (MQ=255)

         caccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACga                                                                               >  1:2362958/1‑93 (MQ=255)

          accaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGGGGGATGATTGGCGCACCGACCACCACCCCCAGCGAATAATAGgag                                                                              >  1:2217798/1‑93 (MQ=255)

           ccaGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACgaga                                                                             <  2:437390/93‑1 (MQ=255)

             aGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATc                                                                           <  2:3080180/93‑1 (MQ=255)

             aGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATc                                                                           <  2:5373892/93‑1 (MQ=255)

             aGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATc                                                                           <  2:1283113/93‑1 (MQ=255)

             aGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGGGCGATGATTGGCGCACCGACCACCACCCCCAGCCCATAATACGAGATc                                                                           >  2:8280753/1‑93 (MQ=255)

              gAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCa                                                                          >  2:6588196/1‑93 (MQ=255)

               aaaCAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCat                                                                         <  1:5632213/93‑1 (MQ=255)

                aacaacAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCata                                                                        <  1:991782/93‑1 (MQ=255)

                aacaacAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCata                                                                        >  1:7477841/1‑93 (MQ=255)

                aacaacAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCata                                                                        <  2:7330351/93‑1 (MQ=255)

                aacaacAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCata                                                                        <  1:6022143/93‑1 (MQ=255)

                aacaacAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCata                                                                        <  1:8842770/93‑1 (MQ=255)

                 acaacaAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCatat                                                                       >  1:6883298/1‑93 (MQ=255)

                 acaacaAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCatat                                                                       >  1:4226679/1‑93 (MQ=255)

                 acaacaAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGGGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCatat                                                                       >  2:7438863/1‑93 (MQ=255)

                  caacaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATg                                                                      <  1:8184189/93‑1 (MQ=255)

                  caacaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATg                                                                      >  1:4926548/1‑93 (MQ=255)

                    acaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGcc                                                                    >  2:7996898/1‑93 (MQ=255)

                    acaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGcc                                                                    >  1:5139069/1‑93 (MQ=255)

                    acaaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGcc                                                                    <  1:3697696/93‑1 (MQ=255)

                     caaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGccc                                                                   <  2:6377008/93‑1 (MQ=255)

                     caaGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGccc                                                                   <  2:6279346/93‑1 (MQ=255)

                        gATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCcggc                                                                >  1:7288084/1‑93 (MQ=255)

                        gATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCcggc                                                                <  1:5545179/93‑1 (MQ=255)

                          tatGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCcggcgg                                                              >  2:235846/1‑93 (MQ=255)

                          tatGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCcggcgg                                                              >  1:4553527/1‑93 (MQ=255)

                            tGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCa                                                            >  1:9019046/1‑93 (MQ=255)

                               ttGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGa                                                         >  1:9727269/1‑93 (MQ=255)

                                tgagtgagTAGCGGCTGGAAAAGAGTGCGATGATTGGAGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCCGCTGCATGaa                                                        >  2:6497545/1‑93 (MQ=255)

                                 gagtgagtAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATCATACGAGATCATATGCCCGGCGGCACGAAt                                                       >  1:694570/1‑93 (MQ=255)

                                 gagtgagtAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAAt                                                       <  2:9586303/93‑1 (MQ=255)

                                  agtgagtAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATc                                                      <  2:9346633/93‑1 (MQ=255)

                                     gagtAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGaa                                                   <  1:166935/93‑1 (MQ=255)

                                      agtAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGaaa                                                  <  2:1929468/93‑1 (MQ=255)

                                         aGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTc                                               <  1:3852386/93‑1 (MQ=255)

                                          gCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTcc                                              <  2:7577325/93‑1 (MQ=255)

                                            ggCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTa                                            <  1:1843649/93‑1 (MQ=255)

                                            ggCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTa                                            >  2:6837958/1‑93 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGCATCGAAATTACTACGTTa                                       >  2:4716442/1‑93 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTGCGTTa                                       >  2:8745355/1‑93 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTa                                       <  1:1341069/93‑1 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTa                                       <  1:2151689/93‑1 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTa                                       <  1:2901853/93‑1 (MQ=255)

                                                 gAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTa                                       <  1:77459/93‑1 (MQ=255)

                                                  aaaaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTAt                                      >  1:6254808/1‑93 (MQ=255)

                                                  aaaaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTAt                                      >  1:4849120/1‑93 (MQ=255)

                                                   aaaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGACTCGAAATTCCTACGTTATg                                     >  1:7878405/1‑93 (MQ=255)

                                                   aaaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATg                                     >  2:6669683/1‑93 (MQ=255)

                                                    aaGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgc                                    >  2:6320270/1‑93 (MQ=255)

                                                      gagTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgcgc                                  >  1:4243462/1‑93 (MQ=255)

                                                      gagTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATgcgc                                  <  1:7173394/93‑1 (MQ=255)

                                                       agTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGcc                                 <  2:5354589/93‑1 (MQ=255)

                                                       agTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGcc                                 <  2:3579575/93‑1 (MQ=255)

                                                       agTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGcc                                 >  1:8379759/1‑93 (MQ=255)

                                                        gTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCa                                <  1:4652280/93‑1 (MQ=255)

                                                        gTGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCa                                <  1:1561457/93‑1 (MQ=255)

                                                         tGCGATGATTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAg                               >  1:556860/1‑93 (MQ=255)

                                                          gCGATGATTGACGCACAGACAACCACCACCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGCAATTCCTACGTTATGCGCCAGc                              >  2:8589895/1‑93 (MQ=255)

                                                            gatgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTc                            >  2:3382563/1‑93 (MQ=255)

                                                            gatgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTc                            >  1:8057949/1‑93 (MQ=255)

                                                             atgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTcc                           <  2:8077420/93‑1 (MQ=255)

                                                             atgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTcc                           <  2:9867687/93‑1 (MQ=255)

                                                              tgatTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCg                          <  2:5080315/93‑1 (MQ=255)

                                                                atTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTg                        <  1:7859191/93‑1 (MQ=255)

                                                                atTGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTg                        <  1:851631/93‑1 (MQ=255)

                                                                 ttGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGa                       >  1:5828674/1‑93 (MQ=255)

                                                                 ttGGCGCACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGa                       <  2:4889637/93‑1 (MQ=255)

                                                                    gcgcACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCa                    <  2:6879656/93‑1 (MQ=255)

                                                                     cgcACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCAc                   <  2:5616134/93‑1 (MQ=255)

                                                                       cACCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGc                 <  2:5576323/93‑1 (MQ=255)

                                                                        aCCGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGcc                >  1:9803408/1‑93 (MQ=255)

                                                                          cGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCa              <  2:4120793/93‑1 (MQ=255)

                                                                          cGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCa              <  1:3092936/93‑1 (MQ=255)

                                                                          cGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCa              <  1:9469903/93‑1 (MQ=255)

                                                                          cGACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCa              <  2:6020326/93‑1 (MQ=255)

                                                                           gACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCAt             >  2:8943930/1‑93 (MQ=255)

                                                                           gACCACCACCCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCAAt             >  1:7026282/1‑93 (MQ=255)

                                                                            accaccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATa            >  1:3469718/1‑93 (MQ=255)

                                                                             ccaccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATaa           <  2:6771445/93‑1 (MQ=255)

                                                                              caccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCcatcat          >  2:8153498/1‑93 (MQ=255)

                                                                              caccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAAt          <  1:4033743/93‑1 (MQ=255)

                                                                              caccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAAt          >  1:6552007/1‑93 (MQ=255)

                                                                               accaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATg         >  1:1971374/1‑93 (MQ=255)

                                                                               accaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATg         >  2:8390969/1‑93 (MQ=255)

                                                                                ccaccCCCAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGc        >  1:9636982/1‑93 (MQ=255)

                                                                                   cccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaa     >  2:552236/1‑93 (MQ=255)

                                                                                    ccccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCaaa    <  2:2407090/93‑1 (MQ=255)

                                                                                     cccAGCGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAAt   >  1:9301224/1‑93 (MQ=255)

GCACAACGCCACCAGAAACAACAAGATATGTTTGAGTGAGTAGCGGCTGGAAAAGAGTGCGATGATTGGCGCACCGACCACCACCCCCAGTGCATAATACGAGATCATATGCCCGGCGGCAGGAATCGAAATTCCTACGTTATGCGCCAGCTCCGTGAGCACGCCCATAATGCCAAATT  >  NC_000913/411472‑411650

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: `‑`