BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,219,434	A→G	A211A (GCT→GCC)	arp ←	Arp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,219,434	0	A	G	99.2%	767.6 / ‑4.2	236	A211A (GCT→GCC)	arp	Arp
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (125/109); total (126/110)
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 263 total aligned reads displayed.

CACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGC  >  NC_000913/4219344‑4219523

 accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGgcggca                                                                                        >  1:804381/1‑93 (MQ=255)

 accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGgcggca                                                                                        >  1:1029489/1‑93 (MQ=255)

  ccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaa                                                                                       >  2:2710305/1‑93 (MQ=255)

  ccaGTAATATGGCTGCCATAAAACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCCgcaa                                                                                       >  2:803057/1‑93 (MQ=255)

   caGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaac                                                                                      >  1:733676/1‑93 (MQ=255)

   caGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaac                                                                                      >  1:2717241/1‑93 (MQ=255)

    aGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaacg                                                                                     >  2:3471070/1‑93 (MQ=255)

       aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaacgcaa                                                                                  <  2:1778584/93‑1 (MQ=255)

       aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaacgcaa                                                                                  <  1:1506454/93‑1 (MQ=255)

       aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGgcaacgcaa                                                                                  <  2:339961/93‑1 (MQ=255)

         tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGcaacgcaacg                                                                                >  1:3052905/1‑93 (MQ=255)

         tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGcaacgcaacg                                                                                >  2:3372423/1‑93 (MQ=255)

          atGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCaacgcaacgc                                                                               <  1:3353295/93‑1 (MQ=255)

             gCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGcccc                                                                            <  1:3191077/93‑1 (MQ=255)

              cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCt                                                                           <  1:2149509/93‑1 (MQ=255)

                  cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAgg                                                                       >  1:2326782/1‑93 (MQ=255)

                  cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAgg                                                                       >  1:3124534/1‑93 (MQ=255)

                   aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGa                                                                      <  1:2414099/93‑1 (MQ=255)

                    tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGAc                                                                     >  1:2467246/1‑93 (MQ=255)

                    tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGAc                                                                     >  1:36572/1‑93 (MQ=255)

                    tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGAc                                                                     <  1:4174215/93‑1 (MQ=255)

                     aaCACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACg                                                                    <  1:4291214/93‑1 (MQ=255)

                      acacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACgg                                                                   >  1:2582657/1‑93 (MQ=255)

                       cacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGc                                                                  <  2:1542106/93‑1 (MQ=255)

                        acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCt                                                                 <  1:2283633/93‑1 (MQ=255)

                        acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCt                                                                 <  2:3549764/93‑1 (MQ=255)

                        acTGATGCCATGTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCGCCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCt                                                                 <  1:109023/93‑1 (MQ=255)

                         cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTg                                                                <  1:1280372/93‑1 (MQ=255)

                         cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTg                                                                >  2:588770/1‑93 (MQ=255)

                         cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTg                                                                >  1:4410230/1‑93 (MQ=255)

                          tGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTgg                                                               >  1:2945391/1‑93 (MQ=255)

                           gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGa                                                              >  1:4371897/1‑93 (MQ=255)

                            aTGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGaa                                                             >  2:4289895/1‑93 (MQ=255)

                               ccATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAAcc                                                          <  1:970217/93‑1 (MQ=255)

                                 aTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCaa                                                        <  1:2200304/93‑1 (MQ=255)

                                  tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGACACCAAc                                                       >  1:3874873/1‑93 (MQ=255)

                                     gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGcc                                                    <  2:1890480/93‑1 (MQ=255)

                                     gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGcc                                                    <  2:243037/93‑1 (MQ=255)

                                      aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCa                                                   >  2:2992952/1‑93 (MQ=255)

                                       ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaa                                                  <  2:1144850/93‑1 (MQ=255)

                                       ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaa                                                  <  2:2526639/93‑1 (MQ=255)

                                        tacATCTTGTCAAAATATGCTTTCATACCTAGCCCGAAATCCCCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaaa                                                 <  1:1934233/91‑1 (MQ=255)

                                        tCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaaa                                                 <  1:2823764/93‑1 (MQ=255)

                                          cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGt                                               >  1:3944556/1‑93 (MQ=255)

                                          cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGt                                               <  1:4733255/93‑1 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAgac                                              >  2:1478601/1‑91 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                              <  2:2542837/93‑1 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                              <  1:2029223/93‑1 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                              <  1:1199671/93‑1 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                              <  1:3521519/93‑1 (MQ=255)

                                           aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                              <  1:4042649/93‑1 (MQ=255)

                                            tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             <  2:3788031/93‑1 (MQ=255)

                                            tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             <  2:4363353/93‑1 (MQ=255)

                                            tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             >  1:566886/1‑93 (MQ=255)

                                            tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                             >  1:3501465/1‑93 (MQ=255)

                                             cTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTg                                            >  1:1424065/1‑93 (MQ=255)

                                              ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTgg                                           <  2:704522/93‑1 (MQ=255)

                                              ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTgg                                           >  2:948954/1‑93 (MQ=255)

                                              ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGACCCTCAGGACGGCTGGAAATCAACGCCAAAGGCTgg                                           >  2:412971/1‑93 (MQ=255)

                                               tGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGt                                          <  2:4481567/93‑1 (MQ=255)

                                                 tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttt                                        <  1:3977957/93‑1 (MQ=255)

                                                   aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                      >  2:3965469/1‑93 (MQ=255)

                                                     aaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGc                                    >  2:3192146/1‑93 (MQ=255)

                                                     aaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGc                                    <  1:2260906/93‑1 (MQ=255)

                                                      atatGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCa                                   >  1:1792993/1‑93 (MQ=255)

                                                         tGCTTTCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttt                                >  1:1888158/1‑93 (MQ=255)

                                                            tttCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttttt                             >  1:488168/1‑93 (MQ=255)

                                                            tttCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttttt                             >  1:4659658/1‑93 (MQ=255)

                                                            tttCATAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttttt                             <  1:4139000/93‑1 (MQ=255)

                                                                aTAACTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgttg                         >  2:3944886/1‑93 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCTCCTCAGGACGGCTGGAAACCATCGCCAAAGTCTGGTTTTTGCATTCATTGTTGaa                       >  2:4214490/1‑93 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                       >  1:3146023/1‑93 (MQ=255)

                                                                  aaCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                       >  1:391061/1‑93 (MQ=255)

                                                                   aCTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                      >  1:3621369/1‑93 (MQ=255)

                                                                    cTAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGt                     >  1:4352529/1‑93 (MQ=255)

                                                                     tAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAGCCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTg                    >  1:2667666/1‑93 (MQ=255)

                                                                     tAGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTg                    >  2:1664238/1‑93 (MQ=255)

                                                                      aGCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGc                   <  2:2373482/93‑1 (MQ=255)

                                                                       gCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                  >  2:2838584/1‑93 (MQ=255)

                                                                       gCCCGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                  >  2:2956170/1‑93 (MQ=255)

                                                                         ccGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTTCGTg                >  2:3418680/1‑93 (MQ=255)

                                                                          cGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga               <  2:308838/93‑1 (MQ=255)

                                                                          cGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga               <  1:1446987/93‑1 (MQ=255)

                                                                          cGAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga               >  2:1333595/1‑93 (MQ=255)

                                                                           gAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTTCGTgag              <  2:1456011/93‑1 (MQ=255)

                                                                           gAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgag              >  1:3383047/1‑93 (MQ=255)

                                                                           gAAATCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgag              >  2:2776546/1‑93 (MQ=255)

                                                                            aaaTCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGATTTTGCATTTTTTGTTGAAGTGCGTgaga             >  1:3723888/1‑93 (MQ=255)

                                                                             aaTCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGaa            >  1:3313656/1‑93 (MQ=255)

                                                                             aaTCACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGaa            >  2:4670774/1‑93 (MQ=255)

                                                                                cACCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtt         <  1:3888269/93‑1 (MQ=255)

                                                                                 aCCGCAGGCGGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt        <  2:877594/93‑1 (MQ=255)

                                                                                    gcaggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc     <  2:3216268/93‑1 (MQ=255)

                                                                                    gcaggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc     <  1:4235287/93‑1 (MQ=255)

                                                                                     caggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCa    <  2:2029507/93‑1 (MQ=255)

                                                                                     caggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCa    <  2:1196667/93‑1 (MQ=255)

                                                                                      aggcggCAACGCAAGGGCCCTCAGGACGGCTGGAAACCAAGGCCAAAGTCTGGTTTTTGTATTTTTTGTTGAAGTGCGTGAGAAATTTGTCAg   >  2:4534895/1‑93 (MQ=255)

                                                                                      aggcggCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAg   >  2:1857094/1‑93 (MQ=255)

                                                                                       ggcggcAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGc  <  2:2363656/93‑1 (MQ=255)

                                                                                       ggcggcAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGc  <  2:1029489/93‑1 (MQ=255)

CACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGC  >  NC_000913/4219344‑4219523

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